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METHOD AND SYSTEM OF REGION-BASED IMAGE CODING WITH 

Dynamic Streaming of Code Blocks 



5 Field of the Invention 

The present invention relates generally to image coding, and more panicularly to the 
compression and streaming of scalable and content-based, randomly accessible digital 
still images. 

10 Background of the Invention 

With the rapid growth of the internet and multimedia applications, there is a great 
demand for new image coding tools that that will provide for high quality processing 
capability, an efficient internal architecture, and flexibility in terms of future 
technological advances. This is a challenge that has been put forth before the JPEG 
15 2000 Committee. Although the main focus should be to provide state-of-the-art 
compression performance, JPEG 2000 should also offer unprecedented content based 
accessibility of the compressed format to support applications of various needs. It is 
highly preferred and advantageous that image content be accessed, manipulated, 
exchanged, and stored in compact form. 

20 In order for JPEG 2000 to be the standard coding foundation of new generation image 
processing systems, it must provide for efficiency in coding, different types of still 
images (bi-level, gray-level, color) with different characteristics (natural images, 
scientific, medical imagery, rendered graphics, text, etc.) within a unified system. In 
addition to providing low bit rate operation with quality performance superiority to 

25 existing standards, this new system should include many modem features as listed in 
the JPEG 2000 requirement document. 

The open architecture and the set of algorithms presented in this document are based 
on Digital Accelerator Corporation's (DAC) Region based Image Coding System 
(RICS). RICS has not only achieved a rate distortion performance competitive with 
30 best known compression techniques, but also demonstrates a high degree of openness 
and flexibility that will accommodate most well known algorithms as well as new yet 
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to be implemented. The features supponed by RICS covers almost all those listed in 
JPEG 2000 Requirements document. 



Generality 

Generality is a primary concern in the architectural design of RICS. The RICS 
5 architecture attempts to be an integrated platform that supports and facilitates a 
variety of applications that may have different characteristics and requirements. For 
example, providing efficient lossless and lossy compression in a single code stream; 
efficient processing of compound documents containing both bi-level (text) and color 
imagery; progressive transmission by pixel accuracy and/or by resolution; random 
10 access of arbitrary shaped regions; and so on. 

Openness 

We also understand that the new JPEG 2000 standard is intended to be a dynamic, 
rather than static, suite of coding tools that support the new generation imagery 
applications and, at the same time, keeping abreast with the progress of technology. 
15 The mathematical foundation for those leading candidates of new image coding 
methods (most noticeably the various types of multi-resolution analysis techniques, 
such as wavelet transforms) is relatively young and still under intensive investigation. 
If an architectural design is based on or restricted to one or several particular existing 
coding methods, it may become outdated very quickly. 

20 In attempting to produce a flexible and open platform, the RICS architecture 
organizes the image coding process into a set of functionally separable modules, such 
that each module can be developed and optimized individually. Furthermore, all 
modules in the system are designed to be functionally orthogonal with each other, in 
the sense that a new algorithm, without affecting the functionality of other modules, 

25 can effectively replace the base algorithm of any specific module. 

This open architecture will not only be able to accommodate future new algorithms, 
but also makes compatibility with other standards an easy and natural extension of 
many concepts, as will be explained later in this document. 
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Accessibility 

Content based accessibility is becoming an imponant feature in supporting 
applications such as multimedia database query, internet server-client interaction, 
content production, remote diagnostics, and interactive entertainment. The content- 
5 based accessibility requires that semantically meaningful visual objects be used as 
basis for image data representation, explanation, manipulation, and exchange. With 
images being represented in compressed format, it is desirable to perform retrieval 
operations directly in the code space without requiring image reconstruction. In fact, 
any search algorithms that require image reconstruction will be infeasible from a 
:0 practical viewpoint because of the huge amount of images in most image databases. 
The previous JPEG and many existing coding techniques focus primarily on the issue 
of compression ratio, paying minimal attention to the need of content based image 
retrieval. 

RICS, by its name, has a fundamental consideration to various types of regions in 
15 images. The regions referred to as ROI (Region of Interest) are usually user-specified 
primitive geometric shapes, such as rectangles or circles. The regions that define the 
visual objects are usually of arbitrary shapes. Regions can also be generated as the 
result of certain mathematical operations or transform properties (e.g. significance of 
transformation coefficients) used to partition the image into disjoint regions of various 
20 (e.g. tiles, hierarchical blocks, or arbitrary shapes). 

In designing the RICS system, DAC considered carefully the distinction between the 
concept of an object and that of a region. RICS is open to very general definition of 
region. A region is a 2D spatial identity with pure syntactic contribution to a code 
stream. In contrast, an object may contain semantic information. Therefore, the 
25 region is perceived as a more elementary and reliable description than the object. 
RICS is region based, not object based. RICS supports a rich set of region types, 
from primitive geometrical shapes, tiles, hierarchical blocks, to most general arbitrary 
shapes. 

The region based coding strategy effectively supports the content-based accessibility 
30 of imagery data. Specifically, this ability enables random access to the code stream as 
well as providing a processing channel for user defined regions of interest or tile 
based techniques. Supporting MPEG-4 object based accessibility is one of the main 

3 
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Objectives of the RICS design. Furthermore, region-based coding provides a natural 
bridge for 'transcodability' with JPEG and JBIG. 

Scalability 

Many applications require image data that is available at different resolutions or 
5 qualities for decoding. For example, in a progressive transmission process, the bit 
rate control mechanism should allow the image data to be transmitted in certain 
priority order, and be able to tmncate the remaining data flexibly, either upon request 
from the receiving terminal or upon channel limitation. Scalable image coding 
involves generating a coded representation (code stream) in a manner that facilitates 
10 the reconstnjction of the image at multiple resolutions or quality levels by scalable 
coding. 

Ideally, the control of scalability should be centralized in a single module as the last 
stage on the encoder side right before the code stream is fed to the communication 
channel. Furthermore, it is desirable that the scalability control module can 
15 completely handle the required processing locally, without any further request 
propagating back to any previous stages in the encoder. In this way the scalability 
control module avoids the need for multi-pass computation. 

The RICS architecture is designed to support three types of scalability: scalability in 
terms of pixel precision, spatial resolution, and regions. 

20 Compactness 

There is no doubt that the new image coding standard should offer a higher 
compression performance than the fonmer JPEG, especially at the low bit rate end. 
Integrating the compactness, scalability, and accessibility into a general purpose, 
flexible, and open architecture represents a challenge for JPEG 2000. The RICS is 
25 designed to provide a solution. The basic idea of region based coding is as follows: 
• The input image data, after certain transformation, becomes a set of image 
primitives. This set of primitives can be wavelet transform coefficients. DCT 
coefficients, other transform coefficients, or even raw image data. 
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• The image primitives are grouped into regions. Region definition can come from 
user defined ROl, from other application modules, or from certain automatic 
segmentation algorithms running in the primitive space. 

• Each region contains one or more independent coding units (ICU). The primitives 
5 in an ICU are encoded and decoded independently, without reference to primitives 

of any other ICU. This procedure is called the intra-region coding. The outcome 
of an ICU operation is a code block. 

• A multiplexer (MUX) is employed to integrate the code blocks into the final code 
stream. 

10 A Brief Description of the Drawings 

Figure 1 is a simplified RICS block diagram. 
Figure 2 is a detailed RICS block diagram. 
Figure 3 is subband decomposition schemes. 
Figure 4 is a Wavelet Filter Library. 
15 Figure 5 is a performance comparison of YUV and standardized color systems. 
Figure 6 is geometric shapes supponed by RICS. 
Figure 7 is a hierarchical partitioning. 
Figure 8 is the relationship between regions and objects. 

Figure 9 is threshold masks obtained from Level 1 Lena Wavelet Decomposition. 
20 Figure 10 is a Level 1 common mask obtained by thresholding the combined data set. 

Figure 11 is individual region masks obtained by separating raw the common mask. 

Figure 12 is individual region masks obtained by separating raw common mask. 

Figure 13 is a spectral magnitude of zigzag spectrum as a function of position. 

Figure 14 is common mask spectrum filter sizes and captured coefficient numbers. 
25 Figure 15 is a coefficient banding concept used to quantize the common mask 

spectrum. 

Figure 16 is Quantization Band Sizes for Common Mask Spectrum of size 128 by 
128. 

Figure 17 is Spectral Quantization Band Sizes for Other Common Mask Sizes. 
30 Figure 18 is a Mask Overhead for Grayscale Images. 

Figure 19 is a low pass Common Mask Approximation for general coefficient 
classification. 

5 
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Figure 20 is a Mask Quantization following the Spectral Content for Bit Allocation. 
Figure 21 is Translated Common Masks for Remaining Resolution Levels. 
Figure 22 is the three types of ICU structures. 
Figure 23 is a Pavement of a region using type-1 ICUs. 
5 Figure 24 is a Embedded quad tree flowchart. 
Figure 25 is VLC Codes for EQW. 
Figure 26 is the transcode with JBIG. 

Figure 27 is the 'trans-out' and 'irans-in' modes for transcoding with JBIG. 

Figure 28 is Wavelet Mallot Decomposition for Lossy Color Transform Data. 
10 Figure 29 is Level Priority Processing Order for Each Channel (Lossy Case). 

Figure 30 is Level Priority Processing Order (Lossy). 

Figure 31 is Level Priority Color Interleave Processing Order (Lossy). 

Figure 32 is Level Priority Color Processing Order (Lossless). 

Figure 33 is Level Priority Color Interleave Processing Order (Lossless). 
15 Figure 34 is Typical MUX List Data Structure. 

Figure 35 is SNR Progressive Bit Budget Distribution Scheme. 

Figure 36 is MUX List Overhead for Y-Channel 8-Level Decomposition. 

Figure 37 is MUX List Overhead for Sqrjare Images of Various Sizes. 

Figure 38 is MUX List Overhead for Square Images of Various Sizes. 
20 Figure 39 is General Region Level Priority Color Processing Order (Lossy). 

Figure 40 is Region Level Priority Color Processing Order (Lossy). 

Figure 41 is Color Interleave Region Level Priority Processing Order (Lossy). 

Figure 42 is Region Level Priority Color Processing Order (Lossless). 
Figure 43 is Transparent Region Level Priority Color Processing Order (Lossy). 
25 Figure 44 is Transparent Region Level Priority Color Processing Order (Lossless). 
Figure 45 is a graph representation of the transparent region level in color 
mode for the 4 DCT region channels. 

Figure 46 is a graph representation of the region priority level in color 
mode for the 4 DCT channels. 
30 Figure 47 is a graph representation of the absolute region priority level 
color mode over the 4 DCT region channels. 

Figure 48 is a graph representation of scaled region priority level color 
mode of the 4 DCT channels. 

Figure 49 is a graph representation of the region percentage priority level 
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in color mode over the 4 DCT region channels. 

Figure 50 is a graph representaiion of the mixed processing multiplexer mode 
of operation over the 4 DCT channels 
Figure 51 is a diagram of the basic lossless header structure. 
5 Figure 52 is a diagram of the basic lossy header structure. 
Figure 53 is a diagram of the basic header tag structure. 
Figure 54 is a block diagram of the basic bit stream syntax for normal modes 
of operation in the lossy case. 

Figure 55 is a block diagram of basic bit stream syntax for region modes of 
10 operation. 

Figure 56 is a block diagram of the basic bit stream syntax for mixed modes 
of operation. 

Figure 57 is a block diagram of the JPEG transcode/decode method 
Figure 58 is a block diagram of the JBIG transcode/decode method 
15 Figure 59 is a block diagram of the modified post filtering procedure 
Figure 60 are three gradations of post processing in digital still images.. 
Figure 61 is a representation Edge areas where the de-ringing filtering is 
selectively applies. 

Figure 62 is a table of Test flmage: Camera (256x256 grayscale) Test Image: 
20 hk (256x256 color) 

Figure 63 is a table of Test Image hk.raw results 

Figure 64 is a table of Test Image camera.raw (grayscale) results 

Figure 65 is a table of Test Image hk.raw (color image)results 

Figure 66 is a table of Camera raw (part a) results 
25 Figure 67 is a table of Camera.raw (part a) results 

Figure 68 is a table of HK raw results 

Figure 69 is a table of: Quality versus iterations. 

Figure 70 is a table of: Quality versus iterations.As discussed. 
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A Detailed Description of the Drawings 
The System Architecture 

The architecture of the RICS system can be described as the scheme of dynamic 
5 streaming of code blocks. In this design, the image primitives of ICUs generate unit 
code blocks. A code block is a logically independent unit of coding and decoding 
which does not rely on information contained in any other blocks. Compression is 
achieved in each block coder, and the coding efficiency of the block coders 
determines the overall efficiency of a RICS system. The openness of the system is 

10 reflected in the different coding algorithms that can be used to produce different code 
blocks. The system uses a multiplex structure (a MUX) to assemble the code blocks 
into the code stream and realize the bit rate allocation. In short, the RICS architecture 
allows for flexibility in the areas of compactness and openness to the block coders 
and, at the same time, allowing the multiplexer to handle the various schemes of 

15 scalability and random accessibility to the code stream. 

It should be noted that the block used in RICS is a logical unit rather than a geometric 
concept. A code block may correspond to an 8x8 tile (in the case of JPEG mode 
coding), a pyramid data structure in zero tree like coding schemes, a rectangular area 
in block based coding schemes, or an arbitrary shaped area in the raw image buffer. 
20 The independence of encoding and decoding is the primary requirement of a code 
block. A code block is the outcome of an intra-region coding. 

Figure 1 illustrates the simplified RICS coding architecture. A more detailed diagram 
is shown in Figure 2. A detailed description of the function of each module is given 
25 in this document. Because the RICS is intended to be an open architecture, DAC does 
not limit each module to any specific algorithm. Instead, any individual algorithm 
can be placed into a module as long as it meets the functional requirement of that 
module. Supporting algorithms are included for certain modules, which reflect 
preferred operation of the overall system. 
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Types of Input Imape 

As shown in Figure 2, the RICS architecture supports the coding of three types of 
image data: grayscale, color, and bi-level images. 

Transformations 

5 In a typical multi-resolution coding scheme, an image is transformed via a multi- 
resolution decomposition process. In the proposed architecture, transforms such as 
KL, wavelet, wavelet package, lifting scheme, etc. can be placed in the transformation 
module. These transforms produce a set of decomposition coefficients {Cij} at 
different resolution levels and in different spatial orientations. 

10 The RICS architecture also supports DCT or windowed Fourier transform as a 
transformation technique. This is mainly for the transcodability with JPEG. It should 
be noted that the Fourier based transforms have been studied for more than a century; 
its theory is relatively complete and its mathematical and physical properties are well 
understood. Particularly, its translation, scaling, and rotation properties may be very 

15 useful for content based retrieval computations. On the other hand, wavelet 
transforms are relatively new, and many of its properties require further 
investigations. As a result, the support of DCT may have an impact beyond the sole 
backward compatibility consideration. 

RICS also allows the NULL transformation (that is, no transformation is applied at 
20 all). In this case, an identity transformation is applied to the raw image data as the 
transformation step. The NULL transformation is useful in several instances. For 
example, it is usually not beneficial to apply DCT or wavelet transforms to bi-level 
images (text) for compression purpose. Another example is the residual images in 
video coding. The information in a residual image is the difference between video 
25 blocks and has a high frequency spectral content already. It is highly questionable 
whether it is beneficial at all to apply another mathematical decomposition (DCT or 
wavelet) to this type of data for the purpose of compact coding. 

Region Definition 

TTie function of this module is to partition the coefficients produced by the 
30 transformation module into a number of spatial regions. The RICS supports three 
region schemes. 

9 
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1. Automatic partition based on the preordering of the coefficients. 

2. Partition based on user defined ROI or object related shapes. 

3. Partition based on tiling. 

The first partition scheme is also referred to as the region hierarchy formation process. 
Consider the transformation coefficients as a set. The region hierarchy formation 
process partitions the set into a number of hierarchically disjoint subsets according to 
certain definitions of significance. In providing a very general partitioning technique 
captured in a very general region based control architecture, RICS can perform highly 
flexible progressive transmission modes of operation that depend on data priorities set 
up for the code stream. 

Scheme 2 deals with user specified ROI, typically primitive geometrical shapes, such 
as rectangles or circles, as well as object related shapes. Object related shapes could 
come from a variety of sources, such as user input, motion analysis in video 
compression, etc. 

Scheme 3 is a simple partition and requires minimal for shape coding. This scheme is 
essential in the JPEG mode. Some wavelet based compression techniques utilize this 
scheme to explore coding efficiency. This scheme offers very little support for 
content based accessibility. 

Hierarchically disjoint regions can be used in combination with user defined ROI in 
still image processing and objects in video processing. However providing a fair user 
partition for detail information is difficult in still image compression. But automatic 
partitioning or preordering techniques can be performed to control user selections in 
both arbitrary and block based modes of operation. The research presented in this 
document introduces a new multi-level control architecture for advanced multi-level 
image processing. The ability to perform a host of partitioning and preordering 
techniques in both normal and region based modes of operation is encompassed in 
this architecture. Both arbitrary and automatic region formation schemes are handled 
in the same manner at a high level. 

Region Shape Coding 

Three types of region shape are supported in RICS: tiling, primitive geometric shapes, 
and arbitrary shapes. 

10 
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1. Tiling requires only a small set of parameters to describe the configuration 
especially in sequence based processing modes. However the sequences can be 
organized and presented in many ways in packing them into the final code stream. 

2. Primitive geometrical shapes can also be coded efficiently. For example, a 
5 rectangle can be defined by four integers (x mii^ mii) and (j:ma«,^maj) , or 

ixmin,ymi^ and (width, height) etc. 

3. Arbitrary shapes are more difficult and costly to encode. Partitions produced 
automatically by image analysis algorithms may contain many small regions. 
Specifically, the auto-region detection routines presented here produces 

10 hierarchically organized data partitions. This presents a highly challenging 
problem for shape coding. RICS provides two practical solutions to this problem. 
One is a quad tree based hierarchical region description. This mode of operation 
follows bit level ordering at different resolution levels (see embedded quad tree 
wavelet (EQW) in Ch. 4.) The other is a DCT coded multi-level region channel 

15 definition followed by preordering the partition (given the restraints of the code 
stream in terms of overhead). Though the coding efficiency of the second solution 
is currently slightly poorer than the first, it does contain a number of attractive 
features: 

• It provides a single representation for multi-level bitmaps. 

20 • From this single representation, region masks can be reproduced at arbitrary 
resolution levels, which is useful for subband coding. 

• It generates curved shapes, which potentially (e.g. in case of large number of 
complex curved regions) could outperform quad tree based region 
descriptions. 

25 • It has several useful properties, such as translation and rotation properties, 

which the quad tree based description does not offer. In fact, a quad tree 
representation will change dramatically if an object is slightly translated in the 
image, making this type of representation not suitable for content based 
retrieval applications. However, the main advantage of the quad tree 

30 representation is its simplicity and most noticeably in block based modes. 



11 
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The region shape code is included in ihe code stream when the region definition is 
determined at the encoder. In the case of decoder specified regions, the region shape 
coding has a whole new meaning. 

Inira-Repion Coding 

The function of the intra-region coding is to arrange Ihe transformation data in an 
arbitrary shaped region into a one dimensional code stream. Regardless of the region 
definition scheme or the region coding technique, this streaming process requires an 
intermediate state where a control architecture can be designed to tailor the region 
channels whether dealing with a bitmap mask, an auto-detection routine or any other 
of numerous classification techniques. At the decoding end. the inverse routine 
generates the same mask to unpack the values from the one dimensional code stream 
and place them back into the correct position in each region. 

Intra-region coding is completed in block coders. Different block coders can be used 
to produce the ID code stream. For JPEG mode, the zigzag scanning/quantization 
routine is called to pack an 8x8 DCT block. For wavelet based coding, both explicit 
quantization and implicit quantization schemes can be used. In particular, embedded 
zero tree and embedded quad tree can be used as implicit quantization schemes. 
Furthermore most implicit quantization schemes are implicitly decodable. In dealing 
with bi-Ievel images, a JBIG routine can be called a block coder. This effect can be 
staged by calling an implicit quantization scheme using one bit plane. Altematively. 
efficient JBIG routines can be called block coders in an embedded coding scheme at 
each of the multiple bit planes. 

The Multiplexer 

The function of the multiplexer is to assemble the code blocks derived from different 
subbands and different regions in proper orders into a single code stream. Due to the 
richness of region definition and block coding schemes, there are plenty of ways to 
pack the code blocks. Different ways to merging the data lead to different 
transmission priorities. For all transmission ordering modes, the final code stream can 
be transmitted progressively, and can be truncated at any desired place. 
The notion of using a multiplexer has been adopted in some standardized processes 
such as MPEG. In contrast, the use of a multiplexer in the design of a still image 
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coding system is rare. The region based coding strategy opens the opportunity to 
systematically explore the syntactic and semantic richness in code stream ordering 
and transmission medium. With the image segmentation techniques improving in the 
future, region shapes will become more accurate in tracking objects in the image. In 
5 that case, the multiplexer will not only work as a syntactic composer, but also impose 
semantic meanings to the code stream. 

Hi^hlipht of Features 

The RICS is a true open architecture. It supports not only the algorithms DAC has 
developed, but also can accommodate most existing well known compression 
10 algorithms. Users may include their own functions that are appropriate to their 
application in a number of modules such as transformation, region definition, region 
shape coding, and intra-region coding all under the MUX control architecture. It is 
also implicitly flexible to new technological advances. 

• DAC's current implementation offers superior low bit rate performance that is 
15 competitive with best existing compression techniques. 

• The richness in region definition in the RICS allows great accessibility, thus 
providing a solid foundation for content based image applications. 

• It covers compression for both continuous tone and bi-level images in a single 
unified architecture. 

20 • It provides lossy and lossless compression in a single, natural, code stream in the 
course of progressive decoding. 

• It provides a variety of progressive transmission modes that allow images to be 
reconstructed with increasing pixel accuracy or spatial resolution using region 
priority modes for user specified ROI or system defined significant areas. 

25 • It supports both fixed rate and fixed size modes. 

• It supports very flexible random access to and processing for regions with 
arbitrary shape. 

• It provides a graceful backward compatibility (or transcodability) with JPEG. 

• The generic region definition in RICS is a very suitable interface for object-based 
30 coding schemes currently under development by MPEG-4. (In fact, as the 

13 
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arbitrary region shape begins to fit dynamic objects witii more accuracy, motion 
estimation is more definite, and consequently allows for more efficient error 
compensation using region-based coding.) 

• Our general region shape definition provides a solid basis for object based 
5 composition and object based information embedding. Multiple objects with 

arbitrary shape are accepted. 

• Our bit stream is robust to bit errors. Each unit structure used by the MUX is 
independently decodable. 

• Incorporating standardized encryption techniques with the RICS system is 
10 straightforward. 

Transforms 

In order to support multiple application needs, provide transcodability, and 
accommodate future growth, the RICS architecture is designed to support three 
15 categories of transformation: the DCT, wavelet transforms, and a special NULL 
transform. 

DCT 

The RICS architecture supports the DCT transforms as defined in the baseline JPEG. 
Wavelet Transform 

20 The RICS architecture supports various kinds of subband decomposition schemes, 
including the three schemes in Figure 3. 

The Null Transform 

In dealing with the compound documents with mixed contents, it is sometimes 
required to encode some areas of the image without any transform. In particular, 
25 regions with bi-level pixels usually do not need to be transformed for coding purpose. 
In coding these regions, the transform stage is bypassed; the effect is simulated by a 
NULL transform stage in the process. 

The present version of RICS uses totally 34 types of wavelet kernels as given. in the 
30 table in Figure 4. In our experiments, the biorthogonal filters (Bior9-15 of the table) 

14 
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seemed to give the best compression. Both low pass filtering and high pass filtering 
are done using convolution. After filtering, the wavelet coefficients are down 
sampled by 2. This process is repeated using the low pass part until the desired 
decomposition level is reached. In inverse wavelet transform, quadtrature mirror 
5 filters (QMF) for low and high pass are used. Then the coefficients are up-sampled 
by 2. 

Color Transform 

In addition to YUV format, the RICS system uses the Karhunen-Loeve Transform for 
color transformation. Following the notation of statistics we term this process as the 
10 color standardization. The KL color transfonn requires more computation than YUV 
transform. Figure 5 shows the PSNR comparison of the two color transforms. More 
test results are also available from DAC. 

Region Definition and Shape Coding 

15 The notion of region processing plays a fundamental role in the operation of RICS. 
The choice of region based coding is motivated by many application needs such as 
content based retrieval, interactive multimedia application, graphics object and image 
composition, and coding of dynamic object in video compression. 
There are two fundamental issues in a region based coding strategy: defining regions 

20 and coding the regions. In the RICS system, region coding is divided into the shape 
coding and the intra-region coding (the content coding). Section 3, describes the 
various schemes for region definition and shape coding. The intra-region coding is 
discussed below. 

From the viewpoint of JPEG 2000, the process of region definition does not have to 
25 be standardized. However, schemes for forming regions that the JPEG 2000 will need 
to support should be defined clearly. The RICS supports three types of region 
definition. 

• Tile based organization 

• primitive geometric representations 
30 • arbitrary shapes 
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Region definition is an optional step: an image can be coded without specifying any 
region (the non-region mode). In this case, the entire image area is considered as a 
single region. 

Tiling: Definitions and Shape Coding 
5 Tiling is perhaps the simplest region definition scheme. In this scheme, the entire 
image area is divided into a number of rectangular blocks. In particular, 8 by 8 tiles 
are used in the JPEG coding mode. Most tiling approaches cannot cover the natural 
shapes, region trends or partial objects in a picture. 

For coding a tiling scheme, a small set of global parameters will be sufficient. The 
10 primary parameter is the size of tiling block and the technology is developed around 
this theme. In case of subband decomposition, tiling block size may vary from one 
subband to another. However, the RICS architecture has a sound operational 
foundation. It can be operated in both tile based or arbitrary processing modes of 
optional arbitrary region formation schemes. 

15 Primitive Geometrical Shapes 

Primitive geometrical shapes are ideal for supporting user interaction; i.e. user 
specified ROI, and for processing compound documents where the text areas can be 
well covered with one or more rectangular boxes. Geometric shapes currently 
supported by RICS are listed in the Table in Figure 6. 

20 

Arbitrary Regions and a Generic Binary Partition Hierarchy 

Let Co be the set of image primitives. Given an ordering relationship generic 
binaries partition of into C.^and C„ is defined by the following conditions. 

(1) C,oUC„ =Co. 

(2) C,onC„ =<D,and (Eq.III.l) 

(3) for any o € C,„and A e C,, , a)^b holds. 



Recursively, each of the new generated subsets can be further partitioned into two 
30 subsets, resulting in a hierarchical structure to represent the panition (see Figure 7) 
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A useful ordering relation popular in the compression community is "more significant 
than": a yb whereby a is more significant than b. The definition of significance can 
be very general. One definition may be used to create an ordering, or several 
5 definitions are combined together to generate an ordering. This can be illustrated by 
the following example shown in Figure 8. Suppose the region masks shown in the 
leftmost image are generated by considering the magnitude of wavelet transform 
coefficients as the measure of significance. The region mask in the middle image is 
an object shape generated by another definition of importance that, for example, can 
10 be the intensive dynamic region in a video frame. The rightmost image shows the 
intersection of the two region masks, which is a new region partitioned according to a 
multiple definition of significance. 

As an example, if the significance is measured by magnitude of wavelet transform 
15 coefficients, then the ordering relation is simply the numeric relation. To create a 
binary partition, a set of thresholds is sufficient. The first threshold, say , separate 
the initial set into two subsets, C,o and C„. Subsequent threshold values can be 
specified to recursively partition the new subsets. 

The above partition scheme produces accurately a number of disjoint subsets which 
20 satisfies the strict ordering relationship of >■ . However, it usually produces many 
small, scattered regions. Consequently, the representation of region shapes for these 
subsets can be potentially expensive. From the representation point of view, fewer 
numbers of smoothly shaped regions are desirable. In the following, a less strict 
partition scheme is defined, which takes the above partition scheme as a special case. 

25 Instead of requiring the condition a >- b to be held for all aeC^g, beC„, it is 
replaced by a looser condition 

-tZ? 6 C, , A V a , for any a € Co . (Eq.III.2) 

Intuitively speaking, this definition defines a partition in which no element in C, 
proceeds any element of C^. In other words, some elements belonging to C, (by 
Eq.III.l) may now be placed in C^, but, no element belonging to C, (by Eq.III.l) 
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can be placed in C, (Eq.III.2). In doing so. the subset will be enlarged, and the 
region shape of C„ will be smoother. At the same time, the enlargement of is 
controlled to such an extent that the loss of classification accuracy is limited to a 
reasonably low level. 

5 Another scheme, which is symmetric to the Eq.III.2, can be defined similarly for the 
following condition. 

10 e Q , fl ^ Z> , for any i) € C, . (Eq.III.3) 

In this case, we ensure that no element in Cq is proceeded by any element in C, . 
It should be pointed out that there is a tradeoff between the enlargement of and the 
space requirement for the region shape representation. It is possible (and indeed quite 
15 often) that there are some isolated scattered points that belong to Co . But, in order to 
exclude those isolated points from C,, we have to sacrifice many elements from C,. 
This will result in a very large Co , a degraded classification. An approximation to the 
scheme of Eq.III.2 is introduced to rectify this problem. A small number of elements 
in Co are allowed to be "misclassified" into C,. Those misciassified elements can be 
20 separated later in the intra-region coding stage that generally follows natural 
processing orders available for further classification of the partitioned elements. 
When such region hierarchy fomiation schemes are applied to image primitive planes, 
they serve as preordering processes that set up a partial ordering among the 
primitives. Then, the problem of coding this ordering becomes a geometry problem 
25 for coding the shapes of the partitioned regions. Transformations such as wavelet are 
known to use a joint spatial/frequency domain methodology. The region hierarchy 
formation scheme mentioned above has an intuitive geometric explanation for this 
class. Transformations such as DCT and Fourier transform are purely frequency 
domain methods. It is interesting to note that when the above partitioning scheme is 
30 applied to a 2D DCT transformation plane, the generated regions are roughly 
circularly shaped bands. The popular zigzag sequence used in JPEG turns out to be a 
simple and reasonable approximation of these bands. Therefore, the region based 
coding strategy accommodates naturally both spatial/frequency domains and pure 
frequency domain transfonmations. 
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Automatic Region Formation in the Wavelet Transform Domain 



DAC has developed an automatic region formation scheme that is used to categorize 
wavelet coefficients in the transform domain according to magnitude. The procedure 
will be outlined in the following sections. Also included are many of the design 
5 details used to develop the first version of the general region classification scheme. 
The advantages and disadvantages of the technique as well as performance and 
overhead issues will be discussed. Finally conclusions will be drawn to summarize 
the results. 

Detecting the Regions 

10 The wavelet transform is used to decompose the original image information into a 
multi-resolution hierarchy of data that is more suitable for compression. The result of 
completing one pass of the 2D wavelet transform is one low pass part (LL|) and three 
high pass parts (HLi, LHi and HH,.) The transformation procedure is repeated using 
the low pass LL part as the starting point for each succeeding pass. The LL part is an 

15 approximation of the original image at a lower resolution. At each resolution level 
the HL parts carry the vertical, the LH parts carry the horizontal and the HH parts 
carry the diagonal edge information. The level of detail information contained in any 
particular orientation is specific to the local response generated in the choice of 
wavelet kernel. The original image is transformed into a unique band pass response 
20 hierarchy of detail information. In any compression process, it is the highest energy 
coefficients that are considered first for making a contribution to the reconstructed 
image. It is this transformation properly that is used to partition the coefficients into 
regions of interest corresponding to high energy locations in the original image space. 

The magnitude of the coefficients in a particular subband is related to the response 
25 that the original image information has to the transform kernel. Sharp edges in the 
original image (or a lower resolution LL part) are indicated by an increased kernel 
response in terms of coefficient magnitude, in the general vicinity of a specified area 
of concern. Other image changes that are not as abrupt will translate into a gradual 
response that is less significant. In those areas where little or no change occurs, the 
30 wavelet transform will indicate little or no response. 

Level 1 Kernel Response Threshold Images 
19 
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DAC's automatic region formation scheme makes use of coefficient magnitude 
information to categorize the data. The starting point for the procedure is at the 
highest resolution level of the wavelet transform hierarchy looking at the three detail 
data sets (HLi, LHj, HHi.) A 2 bit representation of each detail orientation for a 256 
by 256 Lena image is given in Figure 9. These images are obtained by using 
threshold values to categorize the coefficients by magnitude. The wavelet kernel used 
in this case is the standard 9-7 filter set implemented in a lifting scheme. 

A decaying histogram procedure is used to determine the threshold values in each 
case. In this particular analysis the coefficients are partitioned into regions according 
to decreasing levels of magnitude; 10% of coefficients are partitioned into region 1. 
15% into region 2, 25% into region 3 and 50% into region 4. Note that the largest 
coefficients appear darker in the images while the smallest appear lighter. The 
threshold values are calculated automatically for each orientation. The values used to 
generate the partition in each case are given in Eq.III.4 (assume Cj is the magnitude of 
the coefficient under consideration.) 



HL,: 0 s Ci < 2, 2 < Ci < 4, 4 < Ci < 8 and Q ^ 8 

LH,: 0 < C, < 3. 3 < C,- < 5, 5 < Q < 14 and Q > 14 (Eq.III.4) 

HH,: 0 < Ci < 2, 2 < Ci < 3, 3 < Q < 5 and Ci > 5 



The set of threshold images appearing in Figure 9 is motivational in determining a 
region formation scheme based on this approach. It is apparent that wavelet 
coefficients can be classified in this manner to form a mask. However, the large 
overhead required to code each individual mask must be solved in order to make this 
type of classification scheme viable for implementation. 

Formation of the Raw Common Mask 

The first step taken in the development of a technique designed to reduce the mask 
overhead is to consider a common mask approach that can be used to capture the most 
important coefficients between all three orientations at the lowest resolution level. 
Given that the data range is different in each orientation, the data is normalized to the 
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largest range in an absolute value sense. This step is taken in order to compensate for 
the range differences that exist between the data sets, and to ensure that corresponding 
pixels from each orientation can make an equal contribution in the formation of the 
common mask. This step is illustrated in Eq.IlI.S. 

HL", = NewRange(HL,) = ScaleRange(MaxRange(LH,, HL, , HH,)) 

LH', = NewRangetLH,) = ScaleRange(MaxRange(LHi, HL, . HH,)) (Eq.111.5) 

HH", = NewRangetHH,) = ScaleRange(MaxRange(LH,, HL, . HH,)) 



Once the data ranges for the smallest two sets have been scaled to that of the largest, 
the common mask values can be determined. The next step is to generate a new 
normalized data set by taking the maximum value between scaled orientations at each 
15 pixel location. Let C (i) be the new coefficient obtained in this procedure and let 
C HLi . C LHi and C hhi be the individual scaled values from each orientation. The 
new data set is formed using the step illustrated in Eq.III.6. 

C (i) = MAX(ABS(C HLiO)). ABS(C HLi(i)). ABS(C hliO))) (Eq.III.6) 



20 The final step in forming the common mask is to threshold the new coefficients based 
on the histogram decay and the percentage of the total data amount to encapsulate in 
each region category. The same values mentioned earlier are used in this case (10%, 
15%. 25% and 50%.) The raw common mask categorization image appears in Figure 
10. The threshold values used to generate the mask are given in Eq.III.7. 

H',: 0<Ci<5,5<Ci< 11. 11 <Ci<26andCi>26 (Eq.III.7) 



30 Note the level of detail contained in this image. High frequency edge information 
representing the most important coefficients appears in the darkest areas. Low 
frequency changes representing areas of the image that are relatively uniform appear 
in the lightest areas. The combined kernel response from each orientation follows the 
high energy areas that exist in the original image. 
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Generally, the raw size overhead of the common mask obtained in the previous 
section is still too large to be included in the final bit stream for most image 
processing needs (i.e. an overhead of 2 Bpp for this particular 4 region mask). This 
5 section outlines the next step in the automatic region formation scheme. The DCT 
transform is used to reduce the overhead size to an acceptable amount. 
The DCT is not applied directly to the common mask data in its mult-valued form. 
Rather it starts by considering each region contained in the common mask as an 
binary data set such that the sum of the individual data sets is equal to the multi- 
10 valued mask in functional form. The DCT transform applies the transform of the sum 
of functions in the same manner as taking the transform of each individual function 
and summing the results. Furthenmore there are unique spectra in each case that 
should be considered separately for a fully integrated analysis. The individual region 
masks obtained by decomposing the raw common mask for the Lena image appear in 
15 Figure 11. Notice that although there are only three images, the fourth region channel 
(the background) is implicit. 

The 2D fast DCT algorithm is used to transform each binary data set into the 
frequency domain for spectral analysis. The spectra for the upper three regions 
appear in Figure 12. Please note that a log scaling of the coefficients is used in the 
20 images to exaggerate the appearance for display purposes. Notice how the low energy 
content is dispersed through the spectrum in each case. Even at such low energy 
levels, the energy is localized in some areas. 

It is apparent from the spectral results of Figure 12 that the magnitude of the DCT 
coefficients decreases rapidly in moving down each spectrum. The interaction of the 

25 content retained for each partition is used to guide the classification procedure. The 
interrelationships that exist in each binary spectrum must be investigated further. 
A plot of the overall summed spectra in short integer format appears in Figure 13. 
The conventional zigzag ordering is used to re-group the coefficients before 
quantization. Note that only the first half of the data in reflected in the plot. The 

30 outer 8K are omitted since they do not change much from the DC level. 



22 



Copied from 10340491 on 04/01/2005 

BNSDOCID: <W0 00<9571AE I > 



wo 00/49571 PCT/CAOO/00134 
The first consideration observed in the plot is that the coefficient magnitude 
deteriorates at an alarming rate over a relatively small number of coefficients in the 
first portion of the spectrum. The second observation is that most of the spectral 
energy is concentrated in the first 12%- of the data. Together these observations 
5 suggest that implementing a low pass spectral filtering stage, followed by careful 
quantization will reduce the amount of data that must be retained to capture the region 
formation concepts at the lowest level ( in this common mask partitioning scheme). 
Note the apparently large dc offset. Part of this level is introduced by scaling the 
original wavelet data sets to a unified range. The rest of the offset is due to the DCT 
10 transform. 

Low Pass Filtering the Common Mask Spectrum 

The total mask spectrum must be filtered carefully to reduce the excess content that 
has little contribution to the region trends that exist in the retained coefficients. The 
main concern is to determine how to quantify the filter size such that a dynamic 
15 implementation may be obtained that will work accommodate any arbitrary image 
sizes. 

The first step taken towards finding a dynamic solution to this problem comes from 
experimentation. In looking at the result of using zigzag ordering, coefficient 
magnitude and significance deteriorate at some exponential rate. Experimental 
20 evidence suggests that as the size of the original image increases, the number of 
coefficients that must be retained to construct a similar quality mask increase on a 
logarithmic scale. This is an important observation since doubling the original image 
dimensions translates into a log 2 increase in the number of retained spectral 
coefficients. 

25 Assume that p is a general log base to begin our discussion. If the filter size for an N 
by N common mask spectrum has been determined experimentally as <|»n, then what is 
the size of the filter for a 2N by 2N mask spectrum? If the unknown filter size as ^n, 
then the following simple logarithmic relationship can be used to compute the new 
filter size. 

128 

-TTT- =L08(>{it.N/4.:N) 
35 (Eq.III.8) 

lj>2N = •t'N ■ P 
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A number of experiments were conducted to determine a reasonable filter size that 
5 can be used as a base for filtering common mask spectra. It has been determined that 
for masks of size 128 by 128, a filter size of 40 to 45 diagonal rows yields reasonable 
results. The table in Figure 14 give the filter sizes for raw mask sizes using 1.375 for 
p. The logarithmic base P can be used for final calibration in the current design 
implementation. 

10 Notice the filter size specified for a common mask spectrum of 32 by 32 is given as 
32 diagonal rows. This value is selected to calibrate the filter sizes for the larger data 
sets. The results indicate the filter size required for a first cut of the common mask 
coefficients does not have to be exact since the distribution of the spectral content 
may change slightly from one image to the next. However, it must be large enough to 

15 capture the majority of the most important spectral content for each region partition. 

Spectral Quantization Techniques for Filtered Coefficients 

DAC has developed two techniques for quantizing the filtered mask spectrum. The 
first approach is basically a modified uniform quantization procedure that is used to 
20 verify the concepts. It proved useful in confirming the algorithms and developing the 
concepts, but it lacks robustness and flexibility. The second technique is a general 
approach that is much more robust and elegant. It is dynamic since it tracks the 
magnitude of the spectral content in determining how to quantize the coefficients. 

A First Quantization Approach 

25 The first approach combines the filter dimension determination into the quantization 
procedure. A base mask size of 128 by 128 is use as the starting point in this 
approach. Assuming that a logarithmic relationship exists for coefficient importance 
along the diagonal order, a quantization technique can be developed to exploit the 
relationship. 
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The main premise is tliat there are bands of spectral coefficients that can be used to 
guide the quantization procedure. The importance of the coefficients in each band 
decreases along the spectrum. This concept is illustrated in Figure 15. 

The number of diagonal zigzag rows to include in each band is a function of the filter 
5 size and the original common mask size as well as the spectral content. As a first 
approximation, a table look up technique is used to determine the number of rows to 
use in each band. The base band sizes determined experimentally to give reasonable 
partition and reconstruction results for a 128 by 128 mask are given in the table in 
Figure 16. 

10 

Notice that the number of diagonal rows doubles in each step. The number of bits 
selected to quantize each band decreases from one band to the next in the following 
fashion. 

• Band 1:16 bits 
15 • Band 2: 10 bits 

• Band 3: 8 bits 

• Band 4: 8 bits 

A series of band sizes were found experimentally of mask dimensions by powers of 2 
increments. The results are in the table of Figure 17. These values are used to 
20 partition the spectral content for quantization. This information together with the 
number of bits used in each band and the filter size measurement determines the 
common mask overhead of an image. 

The mask overhead based for this quantization scheme is tabulated in the table in 
Figure 18 for the corresponding image sizes. There is an additional overhead header 
25 included in this result for mask reconstruction on the decoder side. Note that the 
mask overhead is the same size for both grayscale and color images since it is 
generated from the intensity information. However, the percentage overhead in the 
color case is reduced by a factor of 3 since it can be distributed over 3 channels. 

After applying the inverse DCT on the quantized coefficients a low pass 
30 approximation of the original common mask is obtained. It is this result that is used 
to guide the region channels in classifying and processing the wavelet coefficients. 

25 
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The result indicated in Figure 19 is obtained by applying this technique to the first 
wavelet transformed level for a 256 by 256 Lena image; From Table 17, the overhead 
for this mask is 1021 bytes. 

Notice the high energy changes have been eliminated by low pass filtering (especially 
5 the changes that occur over small portions of the image). However, the mask does 
succeed in capturing most of the essential parts of the original image. This mask is 
common to all orientation at this level and used to partition the coefficients into 
region categories. 

A Second Quantization Approach 

10 There are some apparent problems with the underlying procedures used for 
quantization in the first approach. The first of these is that the coefficient band 
classification is fairly rigid and not that flexible. The filtering and quantization stages 
are tied together and fully independent. Affecting one parameter has consequences to 
the other. It is true that a logarithmic relationship exists in the spectrum that can be 

15 exploited to determine a reasonable filter size, but the quantization step must be 
tailored from the filtering stage. 

DAC is currently implementing another quantization scheme for mask filtering that is 
much more dynamic and flexible. Quantization is no longer tied to the filtering stage 
and the coefficient banding technique of the previous approach. The current 

20 technique follows the actual content of the specuvm in forming quantization 
categories. In this way, coefficients are grouped into categories according to the 
energy gradient of the spectral distribution. The concept is illustrated in Figure 20. 
The rate of spectral decent together with position information within the spectrum can 
be used to quantize the number of bits used to categorize the data. This process can 

25 be controlled precisely. The procedure can be calibrated as required by studying the 
effects of changing image dimension in terms of an optimal quantization change. 



The previous quantization technique can be used to estimate the mask overhead for 
the new scheme proposed here. Consider the same base conditions presented for the 
banding approach applied to the 256 by 256 Lena image. The spectral filter size in 

26 



Copied from 10340491 on 04/01/2005 

NSDOCID: <WO 0&19571A2.I_> 



wo 00/49571 PCT/CAOO/00134 
this case is 44 diagonal rows. If 4 quantization intervals are used at 16. 12, 8 and 4 
bits each the mask overhead can be estimated. Let Ni be the number of coefficients in 
each interval, bi the number of bits used to quantize this interval, and Om be the total 
mask overhead. A rough estimate of the new packed size is calculated in Eq.III.9. 

3 

= 2 Ni.bi 

' = " (Eq.IIl.9) 
= 6 X 16 + 39 X 12 + 186 X 8 + 804 X 4 

= 5268 bits (- 659 bytes) 
15 

Even though the actual overhead has not yet been determined, from the insight gained 
in using the first quantization approach a saving of 30% is not unreasonable. This 
promises to be a substantial saving in the mask overhead. There are ways to utilize 
20 saved overhead. The first is to reduce the amount of mask information in the code 
stream. This approach will retain the same mask quality leaving extra space for 
refinement information in the code stream. The second way to utilize the saved 
overhead is to retain a better quality mask in keeping the overhead size constant. 

The Common Mask and Multi-resolution Classification 
25 A translation technique is required to apply the region mask at other resolution levels. 
The simplest approach is to look at the significance of the mask coefficients in a small 
block to determine an appropriate value for the next resolution level. This process is 
repeated until a mask for each resolution is obtained. The most logical approach is to 
take the largest coefficient through to the next level. 

Mvfuiju = Max(Mv^2i^ , Mv^jzi^,^ , Mv^jj^jj^, , Mv„j2i^,_2j+,) (Eq-ID.IO) 

35 

Some heuristic approaches have been investigated for resolution translation. One 
approach is to make decisions based on the sum of the individual mask coefficients. 
The upper level masks are affected by either enlarging or shrinking the encompassed 
40 coverage regions. However, there is currently no conclusive evidence confirming the 
validity of this approach and more experimentation is required. 
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It is worth mentioning at this point, that DAC has begun investigating the wavelet 
kernel ROI masl< formation technique of verification model (VM) 3.0 as a method of 
resolution translation for other transformation levels. Preliminary investigations 
indicate this technique introduces a gradient bias for regions of higher classification 
5 level in the translation. A combined hybrid approach may have some advantages. 
The rest of the common mask hierarchy obtained by using the simple 4 to 1 resolution 
translation techniques and observing the most important region as key is given in 
Figure 21. The complete set of region masks is used to classify regions of importance 
in the wavelet transform domain according to coefficient magnitude. 

10 Misclassification and the Common Region Mask 

It should be noted that while it is true that misclassification occurs in this technique, 
there is a tradeoff between code stream overhead and region channel accuracy. The 
original raw common mask size is 32 kBits for a 4 region mask. The basic uniform 
quantization scheme compresses this amount to 8 kBits. The modified scheme 
15 suggests that mask compression of 5 or 6 to 1 can be achieved with a tradeoff between 
mask accuracy and code stream overhead. 

The low pass filtering and spectral quantization steps can be tailored to deliver similar 
quality masks for most images in a dynamic implementation. In addition the partition 
can be studied in regards to how the original grouping is affected by the filtering step. 
20 For the Lena image example the original partition is 10/15/25/50%. After filtering the 
coverage ratios used for bit budget calculation and the code stream distribution (see 
multiplexer in Ch.V.) changed to approximately 9/31/33/27%. There may be some 
advantages in adjusting the original panition such that the data coverage ratios are 
controllable in the final stages. 

25 The question of miscoverage must be addressed in all lossy region formation 
techniques. It is the subsequent modules that are affected by the misclassified 
information. Thus the interaction between the DCT region formation module and the 
subsequent sorting and packing modules must be understood in any attempt to obtain 
the best region channel results for this classification scheme. 

30 
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Subsequent stages in the compression algorithm are already tailored to process 
information in a descending level of importance. Most wavelet implicit quantization 
procedures process the coefficients in bit plane order of importance. What this means 
5 in terms of any region formation scheme is that the information encompassed in each 
region is further partitioned by bit level processing. In the common mask approach, 
further processing causes the most important misclassified coefficients to filter to the 
top at a faster rate than misclassified coefficients of lower importance. Effectively bit 
level processing can be considered as a second partitioning stage that deals with 
10 misclassification. 

There is a tradeoff between the accuracy of the original region partition and the 
additional bit level processing overhead placed on subsequent sorting stages. It is the 
final processing stages that implicitly accomplish the final coefficient partitioning. 
DAC is currently using two bit plane partitioning core technologies for both grayscale 
15 and color images, and binary partitioning for bi-level images is now under 
implementation. 

Bit Plane Sorting for Common Mask Miscoverage 

The concept of bit plane sorting is a well know technique used to organize priority in 
a distributed list of data. ID sorting a list of information generates a map specific to 

20 that list in an optimal fashion. The largest values are mapped first as the list 
transposition routine progresses. The order of decent is the standard order used to 
classify information if considered from the normal processing standpoint where no 
regions are involved at all. As a result, most of the misclassified information is 
considered for further classification and are thus biased to be included first in the final 

25 code stream order. The only difference is that now there is a region description 
overhead as well as the final mapping overhead of sorting the region of interest 
partitions in a hierarchical fashion. 

Currently DAC has not documented the exact overhead introduced by region 
processing in comparison to normal processing modes that do not use regions. But 
30 the final results obtained by using the common mask approach are encouraging to say 
the least. At this level the underlying routines do not care where the information to be 
sorted originated. The ordering technique will run to completion whether region 
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channels are involved or not. The original misclassified coefficients generated in the 
first cut will filter through to the final ordering set by the sorting stage. There is an 
important tradeoff here between the soning cost and the region overhead that must be 
considered in the final design implementation. 

5 Bit Plane Ordering for Common Mask Miscoverage 

DAC has developed a quad tree ordering technique (see EQW Ch.4) to track 
coefficient importance in both normal and region based processing modes of 
operation. The normal quad tree approach is to track the leading one position through 
each block under consideration in a recursive fashion. A map is generated in the 
10 process for each entry based on importance such that decoding follows the same 
order. At this basic level the quad tree order is an efficient method of processing the 
multi-resolution decomposition image information in an effective manner. 

The tracking or positioning overhead introduced by mapping the data is generally 
smaller than in the ID sorting case. Retaining vertical and horizontal positioning 

15 information is highly advantageous (in most image processing applications) since the 
underlying technology is generally consistent in this medium. The ID sorting 
approach sacrifices the natural vertical correlation that may exit (although it can be 
partially recovered by bit level entropy coding in the soning routine.) However each 
of the two soning techniques can be used from an operational standpoint to produce a 

20 similar net effect. One may produce slightly better results than the next, but the 
control architecture is consistent. Each final ordering routine can be tuned to meet 
specific processing needs. The ID sorting approach may be a useful candidate for bi- 
level image processing, and DAC is currently investigating this possibility. 

The quad tree approach has currently been modified to operate in arbitrary region 
25 processing modes. The ID sorting routine tracks the importance in a list of 
information specific to the block under consideration. That is the region partition is 
distributed in ID hierarchical lists of information. The original preordering of the 
coefficients is captured in the common mask procedure. Thus the decoder map 
already exists for the inverse ordering. 

30 Initially DAC's region processing quad tree begins by building an information map 
for each region channel to guide the classification. The core engine has been 
modified to drops blocks not contributing to the final ordering in each region 
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category. Currently the overhead comparisons for each sorting case have not been 
documented. The results will be available from DAC once the region based quad tree 
design implementation is finalized (currently under implementation at DAC.) 

In each case, the misclassified coefficients obtained from the original preorder mask 
5 partition are considered again for ordering in the sorting stages that follow. 
Effectively, the most imponant misclassified coefficients still have the opportunity to 
make important contributions to the final code stream as well as improving the 
reconstruction quality measurements in the decoded image. The difference between 
this approach and most efficient classification approaches is that the original 

10 processing order used of the efficient scheme is somewhat distributed through each 
region channels. However the misclassification in each region channel is forgiven to 
a cenain extend the subsequent bit level sorting module. One of the prices that must 
be paid for any automatic region classification scheme is that the optimal processing 
orders that exist in standard approaches must be considered in a different context. 

15 The original ordering concepts used for standard procedures must be extended to 
include other ordering techniques specifically tailored for operation in arbitrary region 
processing. 

General Region Processing Concerns 

There is one important hurdle that must be overcome in the design of any arbitrary 
20 region processing technology. There are times when rigid tile or block control is the 
optimal way to process image information. However, it is difficult to achieve 
arbitrary accuracy in a strict processing order. DAC's ID and 2D sorting techniques 
address this issue in design implementation. In our opinion arbitrary region 
processing channels can be defined in a still image environment. Careful study and 
25 analysis suggests that this is a viable approach not only for the DCT common mask 
approach but also for other arbitrary region formations techniques that have not yet 
been developed. Furthermore region classification may have interesting properties 
that can transferred to video processing environment where the region of interest 
concept and the object motion vectors can be considered together in a unified 
30 approach to region of interest processing. 

DAC's region technology in its original form was designed to include region 
processing capabilities from the start. All internal modules can be operated in both 
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normal and region based modes of operation. The MUX architecture (introduced in 
Chapter V) is designed to meet operational requirements in both nomial and region 
processing modes. It is possible to extend the notion of block based image processing 
to a level where a block of information is interpreted in an arbitrary way that depends 
5 only on the underlying core processing technology. In this way, blocks can be 
considered as individual processing units each fully capable of making a contribution 
to the final code stream. Furthermore control architecture can be developed to fully 
exploit the region processing capability. All region formation schemes, whether 
arbitrarily defined user regions of interest, automatic region formation schemes or any 
10 other specially designed region classification technique, can be processed in a 
common framework. The same overall processing architecture can be tailored to 
operate in an optimal fashion for both normal and region processing modes of 
operation. 

Additional Processing Concerns 

15 One of the important concerns in the JPEG-2000 community regarding the 
development of a new compression standard is that it must be able to support ROI 
processing. There are a number of techniques that have been tested with some 
success. Some of these techniques are listed below. 

• Sequence based mode. When this mode of operation is used, each set of data 
encapsulated by a specific region is treated as a separate sequence. The sequences 
are coded independently. 

• Scaling based mode. In this mode of operation, the magnitude of each ROI 
coefficient is increased using a predetenmined shift factor. If multiple ROIs are 
used with increasing levels of imponance, multiple shift factors distinguish each 
ROI. The effect of this scaling technique is to force the ROI coefficients to be 
encoded first. 

• ROI from the middle. In this mode of operation, a concern of coding an ROI in 
the middle of the bit stream is addressed. This is useful for progressive 
transmission of images where the user is given the opportunity to focus on a 
specific ROI before the complete set of image data arrives at the decoder. This is 
important for medical image processing. 
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• ROI without sending mask information. This mode of operation is a special case 
for scaling based mode. The coefficients for each region are scaled such thai the 
magnitude of each ROI exceeds that of the ROI of less importance. The decoder 
knows which ROI the coefficients belong to based on the magnitude. 

5 The ROI coding technique that DAC has developed generally falls into the sequence 
based category. But internal modules have been tailored to operate in both block 
based and arbitrary block sized processing units. DAC is currently studying the 
different mask and region formation techniques to determine the effect each ordering 
and partitioning approach has on the overall organizational structure required for 
10 region channel processing. The current implementation includes a scalable ROI mode 
of operation implemented without actually shifting the ROI coefficients encompassed 
by the classification. Data scaling is handled in the MUX control architecture. 
Internally, the scaled version is a special case that can be implemented in MUX 
control. 

15 Preliminary study has begun for the last of the four categories cited above. ROI from 
the middle implementations can be addressed by using resolution a progressive 
processing order. This type of ordering is supported in the MUX architecture. DAC 
is currently considering a design implementation for client side region of interest 
requests in a resolution progressive architecture for client-server region interaction. 

20 ROI without mask information makes sense in theory. However it places a huge tax 
on the subsequent sorting and classification stages. Generally speaking, the overhead 
cost introduced by additional sorting is comparable to the size of a bitmap mask that 
could have been used to classify the data initially. The shift introduced to each ROI 
must be processed in the sorting stage. Eventually sorting must traverse from top to 

25 bottom. Coefficient scaling under tight control has some definite benefits in many 
cases. However it would be highly inefficient in a general region processing 
architecture because of the high sorting costs as a result of traversing all region 
channels. In a region formation scheme of 4 or more regions, the total bit level 
difference between all ROI is too large to process entirely. On the other hand, if the 

30 resolution of the overall partition is reduced (less bit levels), less room is available to 
form an accurate region classification. There is some common ground that must be 
explored in order to determine the optimal processing mix. 
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Intra-Region Coding 

The segregation of image primitives into regions provides a basis for accessing the 
image contents. Compact representation of image primitives is achieved via inira- 
5 region coding. 

Independent coding units (ICUs) are the building blocks for intra-region coding. In 

RICS, intra-region coding with ICUs is designed with a set of objectives. 

(1) Full data independence. From the notion of ICU, the encoding and decoding of a 

region is done without referring to the data of any other regions. 
10 (2) Allowance for multiple coding schemes. Each ICU may be coded using different 

coding schemes. New coding schemes can be added to the system (modular 

openness.) 

(3) High coding efficiency. Depending on the characteristics of the primitives in a 
given ICU, one or more efficient coding schemes can be selected to produce a 

15 code block with high compactness. 

(4) JPEG and JBIG transcodability. It is preferred to have a single JPEG 2000 
coding platfonm that also accommodates JPEG and JBIG modes of coding. 

(5) Scalability of code stream. 

(6) Error resilience. The ICUs are natural blocks for bit error localization. 

20 (7) Parallelism. There is no data dependence between ICUs and so the encoding and 
decoding of all ICUs can be performed in parallel. 

Intra-region coding in RICS involves the following steps. 

(1) Choosing a coding scheme. 

(2) Determining ICU structure. 
25 (3) Coding. 

(4) Pre-packing. 

In the rest of this section the first three steps are described, 
Choosing a Coding Scheme 
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For coding an ICU, the current RICS design employs six categories of coding 
schemes. 

• Zigzag DCT Coefficients Packing (JPEG baseline scheme) 

• Embedded Quad tree (or similar schemes) 



There is no general restriction on what schemes have to be used for coding a given 
10 ICU. However, certain preferred embodiments may apply. For example, if DCT is 
used at the transform step, then the zigzag DCT coefficient packing might be 
preferred (but not restricted to). 

Choosing ICU Structures 

While being a logic concept, an ICU still has a geometric structure. The ICU 
15 structure is directly relevant to the coding scheme chosen. Currently, the RICS 
recognizes three types of ICU structures (Figure 22). 

Type 1 ICU Structure: 8x8 Blocks 

This type of ICU is designed exclusively for use in JPEG mode, although it is not 
necessary that DCT transform coefficients be coded this way. 
20 If there is no region defined in the image and the JPEG mode is chosen, the entire 
image is paved with type-1 ICUs. 

If a region R is to be coded using JPEG mode, the set of ICUs that covers R form a 
unique, minimal pavement for/?, in the following procedure. 

Siep 1. Let j'min be the first row from the top that has at least one pixel in R, y nu> the 
25 last row, jcmin the first column from the left, and Xtmx the last column. 

Step 2. Let (j: min, j min) and (j: m..,^ m.) be respectively the top-left and the 



bottom-right comers of the bounding rectangle of R. 

Step 3. Sianing from the top left comer, pave the bounding box completely with 
type-1 ICUs. 



5 



Embedded Zero tree (or similar schemes) 



ID Progressive Sorting Schemes 
Block Based Quantization 
JBIG Coding Algorithms 
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Step 4. Remove those ICUs that cover no pixels in R. 

After this procedure, the remaining ICUs form the minimal pavement of R (Figure 
22). 

Type 2 ICU Structure: Rectangles 

5 Type-2 ICUs are rectangles. Except for the JPEG and embedded zero tree schemes, 
any other inira-region coding scheme uses type-2 ICUs to pave and code a region. 
There is no general restriction on the dimensions of a type-2 ICU. It is usually 
defined by the selected coding scheme. For example, both embedded quad tree and 
explicit block based quantization can use arbitrary sized rectangles, with some 
10 preferred embodiments (such as the 64x64 blocks used in EBCOT of VM 3.0 (A)). It 
should be noted that once a region is paved by type-2 ICUs. different coding schemes 
can be used for each ICU. 

For subband decompositions, no type-2 ICU is allowed to cross subband borders. 
Essentially, this means that type-2 ICUs support only various types of intra-band 
15 coding methods. 

Type-3 ICU Structure: Pyramids 

The third type of ICU is the pyramid structure, which is used by various inter-band 
coding methods, such as the embedded zero tree. In the current RICS design, there is 
one limitation on the region definition procedure for type-3 ICUs: a region rnust be 
20 specified in the top-down manner, from lower resolution to higher resolution in the 
decomposition pyramid. That is, for a given region definition R in LL, every element 
in ^ defines a set of type-3 ICUs (i.e. in three spatial orientations, respectively). 

Let C"(/,/) be a wavelet transform coefficient at spatial location (i,j) in the LL 
subband, C; (m,n) be a wavelet coefficient at spatial location (m,/j) in the subband 
25 at level / (1=1.2....) in orientation r (r=l,2,3). Assuming that a 5-level wavelet - 
decomposition is used, the three type-3 ICUs defined by C"- (/./) can be represented 
in the following fashion. 

In orientation 3 : 

[ClU.j)] U {C\{2i,Zj).C\(2i + \.2j),C\{2i.2j + l),C\(2i + \,2j + \)) U ... 
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In orientation 2: 

{Cs'O'.y)) U {C',{2i,Zj),Cl{Zi + l,2j),Cl{2J.Zj + \),C'j2J + \,2J + i)] U ... 
In orientation 3: 

{ClU.J)} U {Cl(2i,2j),CU^ + l2J),Cl(2j,2j + \),Cli2i + l2J + D} U ... 

5 Notice that the primitive is not included in any the three ICUs. Instead, the 

LL subband is coded differently, using the lype-2 ICUs. 

Coding Type-l ICUs 

Type-1 ICUs are defined exclusively for coding 8x8 DCT coefficient blocks. Coding 
of type-1 ICUs follows the JPEG baseline algorithms. 

10 Coding Type-2 ICUs 

Once a region is decomposed into a set of non-overlapping type-2 ICUs. each ICU is 
coded independently, and the collection of codes for all ICUs is the coded 
representation of that region. In this subsection we describe three techniques that may 
be used for coding type-2 ICUs. 

15 Using Embedded Quad tree 

Embedded Quad tree Wavelet (EQW), Figure 24 is an efficient and fast method for 
type-2 ICU coding. This technique implements an embedded progressive sorting 
scheme in a quad tree-like structure. In contrast to inter-band coding methods, the 
EQW explores the intra-band self-similarity of the wavelet decomposition 

20 coefficients. The EQW-produced code-stream realizes scalability by pixel-precision 
(the scalability by spatial resolution is realized by the multiplexer.) 

Depth-First and Breadth-First Quad trees 

The EQW method can be used for both lossless and lossy coding. In lossless coding, 
the primitives to be encoded in the ICUs are coefficients of a reversible wavelel 
25 transform. In lossy coding, after the wavelet transform, each transform coefficient is 
represented in a fixed-point binary format - typically with less than 16 bits - and is 
treated as an integer. 
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In each ICU, the maximum coefficient magnitude M is detennined. Then a value N is 
found which satisfies the condition 2^ <M < 2'^*' . The EQW works in a bit-plane 
manner: the initial bit-plane is set to 2". followed by 2'^'\2'^'^ ... and so on. A 
binary significance map is produced for every bit-plane by comparing coefficients in 
5 power of 2 increments. 

Figure IV.3 illustrates the EQW encoding process on a single bit-plane. Two working 
lists are used. One is called the list of significant primitives (LSP). Each entry in 
LSP corresponds to an individual primitive in the ICU. The LSP is initialized as an 
empty list. Another list is called the list of insignificant blocks (LIB). Each entry in 
10 the LIB is registered with the coordinates of the top-left comer of the block, together 
with the width and height infonmation. Each LIB entry represents an individual 
primitive of width and height equal to I. At the highest bit-plane 2^, the ICU to be 
encoded is put into the LIB. A temporary list of insignificant blocks (TLIB) is used 
for refreshing the LIB after each bit-plane coding pass is completed. 
15 There are two methods used to add the four sub-blocks into the LIB that occur as a 
result of the quad tree decomposition. The first method, known as depth-first quad 
tree coding, is to insert the four sub-blocks into LIB at the position of their parent 
block. In this case, the four child blocks are evaluated immediately after the parent 
block. This rule is applied recursively until no more subdivision is possible. This 
20 means that control returns to the next entry in the LIB only after the present entry is 
completely encoded up to the highest resolution level. The second method, known as 
breadth-first quad tree coding, is to add the four sub-blocks under consideration to the 
end of LIB. In using the breadth-first process, all parent blocks at the same level will 
be processed first, followed by their respective children blocks. 
25 After all entries in the LIB have been processed, the entries in TLIB can be reordered 
according to certain pre-defined set of rules (this is an optional step.) Experimental 
evidence suggests that a higher PSNR can be achieved by using an effective 
reordering scheme. Then the LIB is replaced with the TLIB and the coding resumes 
at the next bit-plane. 

30 The next step is the refinement pass. The N"" bit is output for entries in the LSP. 
Then, the process resumes using the new LIB and a new bit-plane set to Error! 
Objects cannot be created from editing field codes.. 
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EQW with VLC 

In order to use variable length coding (VLC), the standard EQW algorithm of the 
previous section is modified as follows. 

Step 1. Initialization: output n satisfying the inequality 2" < max{jC,. |} < 2"*' , set the 
5 list of significant primitives (LSP) as an empty list. Put the ICU to be coded into the 
LIB. 

Step 2. Bit-Plane Sorting: for each entry of the LIB, perform quad tree coding. 

If all primitives within the block are insignificant, output a 0 bit and add this block to 

the temporary LIB (TLIB). 

10 Otherwise, look at the significance of the four sub-blocks. Then according to the 
VLC table, output the corresponding VLC sequence (see the table in Figure 25). If 
significant, sub-blocks that are not a single primitive are inserted into the LIB at their 
parent position. If insignificant, they are added to the TLIB. If the sub-blocks are 
single primitives, they are added to the LSP. Sign bits are output only if units are 

15 significant. If insignificant, the units are added to the TLIB. 
Step 3. Reordering (optional). 

Step 4. Refinement: for each entry in the LSP, except those included in the last 
sorting pass, (i.e. with same n), output the nlh most significant bit of |Cj,| . 
Step 5. Quantization Update: decrement n by 1 and go to step 2. 

20 Some other researchers have also proposed the embedded coding methods where quad 
tree is used to explore the intra-band similarity However, some major 
differences exist in how the quad tree operates. The work of [2] has suggested that 
empirically the quad tree decomposition should stop at size 16x16. However (as a 
result of our experience and experimentation) this size may not always be the optimal 

25 choice. Even without VLC, the efficiency of DAC's EQW is comparable with EZW. 
In fact, for a 2x2 block, if it is insignificant, one zero can represent four zeros. 
However, when a 2x2 block is significant, quad tree coding is inefficient. 



"' F. Kossemini et al. "Embedded Quadtree Prediciive Codec", ISO/1EC/JTC1/SC29/WG1/N667. 

Taubman. D. et al. "EBCOT: Embedded Block coding with optimized tmncation", ISO/lEC/ 
JTC1/SC29/ WGl/ N1020R. 
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Using Block-Based Quantization 

The embedded quad tree (EQW. EQPC, etc.) is a special case of the more general 
block-based quantization techniques. Other block-based quantization algorithms, 
such as EBCOT. can also be very efficient for type-2 ICU coding. 

Using JBIG Algorithms 

It is noted that in all implicit quantization schemes using progressive bit-plane coding 
methods, including the embedded zero tree and the embedded quad tree approaches, 
at each bit-plane, the ICU coder is actually dealing with a bi-level image (the 
significance map.) Therefore, by adopting those efficient JBIG algorithms as ICU 
coders in a JPEG 2000 system, one can realize two modes of interplay with JBIG. In 
the following EQW is used as an example to illustrate this idea. 

In the first mode, the transcode mode, the JBIG ICU coder and EQW ICU coder 
provide a two-way transcode between JPEG 2000 and JBIG (Figure 26). In this 
transcode, the EQW procedure is called for bi-level type-2 ICUs, with the total 
number of bit-planes being set to one. Note that the sign bit coding in the EQW 
algorithm should be skipped. 

In the second mode, the embedding mode (Figure 27), the proper JBIG routines are 
called as the bit-plane coder in the EQW routine for certain bit-planes and for certain 
decomposition sizes (the quad tree decomposition stops at this particular size and the 
coding is handed over to the JBIG routine). Since the algorithms for bi-level data 
coding are also evolving with the compression technology, having an embedding 
JBIG mode reserved in the JPEG 2000 system can make the Standard evolve with its 
sister technologies. 

Coding Tvpe-3 ICUs 

Type-3 ICUs are defined to suppon various inter-band coding methods for wavelet 
transfonm coefficients. In particular, the methods of EZW, SPIHT, etc.. are known to 
be efficient approaches for coding the type-3 ICUs. In testing these approaches, it is 
noted that the standard versions must be modified slightly to make them fit into the 
RICS architecture. Normally, these algorithms implement the ICU coding and 
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multiplexing steps into an integrated procedure. A simple modification is required to 
bridge the existing architecture to the RICS architecture. The two steps must be 
separated, which is straightforward. Since the multiplexing stage in RICS can 
effectively simulate the performance of zerotree-based schemes (refer to Ch.V), this 

5 particular split in functionality has virtually no impact on the attainable coding 
efficiency for standard schemes. Instead it adds more flexibility and openness to the 
existing architecture. 

Intra-Repion Coding with Different Types of ICUs 

Generally there is a need for using different types of ICUs for a intra-region coding. 

10 For example, in coding the wavelet transform coefficients, the LL subband is usually 
encoded differently than the other high-pass subbands. Because different sets of 
primitives may possess different statistical characteristics, providing several coding 
scheme choices in the intra-region coding module may offer a higher coding 
efficiency. In practice, the following combinations are useful. 

15 • In a wavelet decomposition, coding the LL subband with type-2 ICUs whereas the 
other subbands with either type-3 or type-2 ICUs. 

• Tiling an image into several regions, with some of the regions paved with type-1 
ICUs and others with type-2 ICUs. 

• Using different coding schemes in different type-2 ICUs. 

20 Intra-Region Coding without Any ICUs 

Intra-region coding can be performed without specifying any ICUs in a particular 
region. In this case, the entire region is considered an ICU. If the natural geometric 
shape of the region fits into any of the three ICU categories, the appropriate ICU 
coder can be used. If the region has a very irregular shape, then the ID coding 

25 method can be an efficient approach. Once the data for a particular region is scanned 
in a certain pre-defined order to form a ID stream, the following ID progressive 
sorting algorithms can be used to produce an embedded code-stream. 

ID Progressive Sorting Algorithms 

The EQW algorithm can be readily extended to 1 D cases where a binary partition is 
30 used instead of a quad tree decomposition: 
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LeiL= {c,} be the ID list to be encoded. LSP the list of significant primitives, LIS 
the list of insignificant subsets, TLIS the temporary list of insignificant subsets. 

Step 1. Initialization: Output n satisfying the inequality 2" < max{|c, |} < 2"*' , set the 
LSP as an empty list. Put the set L into LIS. 

Step 2. Bit-Piane Sorting: For each entry in the LIS, perform the binary partition. 
If all primitives within the subset are insignificant, output a 0 bit and add this subset 
to the TLIS. 

Otherwise, output the two bits that reflect the significance map. For those subsets, if 
they are not single primitives, insert them into the LIS at their parent position if they 
are significant, or add them to the TLIS if they are not. If the subsets are single 
primitives, add them to the LSP and output the sign bit if they are significant, or add 
them to the TLIS if they are insignificant. 

Step 3. Reordering (optional). 

Step 4. Refinement: for each entry in the LSP. except those included in the last 
sorting pass. (i.e. with same n), output the nlh most significant bit of |c,.|; 
Step 5. Quantization Update: decrement n by 1 and go to step 2. 

Processing with a Multiplexer Architecture 

The JPEG-2000 committee has been investigating many emerging compression 
technologies in order to define a new still image compression standard. The emerging 
standard will address both present and near future image compression needs. The task 
of selecting what technologies to include in the new standard is not easy given the rate 
at which technological advances occur in this area. The primary focus of DAC has 
been to develop a still image compression engine that addresses the issues set forth by 
the standards committee. 

The concept of using a multiplexer (MUX) has been adopted in some standardized 
processes of information coding such as MPEG. In contrast, incorporating a MUX as 
an integral part of a still image compression engine is rather unique, and there has 
been very little research conducted in this area. However it will be shown that the 
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MUX concept is exiremely useful in the development of a general purpose still image 
compression engine. 

Background 

One of the concerns in developing a new Standard is the layout design of the encoded 
5 bit stream. A protocol must be developed to guide the new standardized compression 
procedures. Some of the important considerations in this area are listed below. 

• A well defined bit stream. All syntactic and semantic aspects of the encoded 
image format must be defined. The JPEG-2000 standard will encompass a wide 
variety of compression needs. One to one relationships will exist between fields 

10 that exist in the bit stream and functionality that exists in the core compression / 

decompression engine. The number of fields that exist in the bit stream is directly 
related to the overall functionality encompassed by the new standard. 

• A highly degree of data accessibility. The bit stream must be organized in such a 
manner that it need not be decoded in its entirety to interpret the physical meaning 

15 of the compressed data. Highly accessible data requires that the encoded data is 

divided into logical units each corresponding to a specific part of the original 
uncompressed information. The way this was accomplished in the existing JPEG 
standard was to partition the original image into 8x8 blocks. Each block was 
compressed using the DCT. Quantization tables were developed for coding the 

20 coefficients such that the reconstructed image quality was degraded in a standard 

fashion. However, strict block based data accessibility may not be a suitable 
medium for many of the current image processing needs. 

• Dynamic re-orderability. Dynamic re-orderability is a data access issue relating 
to the degree thai encoded information can be reorganized to suite user specific 

25 needs for an image under consideration. This type of data access is an important 

concern for many applications, especially in the medical field. 

• Built in error resilient In any transmission medium there is always the 
possibility of incurring bit errors. The encoded bit stream must have a certain 
degree of error resilience to address this concern. 

30 Using a MUX to organize the encoded bit stream can address these issues. The first 
step is to divide the encoded information into logical groups. Organizing the data in 
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ICUs leads naturally into a high degree of accessibility. Each ICU is an independent 
unit. The size and position information is self-contained within each ICU. One of the 
effects of using a muliplexed bit stream design is that it is inherently error resilient. If 
any logical unit incurs an error, only that unit needs to be recovered. If the encoded 
5 data is divided into an array of logical units in preparation for MUX encoding, 
dynamic re-ordering is much easier to achieve. 

MUX Design Overview 

There are five important considerations to address in the design of the MUX for still 
image compression systems. 

) • Defining a working model. 

• Defining an optimal data processing order. 

• Capturing the data to be compressed in a suitable data structure. 

• Defining a data priority to be used for compression. 

• Packing the compressed data based on data priority and processing order. 

A Working Model for the MUX Discussion 

In preparation for the MUX discussion, consider that a multi-resolution 
decomposition hierarchy of data has been obtained for an image to be compressed. 
Also consider that there are three channels of data. Assume that the original image 
has been transfonned from RGB to YUV for example. In addition to this, assume that 
the U and V channels have been down sampled by a factor of 4:1. This type of 
approach is common for lossy image compression. The data sets appear in Figure 18. 
The diagram illustrates a 4 level multi-resolution decomposition for the Y-channel 
and 3-leveIs for each of the UA'-channels. The wavelet transform MALLAT 
decomposition is used here for convenience. 

This is by no means the only model that exists for the MUX design. DAC makes the 
distinction between region and non-region processing modes of operation in their 
design implementation. A hierarchy of list structures is used to control optimal bit 
budget distribution and flow of both non-region and region bit levels through MUX 
channels. The discussion presented in the next Section begins by introducing the 
concepts in normal processing modes. The normal MUX modes of operation are 
generalized later on to include specialized region and mixed processing modes. 
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Generally speaking, given a wavelet multi-resolution decomposition of data, 
coefficients of a fixed magnitude contained in a lower resolution level are more 
important to the image reconstruction procedure than coefficients of a similar 
5 magnitude at a higher resolution level. This leads to a natural ordering of the 
coefficients beginning at the lowest resolution level. The natural ordering concept is 
an important consideration in the design of the MUX. 

General Color Processing Order (Lossy Case) 

Given the multi-resolution data sets of Figure V.l, a natural processing order exists 
10 for each decomposition channel. This natural order is illustrated in Figure 29. This 
particular order reflects a level priority inherent to each data set. 

Another two orders are obtained by simply interchanging orientation data sets in a 
given level. TTie difference lies only in the convention used for ordering the detail 
information (i.e. LH. HL. HH). The 4-1-1 down sampling relationship that exists in 
15 the color transform domain also exists in the wavelet decomposition data since there 
is one less transform level for the UA^ channels. This is one of the keen design 
considerations that can be exploited in many stages of the MUX implementation and 
encoding / decoding procedures. The data is organized into list structures to cover the 
natural processing orders that exist in the data to be compressed. 

20 Level Priority Processing Orders (Lossy Case) 

The individual orders of Figure V.2 can be combined into a single natural processing 
order for the example under consideration. The new order is obtained by interleaving 
the expressions of Figure 29 using the inherent decomposition level priority of the 
data. The result is illustrated in Figure 30. 

25 Notice that in this particular order, the 4-1-1 color down sampling relationship is 
maintained. In terms of the reconstruction process, the complete level 4 Y-channel 
information (i.e. LL(-HL4-LH4-HH4) is required together with the U and V channel 
low pass information (i.e. LL3 in each case). This is the information the decoder will 
need first to reconstmct the low pass inverse color transfonn image at the lowest 

30 inverse wavelet resolution level. A "level split" on the U / V data channels at the 
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lowest inverse wavelet resolution level is used to exploit the natural relationships that 
exist for both inverse transform spaces. A similar order exists for the lossless case 
where full data sets are used. 

Color Interleave Processing Order (Lossy ) 
5 Another processing order that exists for color images in the wavelet transform domain 
is obtained by interleaving the color channel detail information. The order is 
illustrated in Figure 31. This type of ordering may be suitable for certain applications 
such as those that require the data to be encoded for a progressive download. As in 
the previous case, the down sampling relationships that exist are maintained in both 
10 color / wavelet transform spaces. 

Lossless Color Processing Orders 

As in the lossy case, similar processing orders can be designed for lossless color 
image compression where full data sets exist for all color channels. In this case, 
corresponding orientation information from each wavelet channel is related in the 

15 original inverse color space. These data should be processed and grouped together in 
the wavelet transform domain and eventually the final bit stream. 
The natural processing order for lossless color image processing is illustrated in 
Figure 32. The color interleave order for lossless image processing is illustrated in 
Figure 33. Note that in each case, the full decomposition data set for each channel is 

20 maintained. The inherent relationships that exist in the color transform space are 
maintained in the wavelet transform space for ordering and reconstruction of the 
original image. 

Capturing the Data in a MUX List Structure 

There are a number of candidate algorithms under consideration by the JPEG-2000 
25 committee for optimal quantization of the wavelet coefficients. A common approach 
used to quantize the wavelet coefficients in a progressive manner is to use an optimal 
bit plane ordering technique. DAC has testing a number of bit plane ordering 
techniques in both one and two dimensional schemes. Bit plane ordering techniques 
are generally classified as using implicit quantization for the wavelet data sets. 
30 Coefficients are packed into the final bit stream based on bit level priority. The bit 
plane processing technology will be used to introduce the data structure used in the 
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MUX design. However, the MUX processing discussion that follows is not restricted 
ID bit level processing algorithms. 

The general ordering and processing relationships are useful for introducing the MUX 
control architecture at a basic level. However, the key to the operation and many 
5 benefits of the MUX technique is in the data structure design used to encompass the 
natural ordering concepts. The first step is to define a list structure for each level, 
channel, and orientation of data that exists in the wavelet transform space. A typical 
example of this type of structure is given in Figure 34. 

A general list structure for each level data set can be utilized in many ways. A 

10 lossless "prepack" of multi-resolution hierarchy information is useful for exploiting 
the many relationships that exist in the data taken in MUX list processing context. 
The data and information contained in each list is used to control how the final bit 
stream is organized. All data packed into the final bit stream must conform to this 
structure. New fields may be added, but the basic operation of the structure remains 

15 the same. Most implicit quantization techniques are implicitly decodable. In other 
words, minimal header information is required for each list. The MUX architecture 
organizes the lists or units of information into a bit stream that is scaleable in terms of 
bit precision and resolution and is controllable by the many MUX modes of operation. 
As an example, suppose that EQW (DAC's current two dimensional bit level 

20 processing design) is used to capture the leading ones and refinement bits at each bit 
level in each data set. Following the order given in Figure V.3. a multi-dimensional 
hierarchy of list structures is defined for use in the ensuing data processing stages. 
Within each list (as EQW progresses), the fields in the corresponding list structure are 
updated. The bit stream at each bit level is appended to the MUX "prepack" buffer as 

25 well as its corresponding size being added to the total in the bit packing information 
field. The bit level position where processing begins is put in the high bit information 
field. The total packed list size is placed in the total packed field. Finally if more 
than one internal processing scheme is used (e.g. NULL transform mode for bi-level 
information or other intemal modes of operation), the scheme used for the list is 

30 placed in the scheme field. 
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Each MUX list is fully independent, implicitly contains a full set of statistical 
information available for determining the amount of data to pack for each list and is 
optimal for organizing the final compressed bit stream. 



5 Data Priority and MUX Control 

The information contained in each list sinjciure is used to organize the final bit stream 
based on end user compression requirements. DAC has implemented numerous MUX 
modes o operation. Two common modes of operation are outlined in this Section. 

Signal to Noise Ratio (SNR) Progressive Mode 

10 This mode of MUX places priority on the data such that the final ordering optimizes 
the SNR given a user specified image compression size. For the discussion that 
follows, assume that the user has specified the desired compression ratio in some 
manner. Assume that an optimal bit budget distribution scheme is in place. Thus the 
number of bytes that has been allocated to each color channel is known (more will be 

15 said on how to determine the optimal bit budget distribution for each channel using 
the MUX later). The information fields in the MUX list structures are used to 
determine the amount of data to pack for each list such that the SNR is optimized for 
the specified bit budget. A psuedocode explanation of this technique appears in 
Figure 35. 

20 Fields that refer to the MUX list structure appear in italics while other local variables 
appear in nonmal typeface. Basically the idea is to loop through each list checking to 
see whether the current processing bit level is equal to the bit level of the list under 
consideration. If the two bit levels are equal, the channel bit budget is decreased by 
the amount of data available in this particular list for the bit level in question (from 

25 MUX field pliPackingInfo[cCurBitLeveI]. The MUX field liCurBytesCount is 
incremented by the number of bytes available. Additional processing takes care of the 
remaining bits. The remaining bits will be considered in the channel bit budget on the 
next visit to this paaicular list. The MUX field cCurBiiLevel is decreased by one 
such that it will again be enabled when the bit level in the main loop is decreased by 
30 one. 

This procedure optimizes the SNR since it processes the largest bit levels in each list 
first. Thus information for the largest coefficients throughout individual 
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decomposiiion data sets is sent first according the natural processing order outlined 
earlier. Once the channel loop exits, each MUX list will contain a field indicating the 
amount of data to pack into the final bit stream in each case. 

Resolution Progressive Mode 
5 The resolution progressive mode of bit budget distribution is a simple extension of the 
psuedocode implementation of the SNR mode of the previous Section. In order to 
obtain a resolution progressive mode, the 'Waveket Transform level' for loop is taken 
outside of the mzm BitBudget / BitPlane while loop such that resolution level is given 
priority. In this manner, a lower resolution image can be reconstructed in the inverse 
10 wavelet transform stage. Tht resolution of the ensuing image depends upon the 
number of levels required in the end user requirements. 

Data Packing Using the MUX Lists 

Composing the final bit stream is a simple matter once the bit budget distribution 
scheme has run to completion. The packing technique can follow one of the normal 

15 processing orders outlined earlier (e.g. color level priority or level priority color 
interleave) or another access technique as required by the end user. The amount of 
data to pack for each list has been detemiined in the bit budget distribution stage. If 
the data to be packed for a given list is greater than zero, then the total data size, 
scheme and high processing bit level are packed as header information into the final 

20 bit stream. If the data to be packed for a given list is zero (i.e. not required), then a 
smaller header is packed to indicate a zero length list. Currently a register type 
approach is being developed to handle zero length lists. All lists that have data 
available will set a bit position in the list register so that the decoder can determine 
which lists have made a contribution. Empty list will be flagged with a 0 bit in the list 

25 register. The list register is packed in the final code-stream. This technique will save 
the size of the list header overhead of empty lists. 

Overhead Consideration of the MUX Architecture 

There is a small overhead placed the compressed data by the MUX implementation. 
There are three header fields to pack for each MUX list. 

30 • The total bytes packed. 

• The scheme used to process the list. 
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• The highest bit plane for data in list. 

Note that if there is only one internal scheme used to process the entire image, then 
the scheme field does not need to be packed for each list. 

In order to calculate the cost of using a MUX in terms of header sizes for each list, 
5 consider that a 24 bit color image of size 2048 by 2048 needs to be processed in a 
lossy fashion (i.e. the working MUX model of Figure 28). Assume that the size of the 
smallest wavelet transform level is 8 by 8. Thus there are 8 decomposition levels to 
consider in forming the MUX list hierarchy. Assume that the wavelet decomposition 
data has been converted to a 16 bit integer representation. Also assume that there is 
10 only one internal processing scheme for the lists. Given these initial conditions, a 
worst case analysis is conducted to determine the header size for each list. The results 
are tabulated in Table V.l for the Y-Channel wavelet decomposition. 

From table in Figure 36. the total Y-channel header size for this particular image is 
443 bits. A similar calculation can be conducted on the UA^-Channels (assuming one 

15 less decomposition level for each) to yield 368 bits each. The original raw image size 
is 2048 X 2048 x 3 bytes. Thus the MUX packing overhead in this particular example 
is about 1 header fait for every 854 bits of compressed data. In the lossless case where 
each channels has the same number of decomposition levels (or correspondingly, in 
the 8 bit grayscale case), the overhead is still very small at 1 header bit for every 757 

20 bits of compressed data packed. 

The table in Figure 37 outlines the overhead size for square images with dimensions 
that are a power of 2 beginning at 16 by 16 and ending at 64k by 64k. Note that the 
overhead size if 44 bits in both the lossy and gray scale cases. The reason for this is 
that the wavelet transform is not used on the UA' channels when the original image 
25 dimensions are 16 by 16. The Y-channels is decomposed once. A plot of percentage 
overhead versus image dimension is given in Figure 38. Both lossy and lossless cases 
appear on the same plot. Note that the image size can be quite small while still 
maintaining a relatively low overhead in terms of the total MUX list header sizes. 
Bit Budget control from the MUX Architecture 
30 One of the most difficult tasks that must be addressed in any compression scheme is 
scalability. The intemal procedures must be developed such that the image under 
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consideration can be compressed to a user specified size in an optimal manner. The 
task is compounded for color images since there are three channels to consider. The 
major problem is how to distribute the user specified bit budget between the three 
channels. 

5 Standard color transforms such as YUV or YIQ are normally used to redistribute the 
RGB color information prior to the multi-resolution decomposition stage of the 
compression procedure. These transformations concentrate most of the energy into 
one channel making them better suited for compression. Another color transform that 
is gaining popularity is the KL transform. This particular transform technique is 

10 based upon a principle component analysis (PCA) implemented on the original raw 
color information in order to determine the optimal color redistribution for the three 
RGB channels. Generally speaking, after completing a color transform stage, two of 
the three channels are down sampled by a factor of 4 to 1 for lossy compression. In 
doing this, the amount of data that must be compressed is reduced by a factor of 2 

15 with minimal loss in visual quality for the reconstructed image. 

The easiest bit budget distribution scheme to implement is one that follows the color 
transform down sampling ratios. If the transform under consideration is YUV-411, 
then 4 parts are allocated to the Y-channel for every 1 part allocated to the U / V- 
channels respectively. However, the energy distribution of the wavelet transform data 

20 generally does not follow this strict down sampling guideline. Instead it varies from 
one image to the next. 

DAC has developed a dynamic bit budget control architecture that is based on the 
implicit information contained in each MUX list. In this Section, a bit budget 
allocation procedure will be introduced that can be used for optimal scalability in 

25 normal processing modes of operation. This technique is based implicitly on the 
amount of energy contained in each color channel of the wavelet decomposition data 
sets. A data ratio concept is used together with the prepack information contained in 
the MUX list structures to determine the budget for each channel. 
Compression procedures that process wavelet coefficients in a bit plane order are 

30 implicitly tracking the energy contained in the decomposition. The most significant 
bits (MSBs) of the largest coefficients are packed first followed by refinement bits 
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and the MSBs lhat are significant for coefficients at the next bit level. This process is 
repeated in a recursive fashion until the desired compression size is obtained. 
Let yij , Uij and Vy . be the total data available in a particular orientation level for the 
three color channels. Then the total amount of data available in the Y-channel 
5 decomposition Y, is obtained by summing the individual totals. A similar procedure 
can be followed to obtain U, and V,. 

N 3 

= 5, J "0 (Eq.V.l) 



The next step is to calculate the pack ratios to be used for each channel of the wavelet 
decomposition hierarchy. The pack ratios are determined by taking the ratio of the 
two largest amounts of data to the smallest amount of data. Let the three pack ratios 
be denoted Ry , Ru and R,. The smallest data size will be in either the U or V 
25 channels because of the down sampling step. Note that one of the pack ratios will be 
unity. 



' MIN(U„V,) 
U. 

^" " M1N(U..V.) ^E'l-'^-^) 
V, 

^' " MIN(U,.V,) 



The pack ratios given in Eq.V.2 are used to determine the optimal amount of data to 
allocate for each color channel based on the user specified compressed file size. 
However, any additional overhead introduced into the final bit stream by the MUX 
must be taken into consideration. Let yhij, uhjj and vhij be the individual header sizes 
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in each wavelet channel for a particular orientation level. Then the total header sizes 
for each wavelet channel are obtained by summing the individual totals. 

N 3 
N 3 

,0 V, = uK. ^ -^^^^K (Eq.V.3) 



If doing ROI processing, the packing overhead introduced by the mask must be taken 
into consideration in determining the bit budget for each channel. Lei the total pack 
20 size for the ROI mask be denoted as M,. The logical approach to take is to make an 
adjustment in the bit budget for each channel based upon the pack ratios of Eq.V.2. 
In this manner, the mask overhead is distributed proportionally to each wavelet 
channel. The unit adjustment factor V, is determined from the total mask overhead 
and the total pack ratio. 

U. = (Eq.V.4) 

Rv + Ru + Rv 



The adjustment size for each channel is determined by using the result of Eq.V.4 and 
the pack ratios of Eq.V.2. 

Y- = RvU. 

U. = Ru . u. (Eq.V.5) 



45 A useful set of calculations that comes out of Eq.V.l to Eq.V.5 is the minimum 
compression bits per pixel (BppMin) and the maximum compression bits per pixel 
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(BppMax). These values are calculated below (assume that the original raw image file 
size is Ft bits). 



)PM.. =— • (Y, + U.+ V,+ Yh + Uh + V»+Y. + U.-f V.) 
jpMin =— • (Yh + Uh + Vh + Y.+ U. + V.) 



(Eq.V.6) 



These expression represent bounding compression values available for the image 
under consideration. If all data is packed for each of the wavelet data sets, the result 
15 is BppMax. If just the header infonmation is packed, then the result is BppMin- 

Suppose the user specifies a compressed file size of Fu„r. The unit bit budget Ubb is 
determined from the pack ratios and the total header sizes. 

Uea (Eq.V.7) 
Ry + Ru + Rv 



25 Now the optimal bit budget for each channel Ybb, Ubb and Vbb can be determined 
using Ubb, the individual pack ratios and the overhead parameters. 

30 Ybb = Ubb • Rv + Yt- Y, 

Ubb = Uu,- Ru +Ui,-U. (Eq.V.8) 
Vbb = Vbb. Rv + Vh- V. 

35 

The expressions given in Eq.V.8 will always yield a near optimal bit budget for each 
color decomposition channel. The technique is based implicitly on the energy 
distribution map recorded in MUX lists and the total overhead associated with the 
underlying process. Note that if processing with no ROI. then that variable falls out 
40 of the calculation. 

This concept can be extended to form pack ratios for each subband if more accuracy 
is required. A similar data ratio technique can be implemented to control the bit 
budget in that case. The necessary information is contained in the MUX list 
structures. However for most practical application, the cited distribution technique is 

54 



Copied from 10340491 on 04/01/2005 

3NSCX3CID: <VilO 0O«9571A2.l.> 



wo 00/49571 PCT/CAOO/00134 
sufficiently accurate. The user specified file size can be obtained within several bytes 
. in implementing this technique. In addition, the maximum and minimum attainable 
file sizes are known prior to the final packing stage. 

A simple routine has been developed for distributing the final bit budget in the level 
5 where it is determined it will expire. The idea is quite simple. Given that the bit 
budget will expire on a particular bit plane in a certain transform level, then the bit 
budget is redistributed such that each orientation gets a proportional amount based on 
the amount of data that each orientation can take. The amount of dau each 
orientation gets is determined by using ratios between orientations on the bit level 
10 under consideration. This is an important consideration since the net effect is 
approximately a 1 dB improvement in the PSNR. 

Using the MUX for Region Processing 

So far the focus has been on the design of a MUX control architecture for normal 
image processing modes of operation. The main highlights are the data processing 
15 orders, the MUX list structure, the bit budget control technique, and organizing the 
final bit stream. DAC has extended the MUX concept to include a variety region 
processing modes of operation. 

Region Processing Orders 

The processing orders introduced in the previous Section are specific to the normal 
20 processing modes where one or more ROl are not specified. DAC has implemented a 
region processing design based on the natural extension of the MUX concepts 
introduced in the previous Sections. Both lossy and lossless processing orders are 
again considered. 

General Color Region Level Processing Order 

25 The region level processing order places priority on each ROI in a descending fashion 
followed by the normal level priority scheme of the original MUX design. Inherent in 
the MUX scheme is the general assumption that data of similar magnitude in each 
succeeding region is less important to the overall image reconstruction. A secondary 
priority key is placed on the wavelet transform level. Data of similar magnitude 

30 contained in a lower resolution level is more important to the reconstruction 
procedure than data at a higher resolution level. The general processing order for 
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each color channel (assuming a 4 region priority scheme) is given in Figure 39 where 
R; is the region level for i = 1,2.3,4. 



Repion Level Color Processing Order (Lossy) 

The individual orders of Figure V.IO can be combined into a single processing order. 
5 The new order is obtained by interleaving the expressions making use of the inverse 
color and wavelet transformation relationship that is presented for the normal MUX 
modes of operation above. The result is illustrated in Figure 40. 

The intrinsic 4-1-1 color down sampling relationship is maintained for region 
processing MUX modes. In terms of the reconstruction process, the complete level 4 
10 Y-channel information (i.e. LLi-HL^-LH^-HHJ is required together with the U and V 
channel low pass information (i.e. LL3 in each case) for each region under 
consideration. 

Color Interleave Region Level Processing Order (Lossy) 

The color interleave processing order extends naturally to the MUX region processing 
15 modes. The order is illustrated in Figure 41. This type of ordering may be more 
suitable for certain applications that require the data to be encoded for a progressive 
download based on regions of imponance. As in the previous case, the down 
sampling relationships that exist for the inverse color and the wavelet space are 
maintained. 

20 Region Color Processing Orders (Lossless) 

As in the lossy color region processing MUX modes, similar processing orders can be 
outlined for the lossless case where full data sets exist for all color channels in each 
region. In this case, conesponding decomposition level orientation information in 
each channel is related in the original color transform space for each region under 

25 consideration. These data are processed and grouped together in the wavelet 
transform space and are ultimately put in the final bit stream by the MUX control 
architecture. 

The processing order for lossless color region processing is illustrated in Figure 42. 
Note that in each case, the full decomposition data set for each channel is maintained. 



56 



Copied from 10340491 on 04/01/2005 

WSDOCrD: <WO_004SS71A2.I_> 



wo 00/4957] PCT/CAOO/00134 
The inherent relationships that exist in the inverse color and wavelet space are 
maintained for each region channel. 

DAC has developed both SNR progressive and resolution progressive MUX modes of 
operation for region processing. In addition, several specialized MUX modes were 
5 developed that may be useful in certain applications. The regions can be defined 
automatically by the technique described in Chapter III, or user defined ROI can be 
used to group the data. User defined ROI can have arbitrary shape or can be defined 
in terms of simple geometric elliptical or rectangular region primitives. 
One of the specialized MUX region processing modes can be introduced at this point 
10 since it fits into the context of the ordering discussion. 

Transparent Region Color Processing Orders (Lossy) 

This particular region color processing order is useful for applications that require 
transparent region channels that has little or no influence on the ordering of the data. 
Either of the normal SNR and resolution progressive modes are overlaid with a 

15 complete region channel description. The technique is implemented simply by taking 
the region index to the inner processing loop in the pseudo code implementation 
example of Figure V.8. This ordering technique may be useful in video processing 
applications where a complete region mask description is required together with a 
compressed frame of information. The region description may be used to process 

20 subsequent frames in the sequence. One of the processing orders that falls into this 
MUX mode of operation is illustrated in Figure 43. In this case, the high bit plane of 
each region list is processed first causing the region classification of all coefficients to 
be transparent to the distribution and packing routines. 

Transparent Region Color Processing Orders (Lossless) 
25 Corresponding transparent region processing modes exist for the lossless case. A 
typical example is given in Figure 44. As in the lossy case, the final bit stream is 
organized independent of the constraints of the region channels which may be useful 
in certain cases. 
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MUX List Structure for Region Processing 

The MUX organizational list structure concepts introduced for normal processing 
modes have been extended to include many region processing modes of operation. In 
this case, the total number of lists is a function of the number of ROI. If there are 4 
5 ROI, then there are 4 times as many MUX lists. However, the basic operation of the 
MUX control architecture is similar in each case. 

There are three high level modes used to categorize the MUX architecture. 

• Normal processing mode (ROI disabled). 

• Region processing mode (ROI only). 

10 • Mixed processing mode (normal and ROI enabled). 

The next Section outlines the operation of the bit budget and MUX controls for many 
of the region processing modes developed at DAC. The mixed processing modes are 
briefly discussed later in the Chapter. 

Bit Budget Control for Region Processing MUX Modes 
15 Transparent Region Level Color Mode 

Transparent region processing is mentioned above. In that particular ordering 
example, the region index is placed in the inner most loop in the bit budget 
distribution function. The processing index placement is basically a method of 
incorporating a transparent region layer of processing into the normal MUX modes 
20 outlined earlier in the Chapter. Processing and packing in this mode gives ROI a low 
priority. 

Distributing in this context ensures that the SNR progressive and resolution 
progressive MUX modes operate as they did before. The exception in this case is that 
there is a variable length region mask overhead that must be taken into consideration. 
25 The bit budget distribution functions are calculated as they were for the normal MUX 
processing modes. Based on the file size (or resolution level) requirements, each 
color channel receives a proportional amount of the bit budget based on the ratio 
distribution technique used in the MUX architecture. 

The DCT auto regions and arbitrary or primitives user defined region types can be 
30 operated in this mode. In DAC's internal processing architecture the user must select 
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the region coverage technique used to categorize the wavelet coefficients. DAC has 
implemented the wavelet mask down sampling technique of the JPEG-2000 VM, and 
it can be used in any of the region processing MUX modes. In this manner a common 
mask can be used in each orientation level or the VM mask down sampling technique 
5 can be used for individual orientation level mask coverage, for both automatic or user 
defined ROI in lossy and lossless MUX modes of operation. The number of region 
categories can be selected as 2, 3, or 4 channels. 

In order to test the transparent mode of operation the YUV-41 1 color transform is 
applied to a 24 bit 256 x 256 Glacier park mountain scenery image. The 9-7 kerne! 
implemented in a lifting scheme is used as the wavelet transform. The DCT region 
formation technique of Chapter III is employed to generate the common masks for 
each wavelet transform level. The mean square error (MSE) is measured for a 
number of compression ratios. The result is illustrated in Figure 45 in a plot of MSE 
versus Bpp. The plot shows the break points in MSE for each automatic ROI used in 
the example. The rate of image degradation (i.e. MSE and PSNR) must be controlled 
precisely for each ROI. 

The MUX overhead is calculated based on the mask type and process selections. For 
arbitrary user constructed masks, the mask file can be loaded to guide the MUX. In 
the 4 ROI common mask case, the overhead is 0.5 Bpp in packing the entire mask. In 
the VM mask case, the overhead is 2.0 Bpp. If simple primitives are used to form the 
user mask instead of an arbitrary shape, the overhead is greatly reduced. The arbitrary 
mask overhead can be reduced with the addition of an entropy coding stage. The 
common DCT mask overhead is cited in Chapter III as approximately 1Kb (0.2 Bpp) 
for an 8 bit 256 x 256 image size. These are rough estimates that do not take the 
header sizes into account. However, they do serve as a comparative guideline for the 
current discussion. 

Region Priority Level Color Mode 

In this mode of MUX operation, the pack ratios for each wavelet channel are 
determined as in the normal processing mode. However, there is one important 
difference. In this particular mode of operation pack ratios are determined for each 
region channel. This implies that there are 12 pack ratios for a 4 ROI process. The 
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bit budget distribution is based on the overall wavelet channel pack ratios and the ROI 
data totals in each region channel. The distribution function also takes the mask and 
list header size information into account in the overall calculation. 
This mode of MUX operation is designed to distribute the MSE and the PSNR image 
5 reconstruction measurements in an approximately uniform manner for all region 
channels. A proportional amount of overall bit budget is allocated to each region 
channel based on the region and color channel pack ratios. Figure 46 illustrates the 
result of using this mode of operation for the Glacier park image. Notice how the 
quality of each region channel degrades in comparison to the others. Auto detected 
10 DCT common masks are used to generate the result. 

Absolute Region Priority Level Color Mode 

In this mode of operation, absolute priority is given to the region channels. The 
distribution function is the same in this case in that the bit budget is divided into 12 

15 according to total data ratio technique mentioned earlier. However the bit budget is 
distributed in a biased fashion giving the highest priority to the most important ROI. 
Some regions may not receive a budget based on the compression requirements. The 
bit budget is distributed region by region beginning at the most important ROI. Thus 
complete sections of the original image can be eliminated altogether if the bit budget 

20 is small. If only the important regions of an image need to be saved, this technique 
can be employed to partition the image. 

A plot of MSE versus Bpp for the absolute region priority color MUX mode is given 
in Figure 47 for the Glacier park image. Notice how the regions fade out very fast 
and sharply. This occurs when the region no longer has any budget allocated to it. At 
25 that point, the region is basically invalidated in ternis of any contribution to the 
reconstructed image. That portion of the image basically fades out. Note that the 
amount of quality of the fade out depends on the down sampling technique to translate 
the masks. The VM mask formation technique causes a very gradual fade out to 
occur. In the common mask approach, the effect is much more abrupt. 

30 Scaled Region Priority Level Color Mode 

In this mode of MUX operation degradation rate of each ROI can be controlled. Each 
ROI contributes to the final bit stream. However, instead of having each ROI degrade 
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at a uniform rale (as Section V.5.3.2.), the quality measurements of each succeeding 
can be controlled. The initial bit budget for each region channel is calculated as 
before in that the specified compressed file size is split into 12 bit budgets to be 
spread between each of the 4 region channels. However after the initial split, the 
5 allocation amounts are changed heuristically based on a priority factor that is set for 
each ROI. For example, suppose it is decided that ROI 1 is 50% more important than 
ROl 2. ROI 2 is 30% more important than ROI 3 and ROI 3 is 20% more important 
than ROI 4 (the background). Using this assumption as a starling point, the amount of 
data allocated to each region channel is changed slightly. The net effect is to cause 
10 the quality measurements in each region channel to degrade at slightly different rates. 
A plot illustrating this effect is given in Figure 48. Note how the MSE break points 
for regions 2 and 3 have shifted slightly towards the left. 

The result of using this particular MUX mode illustrates an important property that 
should be available for any ROI processing technique. The problem with many ROI 
15 processing techniques is that it is difficult to control how much each region 
contributes to the final bit stream. The technique outlined here can be used not only 
to implement this effect, but also to control the degradation for each region channel 
based solely on an imponance factor associated with each ROI. 

Region Percentage Priority Level Color Mode 

20 In this mode of operation the user can specify a data percentage for each region based 
on the total amount available in each region channel. The distribution function is 
slightly different in this case. The data totals are detenmined for each region channel 
along with the wavelet color channel totals. There are still 3 bit budget for each 
region channel. In this case however, the amount of data to include for each region 

25 channel can be set by the user as a percentage of the total in each case. 

Currently when this mode of operation is invoked in DAC's compression engine, the 
user can cycle through the operation until the desired size for each region is obtained. 
The data totals are displayed after each run. Note also that in using this mode of 
operation, one or more ROI can be eliminated or faded heavily as in the absolute 

30 region priority MUX mode of Section V.5.3.3. 
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A plot illustrating the result of using this mode of MUX operation is given in Figure 
49. Initially all available data is included for each region followed by equal 
decrements for each region channel. 

Usinp the MUX for Mixed Processing 

5 So far both normal (non-region) and ROI processing modes have been discussed. In 
addition to these modes of MUX operation. DAC has developed mixed mode 
processing capabilities. The modes of operation discussed for both normal and ROI 
processing are extended such that they can be run simultaneously for an image under 
consideration. A wavelet transform level partition is conducted based on the desired 

I number of region and the number of non-region levels. 

Currently non-region levels can be defined for processing lower resolution levels and 
region levels can be defined for process higher resolution levels. However, the 
implementation is not restricted by this distinction. It can be changed to regions over 
non-regions or to an interlaced combination of the two. The parameter that controls 
this distinction is termed the region start level. It can be set to any valid wavelet 
transform level (it is set internally to -1 to disable region processing.) Thus the MUX 
technology can be used in exclusive non-region mode, exclusive region mode or a 
combination mixed mode of operation. 

The ordering techniques outlined for normal and region processing MUX modes will 
not be duplicated in this Section. The same concepts apply for mixed mode 
processing. The bit budget distribution functions work as they did before, but in this 
case they exist simultaneously. In some modes of operation there may be as many as 
15 bit budget definitions used to organize the final bit stream. The method used to 
determine them has not changes. The same ratio technique that has its basis in the 
MUX list structure is used to determine the appropriate allocation in each case. 
There is one additional feature that can be exploited for mixed mode processing. An 
imponance factor can be attached to the non-region levels. The net effect of this 
parameter is to taper the amount of data included for the non-region levels. In the 
current implementation, the non-region levels are processed first. Thus they 
considered first by the distribution function. The importance factor allows the user to 
decrease the amount of data included for the non-region levels by a certain 
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percentage. The delta amount is considered in the bit budget distribution for the 
region levels. 

The region processing overhead is slightly smaller in this case. Depending on the 
region start level parameter, there will be less region header iniformation required in 
5 the bit stream since there are fewer region levels. However, the header overhead for 
the lower resolution levels is quite small anyway. 

One example is given here to illustrate the operation of mixed mode processing (see 
Figure 50). In this particular case SNR progressive mode is used for the upper 3 
levels and scaled region priority on the bottom 3 levels. Thus the bottom level 
10 degradation rates have been adjusted to favor the most important region and 
attenuating in some fashion between the other 3 regions. Notice how the MSE is 
lower in this case with the same region overhead as the previous case of Figure 48. 

DCT Region Formation as a Classification Scheme 

In observing this result, it is apparent that for the same region overhead, a better result 

15 is obtained in mixed mode. There are a number of reasons for this. The first is that 
the region channel coverage is not an exact overlay. All masks used to group or 
categorize data must deal with the down sampling issue at different resolution levels, 
and the tight overhead restraints of the compression channel. As it appears in the 
results presented here the DCT region detection mode can be applied to threshold 

20 wavelet coefficients to form region channels. The amount of data packed for each 
channel can be controlled by the MUX. In addition the initial data split used to create 
the partition can be controlled before the DCT mask procedure begins. 
There is some region migration or intermingling of the coefficients that occurs at the 
boundaries where sharp changes occur in the original image. Some work was 

25 conducted in developing heuristic techniques to decrease the miscoverage in hot 
areas. And there is a benefit there. Furthermore the DCT low pass filtering stage 
affects the original priority bias used to partition the data. There will be a benefit in 
determining the optimal data split. By adjusting the original region partition (e.g. 
equally spaced regions), the plot of MSE versus Bpp can be set to partition the data in 

30 other ways. Inter-resolution coverage is another problem to consider. The mask 
generation technique of the JPEG-2000 VM addresses the miscoverage problem. 
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DAC has completed some initial investigations of the VM down sampling technique 
and further testing is required. 

There is a trade off between accuracy of each region channel and mask overhead 
required for region channel operation and generation of the accuracy in the first place. 
5 The DCT approach shows much promise especially on the highest resolution level 
where the masks are the most accurate. One of the benefits of using this the DCT 
approach is that it can be translated to other transform levels in the frequency domain. 
Or aliematively, the new masks could be generated at a different wavelet level. 

Other Processing Modes 

10 There are many other modes of region processing operation that have not been 
presented in this document. The whole user defined region generation and MUX 
modes were not included. The processing modes introduced in this Chapter can be 
used in the same way to control each region channel. Simple primitives are cheap in 
terms of overhead. One mode of operation supports full mask sets so that arbitrary 
15 mask can be loaded and sent with the data to form the region channels. The mask 
coverage technique can be selected as common masks (a smaller raw size that is up 
sampled at the decoder side) or VM masks, which retain the original image 
dimension. More time is required to study the effects of both techniques as well as 
other region fonnation schemes. 

20 Experimental results for the resolution progressive modes of operation are not 
presented in this document. These modes of operation are not currently available. 
However their implementation is not difficult. These MUX modes are currently 
under development. 

The experimental results presented in this Chapter were obtained using DACs ID bit 
25 level sorting implementation. DAC is currently incorporating the 2D version based 
on EQW into the region processing channel. There may be benefits in retaining both 
techniques. For example, in bi-level image processing. Cun^ently the EQW sorting is 
implemented for normal (non-region) processing MUX modes. Both sorting modes 
are available for lossless compression. There are numerous wavelet kernels and color 
30 transforms for both lossy and lossless compression. In addition, the number of region 
channels can be set with a current maximum of 4. Primitives can be used as desired. 
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The number of primitives is not restricted to 4 with region overlaps given to the most 
important region. 



Bit Stream Syntax 

5 DAC's current bit stream structure is rather dynamic given the different types of 
organizational strategies that exist in the underlying core technology. Normal, region, 
and mixed processing modes in addition to arbitrary wavelet, color, entropy coding 
various sorting stages in both lossy and losses cases. 
Currently the core architecture can be divided into 2 categories. 
10 • Lossless compression 

Full data sets lossless color selection, lossless integer lifting scheme wavelet 
types, lossless sorting and packing stages, with additional entropy coding can be 
selected for completely lossless and slightly lossy (some of the color transforms 
are slightly lossy) image compression. 
15 • Lossy compression 

In the lossy case, color selections, more wavelet kernels / lifting schemes, lossy 
sorting and packing stages, with variable length coding can be selected 

In both lossy and lossless cases region and mixed modes can be used. Currently 
normal operational modes can be realized in an almost transparent fashion for many 
20 current compression schemes as well as our own intemal EQW / 1 D sorting. The 
general form for the lossy and lossless header structures is given in the tables in 
Figures 52 and 53. 

The current architecture is very flexible for introducing new technologies. Many of 
the modules are completely interchangeable. The exact bit stream syntax depends 

25 largely on the mode of operation selected by the user (i.e. regions, no regions, sorting, 
lossy/lossless etc.) Generally normal processing modes can be selected for all 
transform levels, or for any number of lower resolution transform levels. Region 
levels can be defined for ail transform levels, or the higher resolution levels used in 
combination with the lower resolution normal levels. At some future date, the region 

30 and the normal levels can be mixed. 
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Each unit list organized by tiie MUX has a header tag. This header tag carries the list 
size and the high bit plane processing level for the list. The current implementation 
uses 5 bytes per list. However bit packing is currently under implementation for tag 
headers. This will reduce this by a significant amount. Only 2 packing schemes are 
used at the global level so that no additional tag header information is currently 
required. A diagram illustrating the structure of the tag header is given in Figure 53. 
As other core technologies are added to the core engine or the core technology 
advances, other fields may be required in the tag headers. 

The basic structure for normal modes of operation is given in Figure 54. The diagram 
illustrates a lossy pack arrangement for a color image (YUV down sampled color 
transform assumed) according to the packing/processing order given in Figure 30 
above. In this case the code stream consists of the file header, followed by the header 
tags/data. 

The basic structure for region processing modes of operation is given in Figure 55. 
The diagram illustrates a lossy pack arrangement for a color image (YUV down 
sampled color transform assumed) according to the packing/processing order given in 
Figure 40 above. In this case the code stream consists of the file header, region 
header, regions description and finally the header tags/data. 

The basic structure for mixed processing modes of operation is given in Figure 56. In 
this case the first item in the code stream is- the file header followed by the normal 
processing mode header tags/data. Next in line is the region header, region 
description followed by the header tags/data. Notice that there is a region start level 
parameter "r* in this case. Since there is one less transform level for UA^ channels in 
the lossy case, there is a processing shift in the UA^ channels. This makes reference 
to the "level split" notion discussed in Chapter V whereby the natural 4:1:1 
relationship that exists for inverse color/wavelet transforms can be maintained for 
internal processing and the final code stream. 
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Transcodability 

In order to illustrate the generality and openness of RICS, two examples are shown to 
depict how the 'transcode' with bi-level images (JBIG) and JPEG can be handled in 
5 RICS. 

JPEG transcode 

Figure 57 shows how a DCT-based code can be produced in RICS. With this 
transcodability. it will be straightforward to write a small conversion program to 
transcode an old JPEG file into a JPEG 2000 file. Although it is possible to transcode 
10 a JPEG 2000 file into a JPEG file (losing some scalability), it will probably be of very 
little use. 
VJBIG transcode 

Figure 58 shows the execution path illustrating how a bi-level image such as a text 
document can be efficiently handled in the RICS system. The NULL transform is 

15 applied (nothing has to be done). Generally binary text documents contain mostly 
high frequency energy. Multi-resolution decompositions will not necessarily be a 
suitable basis for efficient coding. In fact, current well-known techniques specialized 
for binary images are applied directly in spatial domain. While using a RICS system 
for processing a compound document with mixed grayscale/color/text information, a 

20 set of rectangle shapes suffice to enclose the text regions. This is roughly equivalent 
to run-length coding in existing binary image compression techniques to skip strings 
of zeros. The pixels within each rectangular region are then coded into a one- 
dimensional stream via 1-D algorithms or JBIG routines, as discussed above. 

25 Post Processing 

For reconstructed images with wavelet based coding methods at low bit rate, de- 
ringing is usually a useful post processing procedure to improve (mainly) the visual 
quality. Nonadaptive procedures, which apply the de-ringing filtering to all pixels 
without discrimination, have the following problems. 
30 • While they can successfully remove ringing artifacts, they may also wash out 
details of image. 
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• They usually rely on too many parameters that are to be determined by human 
users. 

• The process is time consuming. 

The RICS system employs an adaptive de-ringing algorithm for post processing. 
Since the ringing artifacts usually appear around edges where sharp changes occur, 
this adaptive filtering process is applied only to edge areas. As the result, the post 
processing removes artifacts around edges and prevents fine details from being 
smoothed out by the filter. Compared with non-adaptive methods, the processing time 
is remarkably reduced since the filtering is applied selectively to pixels. The diagram 
of Figure 59 shows the algorithm of the adaptive post processing filter. First of all, the 
pixels of reconstructed image are classified into edge pixels and non-edge pixels. 
Edge pixels are enlarged into edge regions since the ringing artifacts do not occur 
right at the edges but around the edges. Then the artifact removal filter is applied 
only to the edge areas. 

We did some tests on this adaptive de-ringing procedure. For the filtering stage we 
used the filter described in the JPEG 2000 VM 3.0 (B). Figure 60 shows the result of 
the modified post processing. Artifacts are clearly visible in the first picture especially 
behind the cameraman's back and around the tripods. As we can see from the second 
picture, artifacts are removed, however, the details in the grass field and on the man's 
pants are also removed. In the last picture, the modified filter successfully removes 
artifacts and at the same time retains the details. Figure 61 shows the edge area (in 
white color) that was used in producing the third picture. 

As the result of using this adaptive filtering, processing time is greatly reduced since 
less pixel are processed. The table in Figure 62 shows the performance comparison of 
the VM 3.0 (B) post filtering and the RICS adaptive filtering 

Edge Detection 

To determine if a pixel is an edge pixel, the average power difference between the 
pixel and its neighboring pixels are measured against a threshold (bottom threshold). 
If the average power difference of a pixel is above this threshold, the pixel is marked 
as an edge pixel. The threshold acts as a high pass filter that filters out the pixels with 
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low power levels (Band pass fillers can also be used for edge detection by introducing 
a proper ceiling threshold). 

Edges are thickened to form edge areas after the edge pixels are detected. With the 
right choice of the threshold, this modified filtering process can be applied to 
5 reconstructed images with compression ratios as low as 8:1. Since for most wavelet 
coded images, no artifacts can be visually detected for compression ratio lower than 
8:1, our results suggest that this filter can be applied to most reconstructed pictures 
regardless of its compression ratio. 

Parameters Optimization 

One of the problems with Shen's filter is that there are many parameters associated 
with the filter, which makes the filter difficult to use. Reducing the number of 
parameters will improve the usability of the post processing filter. 
There are three different potential functions implemented in VM 3.0 (B) for 
controlling the post processing: Quadratic Truncation, Huber, and Lorenzian. Our test 
results shown that the potential functions Huber and Lorenzian do not outperform the 
Quadratic Truncation in any of the 62 cases for PSNR measurements. Furthermore, 
visually Lorenzian creates many noticeable ghost images and shadows, and Huber 
does remove artifacts successfully however it blurs the image more than Quadratic 
Truncation. 

Computation Speed 

For evaluating the three potential functions, the following computation needs to be 
calculatenx (n-l)/2 times for every pixel (where n = filter size x 2 -I). 
Lorenzian: 1 float division, 2 float multiplication, 1 float addition and 1 

float 

log operation 

Huber: 1 comparison, 1 float multiplication or 2 float multiplication 

with 1 

addition 

Quadratic Truncation: 1 comparison and 1 float multiplication 
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In our tests, Quadratic Truncation is the fastest and it performs the best for artifact 

removal with the least degradation to image quality. Therefore, the other two 
potential functions will not be used. 

Threshold (y parameter in the potential function) 

> The y parameter is examined using the quadratic truncation potential function (default 
value is 16), given by: 

e = argmin Y,p(X, -Xj) 

. ,{Xi-Xj)>y 
where, p = { , ' 

^ {X,-Xj)\{x,-Xj)<y 

We can see from the equations that only pixels with similar neighboring pixels are 
affected when varying y. When y is decreased, there will be more y^ as the result of p 
and the intensity of similar pixels output image will be more uniformed resulting a 
more blocky looking type of image (less gradual) for regions with similar intensities. 
Experiments were performed with 2 images (one grayscale and one color image). The 
table in Figure 63 shows the result for the color image. 

We can see from the results that the higher the y, the worst the image quality. 
However, the difference in image quality is not very significant. The pictures look 
very similar and there are no significant changes that can be detected. However as 
predicted, the blocky looking clouds (see Figure 60) are seen when y = 8 or below. 
Therefore by setting y lower we can increase the image quality by some small 
amounts without changing the visual quality of image. Therefore y = 12 was chosen to 
be the default value for y parameter in the potential function. 

Filter Length 

Filter length (F) determines of the size of samples collected for pixel estimations. It is 
obvious that the larger the filter length the longer the processing time. The default 
length given by VM3.0 (B) is 9 pixels. Experiments are performed with varying filter 
lengths to examine how the length affects the image quality (MSE and PSNR) and the 
artifact removal ability. Two pictures are used in this experiment (one color and one 
grayscale). The results are shown in the following tables in Figure 64 and 65. 
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As we can see from the experimental results, the small the length, the better the image 
quality and the shorter the processing time. However, artifacts are not completely 
removed for length = 5 (or smaller) at high compression ratios. The default value for 
the filter length is chose to be 7 pixels to increase image quality and to decrease 
5 processing time without degrading the ability for artifact removal. 

Constraint 

The value of constraint (Thl) affects the filter as described in the equations below. 
y=x + c(d,Thl) 
d=x-x 

c(d,ThJ) = sign(cO * Max{0. absW) - Max(0, 2*(abs{d) ■ Thl))) 

From the equations above, we can see that as Thl increases, more pixels will have 'd' 
as the output for c(d,Thl). Therefore, as Thl increases, more estimate value of x will 
be used as the final estimate and the image quality will decrease while the smoothness 
of an image will increase. In order to keep the image quality as high as possible, Thl 
should be kept as low as possible. However, artifacts might not be properly removed 
if Thl is set too low. Three images were tested with different level of constraints in 
the experiment. The test results are given in the tables in Figures 66 to 68. 

These results show that image quality decreases as constraint increases. However, 
some artifacts (very small, can only be seen when the picture is enlarged) are not 
completely removed when C is equal to 4 or less. In order to achieve the best image 
quality while successfully remove artifacts, the value 8 is suggested to be the default 
value. 

Iteration 

It is obvious that processing time is directly proportional to the number of iterations. 
However, artifacts might not be removed successfully if the image is not processed 
enough times. Figure 69 clearly shows that image quality decreases as the number of 
25 iteration increases. However, when image is enlarged, small artifacts can be still 
visible for Rl. In order to achieve the best image quality while successfully remove 
artifacts, 2 iterations is suggested to be the default setting. 



71 



15 



Copied from 10340491 on 04/01/2005 



wo 00/4957] PCT/CAOO/00134 
Mask Shape 

Different shapes of masks are discussed in Shen's paper, however, only the one with 
the shape of a + sign was implemented in VM3.0 (B). Different masks were 
experimented and the shape of mask is best left unchanged since simplifying the mask 
5 degrades the performance of the filter and complicating the mask increases processing 
time greatly. 

Using the Modified Post Processing Filter 

The adaptive de-ringing is implemented based on the post processing filter in VM3.0 
(B). Few parameters have been eliminated and default values have been established to 
10 increase the usability of the post processing filter. 
To use the modified post processing filter: 

Usage: post2 -s width height bpp -i infile -o outf ile [-1 masic lower threshold) 
I-w mask width] [-t thresh] [-f f_length] I-c constraint] [-r iterations] 

Default values will be used for the parameters if the parameters are not specified in 
15 the command line. The default parameters are: 

mask lower threshold: lower threshold for edge deteciion [30] 
mask width: width of mask [12] 

thresh: y parameter in the potential function (for estimation) [12] 
f_lenglh: filler length (7] 
20 constraints: consirainis using in the clipping function (8] 

iterations: number of iterations [2] 

The modified post processing filter perfomis well with the above default parameters, 
however, these parameters can be changed at the command line if the user wishes. 
The post processing filter can now be applied to pictures with 8:1 compression ratio 
25 with very small increase in MSE. The modified post processing filter seems to 
improve image quality for compression ratio beyond 11:1. 

In practice, there are two ways for the user to play with the settings: on the encoder 
side and on the decoder side. 

Decoder end 

30 The decoder will have the full control of the post processing and the of the associate 
parameters. The encoder will have no control on how the images will look at the 
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decoder end. Accordingly, there will be no addition to file header since no post 
processing parameters will be included. 



Encoder end 

When necessary, the encoder can also predetermine the post processing filter 
5 parameters. Post processing filter parameters will be stored in the image header, and 
the image will be restored according to these parameters. In this setting, the user at 
the encoding end knows exactly how the image will look at the decoder end. 
However, adding these parameters in the header will increase the size of the 
compressed image. 

10 Total of 6 parameters will need to be packed in the header: mask lower threshold, 
mask width, estimation threshold, filter length, constraints and number of iterations. 
The table in Figure 70suggests range limits on the parameters to reduce the header 
size. 



In total, three bytes will be needed to include these parameters in the header. 
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WE CLAIM 

1. A region-based method for encoding and decoding digital still images to 

produce a scalable, content accessible compressed bit stream comprising the 
steps: 

decomposing and ordering the raw image data into a hierarchy of multi- 
resolution sub-images; 

determining regions of interest; 

defining a region mask to identify regions of interest; 

encoding region masks for regions of interest 

determining region masks for subsequent levels of resolution; and 

scanning and progressively sorting the region data on the basis of the 
magnitude of the multi-resolution coefficients. 

;. An apparatus for the region-based encoding and decoding of digital still 
images that produces a scalable, content accessible compressed bit stream 
comprising: 

a means of decomposing and ordering the raw image data into a hierarchy of 
multi-resolution sub-images; 

means of determining regions of interest; 

means of defining a region mask to identify regions of interest; 

means of encoding region masks for regions of interest; 

means of determining region masks for subsequent levels of resolution; and 

a means for scanning and progressively sorting the region data on the basis of 
the magnitude of the multi-resolution coefficients. 

A region-based system for encoding and decoding digital still images that 
produces a scalable, content accessible compressed bit stream and comprises 
the steps: 

decomposing and ordering the raw image data into a hierarchy of multi- 
resolution sub-images; 

determining regions of interest; 

defining a region mask to identify regions of interest; 

encoding region masks for regions of interest 
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determining region masks for subsequent levels of resolution; and 

scanning and progressively sorting the region data on the basis of the 
5 magnitude of the multi-resolution coefficients. 

4. A method for encoding and decoding digital still images to produce a scalable, 
content accessible compressed bit stream comprising the steps: 

10 decomposing and ordering the raw image data into a hierarchy of multi- 

resolution sub-images; 

setting an initial threshold of significance and creating a significance index; 

15 determining an initial list of insignificant blocks; 

forming the list of significant coefficients by encoding a significant map using 
a quadtree representation; 

20 recursively reducing the threshold values and repeating the encoding process 

for each threshold value; and 

transmitting refinement bits of significant coefficients. 

^ 5. An apparatus for encoding and decoding of digital still images that produces a 
scalable, content accessible compressed bit stream comprising: 

a means of decomposing and ordering the raw image data into a hierarchy of 
30 multi-resolution sub-images; 

means for setting an initial threshold of significance and creating a 
significance index; 

35 means for determining an initial list of insignificant blocks; 

means of forming the list of significant coefficients by encoding a significant 
map using a quadtree representation; 

40 a means of recursively reducing the threshold values and repeating the 

encoding process; and 

transmitting refinement bits of significant coefficients. 

45 6. A method of decoding digital still images to produce a scalable, content 
accessible compressed bit stream comprising the steps: 

decoding the bitstream header; 
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deiermining the initial threshold values and the array of initial significant 
pixels, insignificant bits and wavelet coefficients; 

decoding the significance maps; 

modifying the significance lists and decoding the refinement bits for each 
threshold level; 

reconstruct the wavelet coefficient array; 
perform the inverse wavelet transform; and 
reconstruct the image. 

A method of transmission of digital signals that creates a scalable, content 
accessible bitstream comprising the steps: 



pack the most significant bits of the largest coefficients first followed by 
refmement bits and the most significant bits that are significant for coefficients 
20 at the next bit level; 

repeat this process in a recursive fashion until the desired compression size is 
obtained; 

calculate the pack ratios to be used for each channel of the wavelet 
25 decomposition hierarchy by taking the ratio of the two largest amounts of data 

to the smallest amount of data; 

determine the optimal amount of data to allocate for each color channel based 
on the user specified compressed file size; and 

30 if performing region of interest processing, consider the packing overhead 

introduced by the mask when determining the bit budget for each channel. 



10 



15 7. 
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Wavelet 
Family 


Wavelet Type ( Filter Length ) 


Haar 


Haar 
(2) 
















Daubechies 


Db4 
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Db6 
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EJYUV EStandardzed Color 




Test Image 



FIGURE 5 





Coding Parameters 


Availability in RICS 


Rectangle/Square 


(Xmin,ymit) and (width, height) etc. 


present 


Circle/Ellipse 


{xa,yo) and r etc. 


present 


Polygon 


n,hci,y\) , Oc2,yz) ... (^,yn) 


Under development 


Cubic parametric curves 


nM^,y\) . ixz,yz) ...ixn,yn) 


Under development 
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Region M ask Dynamic Object Extracted Dynamic Regions 

FIGURE 8 
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Magnitude vs. Diagonal Position . 
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Common Mask Dimension 


Spectral Filter Size 


Included Spectral Cocfllcients 


32x32 


32 


528 


64 X 64 


38 


741 


128 X 128 


44 


990 


256 X 256 


52 


1378 


512x512 


61 


1891 


1 024 X 1024 


71 


2556 
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Band 


Diagonal Rows 
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Mask Spectrum Size 


Band 1 Rows 


Band 2 Rows 


Band 3 Rows 


Band 4 Rows 
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Image Size (pixels) 


Image Size (bytes) 


Mask Overhead (bytes) 


Mask Overhead (%) 
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546 


13.3 
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FIGURE 18 
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ForLL 
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Bi-level 






Image 





EQW ICU coding with 
bit-plane number = 1 



JPEG 2000 
coded bit stream 
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Gray-scale or 
color Image 



Embedded bit-plane quantization. At 
certain bit-planes, call JBIG routines 



JPEG 2000 code with 
embedded JBIG code segments 
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FIGURE 28 



Y-channel: LU-HL,-LH,-Hll,-HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH, 
U-channel: LL3- HLj-LHj-HHj-HLz-LHrHHr HL,-LH,-HH| 
V-channel: LL3- HL3-LH3-HH3-HL2-LH2-HH2-HLrLH,-HH, 



FIGURE 29 



Ylu-Yhu-Ylh4-Yhh4-Uuj-Vuj- 

YHU-YLHj-YHH3-UHU-ULH3-UHH3-VHU-VLHrVHH3- 

Yhl2-Ylh2-Yhh2-Uhu-Ulh2-Uhh2-Vhl2-'Vlh2-Vhh2- 
Y„LrYLHi-Y„„rUHLrULH,-U„„,-VHLi-VLHrVHHi 



FIGURE 30 



Ylu-Yhu<-Ylh<-Yhh4-Uuj-Vuj- 

YHU-UHLvVHU-YLH3-ULH3-VLHrYHH3-UHm-VHH3- 
Yhu-UhU-VhL2-YuI2-UlH2-VlH2-Yhh2-^HH2-VhH2- 
YHm-UHLI-VlIU-YLHI-ULHl-VLHrYHHl-UHHrVHHl 



FIGURE 31 



Yli>.-Ulu-Vli^-Y„u-Ylh4-Yhh4-U„l4-Ulh.-Uhh4-Vhl4-Vlh4-V„h4- 
Yhu-Ylh3-Yhh4-Uhl4-Ulh4-Uhh4-Vhu-Vlh3-Vhh3- 

YhL2-YlH2-YhH2-UhL2-UlH2-UhH2-Vhu-VlH2-VhH2- 
YHLl-YLHI-YHHI-UHLrULHI-UHHI-VHU-VLHl-VHHI 
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Ya4-ULU-VLu-YHL4- UHU-V„L4-yLH^-UuH-.-VLH-|-YHH4-UHW-V^ 
YHU-UHU-VHL3-YLHrUuu-VLH3-Y„„3-U„„rVHH3 
YHL2-UHL2-VHL2-YLH2-ULH2-VLH2-YHHrUHH2-VHH2 
YHLrUHLI-VHLl-YLHI-ULHI-VLHrYHHrUHHI-VHHl 



FIGURE 33 



Data Type: MUXLIST 
Parameters: 

liTotByiesPacked - long integer total bytes packed into the data buffer for this list. 
cScheme - character processing scheme used for data contained in this list. 

cHighBii - character highest bit-level where data processing begins for this list. 

*pucMuxBuff - pointer to unsigned character buffer where data for each 

bil-levei is packed for this list 

Fields for MUX: information for packing after list processing is complete. 

pliBitPackInfo[]6) - pointer to long integer number of bits packed into the data buffer 

at each bit-level for this list 
liCurBytesCount - long integer current byte count used for bit budget distribution 

when packing this list. 
cCurBitLevel - character current bit-level used for packing this list 

cRemainingBits - character remaining bits to be packed at a given bit-level when 

data to be packed is not evenly divisible by 8 for this list. 
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CALCULATE Channel BitBudget // deiermine optimal bit-budget for each color channel. 
INITIALIZE Channel CunreniBitPlane // highest bit plane that exists in each color channel. 
INITIALIZE liCurBytesCoum and cRemainingBits FOR each MUX list 

FOR Each Color Channel // process each channel separately. 

WHILE Channel BitBudget > 0 AND Channel CurrentBitPlane >=0 

FOR Each Wavelet Transform level // beginning at lowest resolution level. 

FOR Each Orientation Set of Data // according to lossy case natural processing order. 
WcCurBitUvel NOT_EQUAL to Channel CurrentBitPlane 

CONTINUE 
ELSE 

SET BitLevelBytes to p/iBitPactln/olChznnei CurrentBitPlane] » 3 

SET RemBits to /)//i5j/PflCik//i/o [Channel CurrentBitPlane] & 7 

IF SumlcRemainingBits, RemBits) >= 8 
INCREMENT BitLevelBytes by 1 
DECREMENT cRemainingBits by 8 - RemBits 

ELSE 

INCREMENT cRemainingBits by RemBits 
ENDIF 

IF Channel BitBudget >= BitLevelBytes 

INCREMENT liCurBylesCount by BitLevelBytes 
DECREMENT Channel BitBudget by BitLevelBytes 
DECREMENT cCurBitLevel by 1 

ELSE 

INCREMENT liCurBytesCoimt by Channel BitBudget 
SET Channel BitBudget to 0 
ENDIF 
ENDIF 

IF Channel BitBudget EQUALS 0 

BREAK 
ENDIF 
END FOR 
END FOR 

DECREMENT Channel CurBilPlane by 1 
END WHILE 
END FOR 
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Orientation 


Max. Data Size (Bytes) 
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Image Dimension 
(Rows & Columns) 


Overhead Bits 
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Percentage List Overhead vs. Image 
pimrnsinn . 
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Y-channel: R,(LU-HL.-LH4-HH4-HL3-LH3-HH3-HLrLH2-HH2- HLrLH,-HH,). 

R2(LL,-HL4-LH,-HH4-HL3-LH3-HHj-HLrLH2-HH2- HL,-LHrHH,), 
R3(LL4-HL4-LH4-HH.-HL3-LH3-HH3-HL2-LH2-HH2- HL,-LHrHH,), 
R4(LL4-HL4-LH4-HH4-HL3-LH3-HH3-HL2-LH2-HH2- HL, -LH , -HH,) . 

U-channel: RjO-Lj- HL3-LH3-HHJ-HL2-LH2-HH2- HL,-LH,-HH,), 
R2(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HLrLH,-HH,), 
R3(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,), 
R4(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,). 

V-channel: R,(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,), 
R2(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LHrHH,). 
R3(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,). 
R4(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HLrLH,-HH,). 



FIGURE 39 



YHU-YLH3-yn»3-UHU-UuU-UHH3-VHlJ-VLH3-VHH3- 

Yhu-Y'ui2-Yhh2-Uhu-Ulh2-Uhh2-Vhl2-Vlh2-Vhh2- 

yHU-YLH|-YHHI-UHU-ULHrUHHI-V„L,-Vui,-V„„,). 

RzCVLU-YHLrYLH^-YHH^-Uuj-Vuj- 

YHU-yLHrYHH3-UHU-ULHrUHH3-VHlJ-VLH3-VHm- 
YhL2-YlH2-YhH2-UhL2-Ulh2-UhK2-VhL2-'VlH2-VhH2- 
YHI.l-YLH.-YH„,-U„u-UL„,-UH„,-VHLrVLH,-V„„,). 

R3 (YlX4- Yhu- Ylh4- YhH4-Ull3- V|jj- 

Yhu-Ylh3-Yhh3-Uhu-Ulh3-Uhh3-Vhl3-Vlh3-Vhh3- 

Yhu-YlH2-YhH2-UhL2-UlH2-UhH2-VhL2-VuI2-Vhh2- 

Yhli-Ylhi-Yhhi-Uhli-Ulhi-Uhhi-Vhli-Vuii-Vhhi). 
R4 fYtu- Yhu- Yu«- Yhh4-Ull3- V|jj- 

Yhu-Yui3-Yhh3-Uhu-Ulh3-Uhh3-Vhl3-Vui3-Vhh3- 

YhL2-YuI2-YhH2-'UhL2-UlH2-UhH2-VhU2-VlH2-Vhh2- 
YHLrYLHl-YHH|-UHL|-ULHI-UHHI-VHLrVLHrV„H,). 
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Rl (YLU-YH^d-YLH4-YHH4-ULU-VLU- 
YHU-UHU-VHU-YL^U-ULH3-VuH3•VHH3-UHH3-VHH3- 
YHL2-UHL2■VHL2-YLHrULH2-VLH2-YHlI^-UHH2-VHH2- 
YHLrU„U-V„u-YLHrULHrVLHl-YHHrUHHl-VHHl). 

R2(Ylu-Y„w-Ylh4-Yhh4-Uuj-Vlu- 

Yhu-Uhu-Vhu-Ylrj-Ulh3-Vlh3-Yhii3-Uhii3-Vhh3- 

YHU-UHL2-'VHU-YLH2-Uui2-VuH2-YnH2-UirH2-VHH2- 

Y„u-U„u-VHLi-Yun-Uu.i-VLH.-YHi.rUHHrVHHi). 
R3(Ylu-Yhw-Ylh4-Yhh4-Uuj-Vuj- 

Y„U-UhL3-V„l3-YlH3-UlH3-VuI3-YhH3-UhH3-VhH3- 
YHU-UHL2-"VHU-YLH2-ULH2-VLH2-YHH2-'-'llHrVHH2- 
Y„u-UHLI-VHLrYLHl-ULH,-VLHrY„HrUHH|-VHHl). 

R4(Ylu-Y„u-Yuh4-Yhh4-Ull3-Vlu- 

Yhij-Uhl3-Vhu-Ylh3-Ulh3-Vlh3-Yhh3-Uhh3-Vhh3- 
Yhl2-Uhl2-Vhl2-Ylh2-Uuh2-Vlh2-Yhh2-Uhh2-Vhh2- 

Y„L|-UHLrVHU-YLHrULHI-VLHI-YHHl-UHHrV„Hl). 



FIGURE 41 



R1(Yix4-Ull4-Vu.4-Yhl4-Ylh4-Yhh4-U„u-Ulh4-Uhh4-Vhl4-Vlh4-Vhii4- 
Yhu-Yuu-Yhh4-Uhu-Ulh4-Uhh4-Vhu-Vlh3-Vhh3- 
Yhl2-Ylh2-Yhh2-Uhu-Ulh2-Uhh2-V)iu-Vuiz-Vhh2- 

Y„LrYLHt-YHH1-UHLrUu!.-UHHrVHLI-VLHrVHH.), 
R2(Yll4-Uu>4-Vlw-Y,u^-Ylh4-YhH4-Uhw-UlH4-UhH4-V„L4-VuI4-VhH4- 

Y„u-Ylh3-Yhh4-Uhl4-Ulh4-Uhh4-Vhu-Vui3-Vhh3- 

YhL2-YlH2-YhH2-UhU-Ulh2-UhH2-VhL2-VlH2-VhH2- 
YHLrYmrYHHl-U,IL|-ULHrUHH.-VHLI-VLHl-VHHl). 

R3(Yix4-Ull4-Vu^-Yhu-Ylh4-Y„h4-Uhu-Ulh4-Uhh4-Vhl4-Vlh4-Vhh4- 
Y„l3-Yuo-Yhh4-Uhu-Ulh4-Uhh4-V„u-Vlh3-Vhh3- 
Yhl2-Ylh2-Yhh2-Uhu-Ulh2-Uhh2-Vhu-"Vlh2-Vhh2- 

YHL)-YLHrY„Hl-UHU-ULHl-UHH.-VHLrVLH,-V„H|), 

R4(Yix4-Ulw-Vll4-Yhlj-Ylh4-Y„h4-Uhu-Ulh4-Uhh4-Vhm-Vlh4-Vhh4- 
Yhl3-Ylh3-Yhh4-Uhl4-Ulh4-Uhh4-Vhu-Vlh3-Vhh3- 
Yhl2-Ylh2-Yhh2-Uhl2-Ulh2-Uhh2-Vhu-Vlh2-'Vhh2- 

YHLrYuil-Y„Hl-UHLrULH)-UHHrVHL.-VLHrVHHl). 
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YLL4-YHl>rYu)4-YHH4-ULL3-VLU- 




Yhu-Ylh3-Yhh3-Uhu-Ulh3-Uhh3 


VhU-VlH3-VhH3- 


YHL2-YLH2-YHH2-UHLrULHrUHH2 


VHU-VLH2-VHHr 


YHLI-YLHi-YHHrUHLI-ULHl-UHHl 


Vhli-Vlhi-Vhhi (Rl. R2. R3 . R4). 
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YLU-ULU-VLU-YHU-YLH4-Y„„4-UHL-.-ULH4-U„H<,-V,n.4.VLH4-V„„4- 

YhU-YlH3-YhH4-UhLJ-UlH4-Uhh4-Vhl3-VlH3-Vhh3- 

YHL2-YLH2-YHH2-Um^-UL„rU„H2-V„u-VL„2-V„„j- 

YHLI-YLHrYHH.-UHLl-ULHrUHHI-VHLI-VLHrVHH, (Rl, R2, R3. R4) 



FIGURE 44 

MSE vs. Bpp for Glacier Park Image 
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FIGURE 45 

MSE vs. Bpp for Glacier Park 
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MSE vs. Bpp for Glacier Park Image 
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FIGURE 47 

MSE vs. Bpp for Glacier Park Image 
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FIGURE 48 



21/35 



Copied from 10340491 on 04/01/2005 



wo 00/49571 



PCT/CAOO/00134 



MSE vs. Region Data Percentages for Glacier Park Image 
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METHOD AND SYSTEM OF REGION-BASED IMAGE CODING WITH 

DYNAMIC Streaming of Code Blocks 



5 Field of the Invention 

The present invention relates generally to image coding, and more particularly to the 
compression and streaming of scalable and content-based, randomly accessible digital 
still images. 

10 Background of the Invention 

With the rapid growth of the internet and multimedia applications, there is a great 
demand for new image coding tools that that will provide for high quality processing 
capability, an efficient internal architecture, and flexibility in terms of future 
technological advances. This is a challenge that has been put forth before the JPEG 

15 2000 Committee. Although the main focus should be to provide state-of-the-art 
compression performance, JPEG 2000 should also offer unprecedented content based 
accessibility of the compressed format to support applications of various needs. It is 
highly preferred and advantageous that image content be accessed, manipulated, 
exchanged, and stored in compact form. 

20 In order for JPEG 2000 to be the standard coding foundation of new generation image 
processing systems, it must provide for efTiciency in coding, different types of still 
images (bi-level, gray-level, color) with different characteristics (natural images, 
scientific, medical imagery, rendered graphics, text, etc.) within a unified system. In 
addition to providing low bit rate operation with quality performance superiority to 

25 existing standards, this new system should include many modem features as listed in 
the JPEG 2000 requirement document. 

The open architecture and the set of algorithms presented in this document are based 
on Digital Accelerator Corporation's (DAC) Region based Image Coding System 
(RJCS). RICS has not only achieved a rate distortion performance competitive with 
30 best known compression techniques, but also demonstrates a high degree of openness 
and flexibility that will accommodate most well known algorithms as well as new yet 
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to be implemented. The features supported by RICS covers almost all those listed in 
JPEG 2000 Requirements document. 

Generality 

Generality is a primary concern in the architectural design of RICS. The RICS 
5 architecture attempts to be an integrated platform that supports and facilitates a 
variety of applications that may have different characteristics and requirements. For 
example, providing efficient lossless and lossy compression in a single code stream; 
efficient processing of compound documents containing both bi-level (text) and color 
imagery; progressive transmission by pixel accuracy and/or by resolution; random 
10 access of arbitrary shaped regions; and so on. 

Openness 

We also understand that the new JPEG 2000 standard is intended to be a dynamic, 
rather than static, suite of coding tools that support the new generation imagery 
applications and, at the same time, keeping abreast with the progress of technology. 
15 The mathematical foundation for those leading candidates of new image coding 
methods (most noticeably the various types of multi-resolution analysis techniques, 
such as wavelet transforms) is relatively young and still under intensive investigation. 
If an architectural design is based on or restriaed to one or several particular existing 
coding methods, it may become outdated very quickly. 
0 In attempting to produce a flexible and open platform, the RICS architecture 
organizes the image coding process into a set of fijnctionally separable modules, such 
that each module can be developed and optimized individually. Furthermore, all 
modules in the system are designed to be functionally orthogonal with each other, in 
the sense that a new algorithm, without affecting the functionality of other modules, 
5 can effectively replace the base algorithm of any specific module. 

This open architecture will not only be able to accommodate future new algorithms, 
but also makes compatibility with other standards an easy and natural extension of 
many concepts, as will be explained later in this document. 



2 
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Accessibility 

Content based accessibility is becoming an important feature in supporting 
applications such as multimedia database query, internet server-client interaction, 
content production, remote diagnostics, and interactive entertainment. The content- 
5 based accessibility requires that semantically meaningful visual objects be used as 
basis for image data representation, explanation, manipulation, and exchange. With 
images being represented in compressed format, it is desirable to perform retrieval 
operations directly in the code space without requiring image reconstruction. In faa, 
any search algorithms that require image reconstrurtion will be infeasible from a 
10 practical viewpoint because of the huge amount of images in most image databases. 
The previous JPEG and many existing coding techniques focus primarily on the issue 
of compression ratio, paying minimal attention to the need of content based image 
retrieval. 

RICS, by its name, has a fundamental consideration to various types of regions in 
15 images. The regions referred to as ROI (Region of Interest) are usually user-specified 
primitive geometric shapes, such as rectangles or circles. The regions that define the 
visual objects are usually of arbitrary shapes. Regions can also be generated as the 
result of certain mathematical operations or transform properties (e.g. significance of 
transformation coefficients) used to partition the image into disjoint regions of various 
20 (e.g. tiles, hierarchical blocks, or arbitrary shapes). 

In designing the RICS system, DAC considered carefully the distinction between the 
concept of an object and that of a region. RICS is open to very general definition of 
region. A region is a 2D spatial identity with pure syntactic contribution to a code 
stream. In contrast, an object may contain semantic information. Therefore, the 
25 region is perceived as a more elementary and reliable description than the object. 
RICS is region based, not object based. RICS supports a rich set of region types, 
from primitive geometrical shapes, tiles, hierarchical blocks, to most general arbitrary 
shapes. 

The region based coding strategy effectively supports the content-based accessibility 
30 of imagery data. Specifically, this ability enables random access to the code stream as 
well as providing a processing channel for user defined regions of interest or tile 
based techniques. Supporting MPEG-4 object based accessibility is one of the main 
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objectives of the RJCS design. Furthermore, region-based coding provides a natural 
bridge for 'transcodability' with JPEG and JBIG. 

Scalability 

Many applications require image data that is available at different resolutions or 
qualities for decoding. For example, in a progressive transmission process, the bit 
rate control mechanism should allow the image data to be transmitted in certain 
priority order, and be able to truncate the remaining data flexibly, either upon request 
from the receiving terminal or upon channel limitation. Scalable image coding 
involves generating a coded representation (code stream) in a manner that facilitates 
the reconstruction of the image at multiple resolutions or quality levels by scalable 
coding. 

Ideally, the control of scalability should be centralized in a single module as the last 
stage on the encoder side right before the code stream is fed to the communication 
channel. Furthermore, it is desirable that the scalability control module can 
completely handle the required processing locally, without any further request 
propagating back to any previous stages in the encoder. In this way the scalability 
control module avoids the need for multi-pass computation. 

The RICS architecture is designed to support three types of scalability: scalability in 
terms of pixel precision, spatial resolution, and regions. 

Compactness 

There is no doubt that the new image coding standard should offer a higher 
compression performance than the former JPEG, especially at the low bit rate end. 
Integrating the compactness, scalability, and accessibility into a general purpose, 
flexible, and open architecture represents a challenge for JPEG 2000. The RICS is 
designed to provide a solution. The basic idea of region based coding is as follows: 

• The input image data, after certain transformation, becomes a set of image 
primitives. This set of primitives can be wavelet transform coefficients, DCT 
coefficients, other transform coefficients, or even raw image data. 
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• The image primitives are grouped into regions. Region definition can come from 
user defined ROI, from other application modules, or from certain automatic 
segmentation algorithms running in the primitive space. 

• Each region contains one or more independent coding units (ICU). The primitives 
in an ICU are encoded and decoded independently, without reference to primitives 
of any other ICU. This procedure is called the intra-region coding. The outcome 
of an ICU operation is a code block. 

• A muhiplexer (MUX) is employed to integrate the code blocks into the final code 
stream. 

A Brief Description of the Drawings 

Figure 1 is a simplified RICS block diagram. 
Figure 2 is a detailed RICS block diagram. 
Figure 3 is subband decomposition schemes. 
Figure 4 is a Wavelet Filter Library. 

Figure 5 is a performance comparison of YUV and standardized color systems. 

Figure 6 is geometric shapes supported by RICS. 

Figure 7 is a hierarchical partitioning. 

Figure 8 is the relationship between regions and objects. 

Figure 9 is threshold masks obtained from Level 1 Lena Wavelet Decomposition. 
Figure 10 is a Level 1 common mask obtained by thresholding the combined data set. 
Figure 1 1 is individual region masks obtained by separating raw the common mask. 
Figure 12 is individual region masks obtained by separating raw common mask. 
Figure 13 is a spectral magnitude of zigzag spectrum as a function of position. 
Figure 14 is common mask spectrum filter sizes and captured coefficient numbers. 
Figure 15 is a coefficient banding concept used to quantize the common mask 
spectrum. 

Figure 16 is Quantization Band Sizes for Common Mask Spectrum of size 128 by 
128. 

Figure 17 is Spectral Quantization Band Sizes for Other Common Mask Sizes. 
Figure 18 is a Mask Overhead for Grayscale Images. 

5 
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Figure 19 is a low pass Common Mask Approximation for general coefficient 
classification. 

Figure 20 is a Mask Quantization following the Spectral Content for Bit Allocation. 
Figure 2] is Translated Common Masks for Remaining Resolution Levels. 
5 Figure 22 is the three types of ICU structures. 

Figure 23 is a Pavement of a region using type-1 ICUs. 
Figure 24 is a Embedded quad tree flowchart. 
Figure 25 is VLC Codes for EQW. 
Figure 26 is the transcode with JBIG. 
10 Figure 27 is the 'trans-out' and 'trans-in' modes for transcoding with JBIG. 
Figure 28 is Wavelet Maliot Decomposition for Lossy Color Transform Data. 
Figure 29 is Level Priority Processing Order for Each Channel (Lossy Case). 
Figure 30 is Level Priority Processing Order (Lossy). 
Figure 3 1 is Level Priority Color Interleave Processing Order (Lossy). 
15 Figure 32 is Level Priority Color Processing Order (Lossless). 

Figure 33 is Level Priority Color Interleave Processing Order (Lossless). 
Figure 34 is Typical MUX List Data Structure. 
Figure 35 is SNR Progressive Bit Budget Distribution Scheme. 
Figure 36 is MUX List Overhead for Y-Channel 8-Level Decomposition. 
20 Figure 37 is MUX List Overhead for Square Images of Various Sizes. 
Figure 38 is MUX List Overhead for Square Images of Various Sizes. 
Figure 39 is General Region Level Priority Color Processing Order (Lossy). 
Figure 40 is Region Level Priority Color Processing Order (Lossy). 
Figure 4 1 is Color Interleave Region Level Priority Processing Order (Lossy). 
25 Figure 42 is Region Level Priority Color Processing Order (Lossless). 

Figure 43 is Transparent Region Level Priority Color Processing Order (Lossy). 
Figure 44 is Transparent Region Level Priority Color Processing Order (Lossless). 
Figure 45 is a graph representation of the transparent region level in color 
mode for the 4 DCT region channels. 
30 Figure 46 is a graph representation of the region priority level in color 
mode for the 4 DCT channels. 

Figure 47 is a graph representation of the absolute region priority level 
color mode over the 4 DCT region channels. 
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Figure 48 is a graph representation of scaled region priority level color 
mode of the 4 DCT channels. 

Figure 49 is a graph representation of the region percentage priority level 
in color mode over the 4 DCT region channels. 

Figure 50 is a graph representation of the mixed processing multiplexer mode 

of operation over the 4 DCT channels 

Figure 5 1 is a diagram of the basic lossless header structure. 

Figure 52 is a diagram of the basic lossy header structure. ^ 

Figure 53 is a diagram of the basic header tag structure. 

Figure 54 is a block diagram of the basic bit stream syntax for normal modes 

of operation in the lossy case. 

Figure 55 is a block diagram of basic bit stream syntax for region modes of 
operation. 

Figure 56 is a block diagram of the basic bit stream syntax for mixed modes 
of operation. 

Figure 57 is a block diagram of the JPEG transcode/decode method 
Figure 58 is a block diagram of the JBIG transcode/decode method 
Figure 59 is a block diagram of the modified post filtering procedure 
Figure 60 are three gradations of post processing in digital still images.. 
Figure 61 is a representation Edge areas where the de-ringing filtering is 
selectively applies. 

Figure 62 is a table of Test flmage: Camera (256x256 grayscale) Test Image: 
hk (256x256 color) 

Figure 63 is a table of Test Image hk.raw results 

Figure 64 is a table of Test Image camera.raw (grayscale) results 

Figure 65 is a table of Test Image hk.raw (color image)results 

Figure 66 is a table of Camera raw (part a) results 

Figure 67 is a table of Camera.raw (part a) results 

Figure 68 is a table of HK raw results 

Figure 69 is a table of Quality versus iterations. 

Figure 70 is a table of Quality versus iterations.As discussed. 

7 
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A Detailed Description of the Drawings 

The System Architecture 

The architecture of the RICS system can be described as the scheme of dynamic 
streaming of code blocks. In this design, the image primitives of ICUs generate unit 
code blocks. A code block is a logically independent unit of coding and decoding 
which does not rely on information contained in any other blocks. Compression is 
achieved in each block coder, and the coding efficiency of the block coders 
determines the overall efficiency of a RICS system. The openness of the system is 
reflected in the different coding algorithms that can be used to produce different code 
blocks. The system uses a multiplex structure (a MUX) to assemble the code blocks 
into the code stream and realize the bit rate allocation. In short, the RICS architecture 
allow^s for flexibility in the areas of compactness and openness to the block coders 
and, at the same time, allowing the multiplexer to handle the various schemes of 
scalability and random accessibility to the code stream. 

It should be noted that the block used in RICS is a logical unit rather than a geometric 
concept. A code block may correspond to an 8x8 tile (in the case of JPEG mode 
coding), a pyramid data structure in zero tree like coding schemes, a rectangular area 
in block based coding schemes, or an arbitrary shaped area in the raw image buffer. 
The independence of encoding and decoding is the primary requirement of a code 
block. A code block is the outcome of an intra-region coding. 

Figure 1 illustrates the simplified RICS coding architecture. A more detailed diagram 
is shown in Figure 2. A detailed description of the function of each module is given 
in this document. Because the RICS is intended to be an open architecture, DAC does 
not limit each module to any specific algorithm. Instead, any individual algorithm 
can be placed into a module as long as it meets the functional requirement of that 
module. Supporting algorithms are included for certain modules, which reflect 
preferred operation of the overall system. 

8 
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Types of Input Image 

As shown in Figure 2, the RICS architecture supports the coding of three types of 
image data: grayscale, color, and bi-level images. 

Transformations 

5 In a typical multi-resolution coding scheme, an image is transformed via a multi- 
resolution decomposition process. In the proposed architecture, transforms such as 
KL, wavelet, wavelet package, lifting scheme, etc. can be placed in the transformation 
module. These transforms produce a set of decomposition coefficients {Cij} at 
different resolution levels and in different spatial orientations. 

10 The RICS architecture also supports DCT or windowed Fourier transform as a 
transformation technique. This is mainly for the transcodability with JPEG. It should 
be noted that the Fourier based transforms have been studied for more than a century; 
its theory is relatively complete and its mathematical and physical properties are well 
understood. Particularly, its translation, scaling, and rotation properties may be very 

15 useful for content based retrieval computations. On the other hand, wavelet 
transforms are relatively new, and many of its properties require further 
investigations. As a result, the support of DCT may have an impact beyond the sole 
backward compatibility consideration. 

RICS also allows the NULL transformation (that is, no transformation is applied at 
20 all). In this case, an identity transformation is applied to the raw image data as the 
transformation step. The NULL transformation is useful in several instances. For 
example, it is usually not beneficial to apply DCT or wavelet transforms to bi-level 
images (text) for compression purpose. Another example is the residual images in 
video coding. The information in a residual image is the difference between video 
25 blocks and has a high frequency spectral content already. It is highly questionable 
whether it is beneficial at all to apply another mathematical decomposition (DCT or 
wavelet) to this type of data for the purpose of compact coding. 
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Region Definition 

The function of this module is to partition the coefficients produced by the 
transformation module into a number of spatial regions. The RICS supports three 
region schemes. 

1 . Automatic partition based on the preordering of the coefficients. 

2. Partition based on user defined ROI or object related shapes. 

3. Partition based on tiling. 

The first partition scheme is also referred to as the region hierarchy formation process. 
Consider the transformation coefTicients as a set. The region hierarchy formation 
process partitions the set into a number of hierarchically disjoint subsets according to 
certain definitions of significance. In providing a very general partitioning technique 
captured in a very general region based control architecture, RICS can perform highly 
flexible progressive transmission modes of operation that depend on data priorities set 
up for the code stream. 

Scheme 2 deals with user specified ROI, typically primitive geometrical shapes, such 
as rectangles or circles, as well as object related shapes. Object related shapes could 
come from a variety of sources, such as user input, motion analysis in video 
compression, etc. 

Scheme 3 is a simple partition and requires minimal for shape coding. This scheme is 
essential in the JPEG mode. Some wavelet based compression techniques utilize this 
scheme to explore coding efficiency. This scheme offers very little support for 
content based accessibility. 

Hierarchically disjoint regions can be used in combination with user defined ROI in 
still image processing and objects in video processing. However providing a fair user 
partition for detail information is difficult in still image compression. But automatic 
partitioning or preordering techniques can be performed to control user selections in 
both arbitrary and block based modes of operation. The research presented in this 
document introduces a new multi-level control architecture for advanced multi-level 
image processing. The ability to perform a host of partitioning and preordering 
techniques in both normal and region based modes of operation is encompassed in 

10 
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this architecture. Both arbitrary and automatic region formation schemes are handled 
in the same manner at a high level. 

Region Shape Coding 

Three types of region shape are supported in RJCS: tiling, primitive geometric shapes, 
and arbitrary shapes. 

1. Tiling requires only a small set of parameters to describe the configuration 
especially in sequence based processing modes. However the sequences can be 
organized and presented in many ways in packing them into the fmal code stream. 

2. Primitive geometrical shapes can also be coded eflficiently. For example, a 
rectangle can be defined by four integers (jTmin, >'min)and (xmix, ymu^, or 
(jCmin, ^min) and (width, height) etc. 

3. Arbitrary shapes are more difficult and costly to encode. Partitions produced 
automatically by image analysis algorithms may contain many small regions. 
Specifically, the auto-region detection routines presented here produces 
hierarchically organized data partitions. This presents a highly challenging 
problem for shape coding. RICS provides two practical solutions to this problem. 
One is a quad tree based hierarchical region description. This mode of operation 
follows bit level ordering at different resolution levels (see embedded quad tree 
wavelet (EQW) in Ch. 4.) The other is a DCT coded multi-level region channel 
definition followed by preordering the partition (given the restraints of the code 
stream in terms of overhead). Though the coding efficiency of the second solution 
is currently slightly poorer than the first, it does contain a number of attractive 
features: 

• It provides a single representation for multi-level bitmaps. 

• From this single representation, region masks can be reproduced at arbitrary 
resolution levels, which is useful for subband coding. 

• It generates curved shapes, which potentially (e.g. in case of large number of 
complex curved regions) could outperform quad tree based region 
descriptions. 

11 
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• It has several useful properties, such as translation and rotation properties, 
which the quad tree based description does not offer. In fact, a quad tree 
representation will change dramatically if an object is slightly translated in the 
image, making this type of representation not suitable for content based 
retrieval applications. However, the main advantage of the quad tree 
representation is its simplicity and most noticeably in block based modes. 
The region shape code is included in the code stream when the region definition is 
determined at the encoder. In the case of decoder specified regions, the region shape 
coding has a whole new meaning. 

Intra-Region Coding 

The flinrtion of the intra-region coding is to arrange the transformation data in an 
arbitrary shaped region into a one dimensional code stream. Regardless of the region 
definition scheme or the region coding technique, this streaming process requires an 
intermediate state where a control architecture can be designed to tailor the region 
channels whether dealing with a bitmap mask, an auto-detection routine or any other 
of numerous classification techniques. At the decoding end, the inverse routine 
generates the same mask to unpack the values fi^om the one dimensional code stream 
and place them back into the correct position in each region. 

Intra-region coding is completed in block coders. Different block coders can be used 
to produce the ID code stream. For JPEG mode, the zigzag scanning/quantization 
routine is called to pack an 8x8 DCT block. For wavelet based coding, both explicit 
quantization and implicit quantization schemes can be used. In particular, embedded 
zero tree and embedded quad tree can be used as implicit quantization schemes. 
Furthermore most implicit quantization schemes are implicitly decodable. In dealing 
with bi-level images, a JBIG routine can be called a block coder. This effect can be 
staged by calling an implicit quantization scheme using one bit plane. Alternatively, 
efficient JBIG routines can be called block coders in an embedded coding scheme at 
each of the multiple bit planes. 

The Multiplexer 

The function of the multiplexer is to assemble the code blocks derived from different 
subbands and different regions in proper orders into a single code stream. Due to the 
12 
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richness of region definition and block coding schemes, there are plenty of ways to 
pack the code blocks. Different ways to merging the data lead to different 
transmission priorities. For all transmission ordering modes, the final code stream can 
be transmitted progressively, and can be truncated at any desired place. 
5 The notion of using a multiplexer has been adopted in some standardized processes 
such as MPEG. In contrast, the use of a multiplexer in the design of a still image 
coding system is rare. The region based coding strategy opens the opportunity to 
systematically explore the syntactic and semantic richness in code stream ordering 
and transmission medium. With the image segmentation techniques improving in the 
10 future, region shapes will become more accurate in tracking objects in the image. In 
that case, the multiplexer will not only work as a syntactic composer, but also impose 
semantic meanings to the code stream. 

Hi ghlight of Features 

The RICS is a true open architecture. It supports not only the algorithms DAC has 
15 developed, but also can accommodate most existing well known compression 
algorithms. Users may include their own functions that are appropriate to their 
application in a number of modules such as transformation, region definition, region 
shape coding, and intra-region coding all under the MUX control architecture. It is 
also implicitly flexible to new technological advances. 
20 • DAC's current implementation offers superior low bit rate performance that is 
competitive with best existing compression techniques. 

• The richness in region definition in the RICS allows great accessibility, thus 
providing a solid foundation for content based image applications. 

• It covers compression for both continuous tone and bi-level images in a single 
25 unified architecture. 

• It provides lossy and lossless compression in a single, natural, code stream in the 
course of progressive decoding. 

• It provides a variety of progressive transmission modes that allow images to be 
reconstructed with increasing pixel accuracy or spatial resolution using region 

30 priority modes for user specified ROI or system defined significant areas. 

13 
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• It supports both fixed rate and fixed size modes. 

• It supports very flexible random access to and processing for regions with 
arbitrary shape. 

• It provides a graceful backward compatibility (or transcodability) with JPEG. 

• The generic region definition in RICS is a very suitable interface for object-based 
coding schemes currently under development by MPEG-4. (In fact, as the 
arbitrary region shape begins to fit dynamic objects with more accuracy, motion 
estimation is more definite, and consequently allows for more efficient error 
compensation using region-based coding.) 

• Our general region shape definition provides a solid basis for object based 
composition and object based information embedding. Multiple objeas with 
arbitrary shape are accepted. 

• Our bit stream is robust to bit errors. Each unit structure used by the MUX is 
independently decodable. 

• Incorporating standardized encryption techniques with the RICS system is 
straightforward. 

Transforms 

In order to support multiple application needs, provide transcodability, and 
accommodate future growth, the RICS architecture is designed to support three 
categories of transformation: the DCT, wavelet transforms, and a special NULL 
transform. 

DCT 

The RICS architecture supports the DCT transforms as defined in the baseline JPEG. 
Wavelet Tran.sfnrm 

The RICS architecture supports various kinds of subband decomposition schemes, 
including the three schemes in Figure 3. 



14 
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The Null Transform 

In dealing with the compound documents with mixed contents, it is sometimes 
required to encode some areas of the image without any transform. In particular, 
regions with bi-level pixels usually do not need to be transformed for coding purpose. 
5 In coding these regions, the transform stage is bypassed; the effect is simulated by a 
NULL transform stage in the process. 

The present version of RICS uses totally 34 types of wavelet kernels as given in the 
table in Figure 4. In our experiments, the biorthogonal filters (Bior9-15 of the table) 

10 seemed to give the best compression. Both low pass filtering and high pass filtering 
are done using convolution. Afler filtering, the wavelet coefficients are down 
sampled by 2. This process is repeated using the low pass part until the desired 
decomposition level is reached. In inverse wavelet transform, quadtrature mirror 
filters (QMF) for low and high pass are used. Then the coefficients are up-sampled 

15 by 2. 

Color Transform 

In addition to YUV format, the RICS system uses the Karhunen-Loeve Transform for 
color transformation. Following the notation of statistics we term this process as the 
color standardization. The KL color transform requires more computation than YUV 
20 transform. Figure 5 shows the PSNR comparison of the two color transforms. More 
test results are also available from DAC. 

Region Definition and Shape Coding 

The notion of region processing plays a fundamental role in the operation of RICS. 
25 The choice of region based coding is motivated by many application needs such as 
content based retrieval, interactive multimedia application, graphics object and image 
composition, and coding of dynamic object in video compression. 

There are two fijndamental issues in a region based coding strategy: defining regions 
and coding the regions. In the RICS system, region coding is divided into the shape 
30 coding and the intra-region coding (the content coding). Section 3, describes the 

15 
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various schemes for region definition and shape coding. The intra-region coding is 
discussed below. 

From the viewpoint of JPEG 2000, the process of region definition does not have to 
be standardized. However, schemes for forming regions that the JPEG 2000 will need 
to suppon should be defined clearly. The RJCS supports three types of region 
definition. 

• Tile based organization 

• primitive geometric representations 

• arbitrary shapes 

Region definition is an optional step: an image can be coded wathout specifying any 
region (the non-region mode). In this case, the entire image area is considered as a 
single region. 

Tiling: Definitions and Shape Coding 

Tiling is perhaps the simplest region definition scheme. In this scheme, the entire 
image area is divided into a number of rectangular blocks. In particular, 8 by 8 tiles 
are used in the JPEG coding mode. Most tiling approaches cannot cover the natural 
shapes, region trends or partial objects in a picture. 

For coding a tiling scheme, a small set of global parameters will be sufficient. The 
primary parameter is the size of tiling block and the technology is developed around 
this theme. In case of subband decomposition, tiling block size may vary from one 
subband to another. However, the RICS architecture has a sound operational 
foundation. It can be operated in both tile based or arbitrary processing modes of 
optional arbitrary region formation schemes. 

Primitive Geometrical Shapes 

Primitive geometrical shapes are ideal for supporting user interaction; i.e. user 
specified ROI, and for processing compound documents where the text areas can be 
well covered with one or more rectangular boxes. Geometric shapes currently 
supported by RICS are listed in the Table in Figure 6. 

16 
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Arbitrary Regions and a Generic Binary Partition Hierarchy 

Let Co be the set of image primitives. Given an ordering relationship >- , generic 
binaries partition of Cointo Ci^and C,, is defined by the following conditions. 

(1) C,puC„ =Q. 

(2) C,ooC„ =<D,and (Eq.m.l) 

(3) for any a e C,o and A e C„ , a>b holds. 

Recursively, each of the new generated subsets can be further partitioned into two 
subsets, resulting in a hierarchical structure to represent the partition (see Figure 7) 

A useful ordering relation popular in the compression community is "more significant 
than": a>-h whereby a is more significant than b. The definition of significance can 
be very general. One definition may be used to create an ordering, or several 
definitions are combined together to generate an ordering. This can be illustrated by 
the following example shown in Figure 8. Suppose the region masks shown in the 
leftmost image are gerierated by considering the magnitude of wavelet transform 
coefficients as the measure of significance. The region mask in the middle image is 
an object shape generated by another definition of importance that, for example, can 
be the intensive dynamic region in a video frame. The rightmost image shows the 
intersection of the two region masks, which is a new region partitioned according to a 
multiple definition of significance. 

As an example, if the significance is measured by magnitude of wavelet transform 
coefficients, then the ordering relation is simply the numeric ">" relation. To create a 
binary partition, a set of thresholds is sufficient. The first threshold, say T^, separate 
the initial set Q into two subsets, C,o and C,,. Subsequent threshold values can be 
specified to recursively partition the new subsets. 

The above partition scheme produces accurately a number of disjoint subsets which 
satisfies the strict ordering relationship of >■ . However, it usually produces many 
small, scattered regions. Consequently, the representation of region shapes for these 
17 
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subsets can be potentially expensive. From the representation point of view, fewer 
numbers of smoothly shaped regions are desirable. In the following, a less strict 
partition scheme is defmed, which takes the above partition scheme as a special case. 

Instead of requiring the condition ayb to be held for all aeC^^, ieC„, it is 
replaced by a looser condition 

->beC^, b>~ a, for any a eCf,. (Eq.in.2) 

Intuitively speaking, this definition defines a partition in which no element in C, 
proceeds any element of . In other words, some elements belonging to C, (by 
Eq.m.l) may now be placed in Co, but, no element belonging to Cq C, (by Eq.m.l) 
can be placed in C, (Eq.in.2). In doing so, the subset Q vnll be enlarged, and the 
region shape of will be smoother. At the same time, the enlargement of Q is 
controlled to such an extent that the loss of classification accuracy is limited to a 
reasonably low level. 

Another scheme, which is symmetric to the Eq.ni.2, can be defined similarly for the 
following condition. 

-.a e Co , a K 6 , for any i> e C, . (Eq.ni.3) 

In this case, we ensure that no element in Co is proceeded by any element in C, . 

It should be pointed out that there is a tradeoff between the enlargement of C^ and the 
space requirement for the region shape representation. It is possible (and indeed quite 
often) that there are some isolated scattered points that belong toCj . But, in order to 
exclude those isolated points fi-om C, , we have to sacrifice many elements fi-om C, . 
This will result in a very large Co , a degraded classification. An approximation to the 
scheme of Eq.III.2 is introduced to rertify this problem. A small number of elements 
in Co are allowed to be "misclassified" into C, . Those misclassified elements can be 
separated later in the intra-region coding stage that generally follows natural 
processing orders available for further classification of the partitioned elements. 

18 
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When such region hierarchy formation schemes are applied to image primitive planes, 
they serve as preordering processes that set up a partial ordering among the 
primitives. Then, the problem of coding this ordering becomes a geometry problem 
for coding the shapes of the partitioned regions. Transformations such as wavelet are 
5 known to use a joint spatial/frequency domain methodology. The region hierarchy 
formation scheme mentioned above has an intuitive geometric explanation for this 
class. Transformations such as DCT and Fourier transform are purely frequency 
domain methods. It is interesting to note that when the above partitioning scheme is 
applied to a 2D DCT transformation plane, the generated regions are roughly 
10 circularly shaped bands. The popular zigzag sequence used in JPEG turns out to be a 
simple and reasonable approximation of these bands. Therefore, the region based 
coding strategy accommodates naturally both spatial/frequency domains and pure 
frequency domain transformations. 

Automatic Region Formation in the Wavelet Transform Domain 

15 DAC has developed an automatic region formation scheme that is used to categorize 
wavelet coefficients in the transform domain according to magnitude. The procedure 
will be outlined in the following sections. Also included are many of the design 
details used to develop the first version of the general region classification scheme. 
The advantages and disadvantages of the technique as well as performance and 

20 overhead issues will be discussed. Finally conclusions will be drawn to summarize 
the results. 

Detecting the Regions 

The wavelet transform is used to decompose the original image information into a 
multi-resolution hierarchy of data that is more suitable for compression. The result of 

25 completing one pass of the 2D wavelet transform is one low pass part (LLi) and three 
high pass parts (HLj, LHi and HHi.) The transformation procedure is repeated using 
the low pass LL part as the starting point for each succeeding pass. The LL part is an 
approximation of the original image at a lower resolution. At each resolution level 
the HL parts carry the vertical, the LH parts carry the horizontal and the HH parts 

30 carry the diagonal edge information. The level of detail information contained in any 
particular orientation is specific to the local response generated in the choice of 
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wavelet kernel. The original image is transformed into a unique band pass response 
hierarchy of detail information. In any compression process, it is the highest energy 
coefficients that are considered first for making a contribution to the reconstructed 
image. It is this transformation property that is used to partition the coefficients into 
regions of interest corresponding to high energy locations in the original image space. 
The magnitude of the coefficients in a particular subband is related to the response 
that the original image information has to the transform kernel. Shaqj edges in the 
original image (or a lower resolution LL part) are indicated by an increased kernel 
response in terms of coefficrent magnitude, in the general vicinity of a specified area 
of concern. Other image changes that are not as abrupt will translate into a gradual 
response that is less significant. In those areas where little or no change occurs, the 
wavelet transform will indicate little or no response. 

Level 1 Kernel Response Threshold Images 

DAC's automatic region formation scheme makes use of coefficient magnitude 
information to categorize the data. The starting point for the procedure is at the 
highest resolution level of the wavelet transform hierarchy looking at the three detail 
data sets (HLi, LHi, HHj.) A 2 bit representation of each detail orientation for a 256 
by 256 Lena image is given in Figure 9. These images are obtained by using 
threshold values to categorize the coefficients by magnitude. The wavelet kernel used 
in this case is the standard 9-7 filter set implemented in a lifting scheme. 

A decaying histogram procedure is used to determine the threshold values in each 
case. In this particular analysis the coefficients are partitioned into regions according 
to decreasing levels of magnitude; 10% of coefficients are partitioned into region 1, 
15% into region 2, 25% into region 3 and 50% into region 4. Note that the largest 
coefficients appear darker in the images while the smallest appear lighter. The 
threshold values are calculated automatically for each orientation. The values used to 
generate the partition in each case are given in Eq.m.4 (assume Q is the magnitude of 
the coefficient under consideration.) 



20 
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HLi: 0 < Ci < 2, 2 < Ci < 4, 4 < Ci < 8 and Ci > 8 

LH,: 0 < Ci < 3, 3 < Ci < 5, 5 < Ci < 14 and Ci ^ 14 (Eq.m.4) 

HH,: 0 < Ci < 2, 2 < Ci < 3, 3 < Q < 5 and Ci ^ 5 

The set of threshold images appearing in Figure 9 is motivational in determining a 
10 region formation scheme based on this approach. It is apparent that wavelet 
coefficients can be classified in this manner to form a mask. However, the large 
overhead required to code each individual mask must be solved in order to make this 
type of classification scheme viable for implementation. 

Formation of the Raw Common Mask 

15 The first step taken in the development of a technique designed to reduce the mask 
overhead is to consider a common mask approach that can be used to capture the most 
important coefficients between all three orientations at the lowest resolution level. 
Given that the data range is different in each orientation, the data is normalized to the 
largest range in an absolute value sense. This step is taken in order to compensate for 

20 the range differences that exist between the data sets, and to ensure that corresponding 
pixels from each orientation can make an equal contribution in the formation of the 
common mask. This step is illustrated in Eq.ni.5. 

HL", = NewRange(HL,) = ScaleRange(MaxRange(LH,. HL, , HH,)) 

LH", = NewRange(LH,) = ScaleRange(MaxRange(LH,, HL, . HH,)) (Eq.m.5) 

HH", = NewRange{HH|) = Sca]eRange(MaxRange(LH|, HL, , HH,)) 

30 Once the data ranges for the smallest two sets have been scaled to that of the largest, 
the common mask values can be determined. The next step is to generate a new 
normalized data set by taking the maximum value between scaled orientations at each 
pixel location. Let C (i) be the new coefficient obtained in this procedure and let 
C HLI , C LHi and C'hhi be the individual scaled values from each orientation. The 

35 new data set is formed using the step illustrated in Eq.in.6. 

C (i) = MAX(ABS(C'HLi(i)), ABS(C'HLi(i)), ABS(C'HLi(i))) (Eq.m.6) 
21 
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The final step in forming the common mask is to threshold the new coefficients based 
on the histogram decay and the percentage of the total data amount to encapsulate in 
each region category. The same values mentioned earlier are used in this case (10%, 
1 5%, 25% and 50% ) The raw common mask categorization image appears in Figure 
5 1 0. The threshold values used to generate the mask are given in Eq.III.7. 

H ,: 0 < Ci < 5, 5 < Ci < 1 1, 11 < Ci < 26 and Ci > 26 (Eq.m.7) 

10 

Note the level of detail contained in this image. High frequency edge information 
representing the most important coefficients appears in the darkest areas. Low 
frequency changes representing areas of the image that are relatively uniform appear 
in the lightest areas. The combined kernel response from each orientation follows the 
15 high energy areas that exist in the original image. 

Common Mas k Conditioning Using the DCT Method 

Generally, the raw size overhead of the common mask obtained in the previous 
section is still too large to be included in the final bit stream for most image 
20 processing needs (i.e. an overhead of 2 Bpp for this particular 4 region mask). This 
seaion outlines the next step in the automatic region formation scheme. The DCT 
transform is used to reduce the overhead size to an acceptable amount. 
The DCT is not applied directly to the common mask data in its mult-valued form. 
Rather it starts by considering each region contained in the common mask as an 
25 binary data set such that the sum of the individual data sets is equal to the multi- 
valued mask in functional form. The DCT transform applies the transform of the sum 
of functions in the same manner as taking the transform of each individual function 
and summing the results. Furthermore there are unique spectra in each case that 
should be considered separately for a fully imegrated analysis. The individual region 
30 masks obtained by decomposing the raw common mask for the Lena image appear in 
Figure 11. Notice that although there are only three images, the fourth region channel 
(the background) is implicit. 

The 2D fast DCT algorithm is used to transform each binary data set into the 
frequency domain for spectral analysis. The spectra for the upper three regions 
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appear in Figure 12. Please note that a log scaling of the coefficients is used in the 
images to exaggerate the appearance for display purposes. Notice how the low energy 
content is dispersed through the spectrum in each case. Even at such low energy 
levels, the energy is localized in some areas. 

It is apparent from the spectral results of Figure 12 that the magnitude of the DCT 
coefficients decreases rapidly in moving down each spectrum. The interaction of the 
content retained for each partition is used to guide the classification procedure. The 
interrelationships that exist in each binary spectrum must be investigated further. 

A plot of the overall summed spectra in short integer format appears in Figure 13. 
The conventional zigzag ordering is used to re-group the coefficients before 
quantization. Note that only the first half of the data in reflected in the plot. The 
outer 8K are omitted since they do not change much from the DC level. 

The first consideration observed in the plot is that the coefficient magnitude 
deteriorates at an alarming rate over a relatively small number of coefficients in the 
first portion of the spectrum. The second observation is that most of the spectral 
energy is concentrated in the first 12% of the data. Together these observations 
suggest that implementing a low pass spectral filtering stage, followed by careful 
quantization will reduce the amount of data that must be retained to capture the region 
formation concepts at the lowest level ( in this common mask partitioning scheme). 
Note the apparently large dc offset. Part of this level is introduced by scaling the 
original wavelet data sets to a unified range. The rest of the offset is due to the DCT 
transform. 

T^w Pass Filtering the Common Mask Spectrum 

The total mask spectrum must be filtered carefully to reduce the excess content that 
has little contribution to the region trends that exist in the retained coefficients. The 
main concern is to determine how to quantify the filter size such that a dynamic 
implementation may be obtained that will work accommodate any arbitrary image 
sizes. 

The first step taken towards finding a dynamic solution to this problem comes from 
experimentation. In looking at the result of using zigzag ordering, coefficient 

23 
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magnitude and significance deteriorate at some exponential rate. Experimental 
evidence suggests that as the size of the original image increases, the number of 
coefficients that must be retained to construct a similar quality mask increase on a 
logarithmic scale. This is an important observation since doubling the original image 
5 dimensions translates into a log 2 increase, in the number of retained spectral 
coefficients. 

Assume that 3 is a general log base to begin our discussion. If the filter size for an N 
by N common mask spectrum has been determined experimentally as ^u, then what is 
the size of the filter for a 2N by 2N mask spectrum? If the unknowm filter size as <j)2N, 
10 then the following simple logarithmic relationship can be used to compute the new 
filter size. 

(Eq.m.8) 



A number of experiments were conducted to determine a reasonable filter size that 
can be used as a base for filtering common mask spectra. It has been determined that 
for masks of size 128 by 128, a filter size of 40 to 45 diagonal rows yields reasonable 
25 results. The table in Figure 14 give the filter sizes for raw mask sizes using 1.375 for 
p. The logarithmic base P can be used for final calibration in the current design 
implementation. 

Notice the filter size specified for a common mask spectrum of 32 by 32 is given as 
32 diagonal rows. This value is selected to calibrate the filter sizes for the larger data 
30 sets. The results indicate the filter size required for a first cut of the common mask 
coefficients does not have to be exact since the distribution of the spectral content 
may change slightly fi-om one image to the next. However, it must be large enough to 
capture the majority of the most important spectral content for each region partition. 
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Spectral Quantization Techniques for Filtered Coefficients 

DAC has developed two techniques for quantizing the filtered mask spectrum. The 
first approach is basically a modified uniform quantization procedure that is used to 
verify the concepts. It proved useful in confirming the algorithms and developing the 
5 concepts, but it lacks robustness and flexibility. The second technique is a general 
approach that is much more robust and elegant. It is dynamic since it tracks the 
magnitude of the spectral content in determining how to quantize the coefficients. 

A First Quantization Approach 

The first approach combines the filter dimension determination into the quantization 
10 procedure. A base mask size of 128 by 128 is use as the starting point in this 
approach. Assuming that a logarithmic relationship exists for coefficient importance 
along the diagonal order, a quantization technique can be developed to exploit the 
relationship. 

The main premise is that there are bands of spectral coefficients that can be used to 
15 guide the quantization procedure. The importance of the coefficients in each band 
decreases along the spectrum. This concept is illustrated in Figure 15. 

The number of diagonal zigzag rows to include in each band is a function of the filter 
size and the original common mask size as well as the spectral content. As a fu-st 
approximation, a table look up technique is used to determine the number of rows to 
20 use in each band. The base band sizes determined experimentally to give reasonable 
partition and reconstruction results for a 128 by 128 mask are given in the table in 



Notice that the number of diagonal rows doubles in each step. The number of bits 
25 selected to quantize each band decreases from one band to the next in the following 
fashion. 



Figure 16. 



Band 1 



16 bits 



Band 2: 



10 bits 



Band 3: 



8 bits 



30 



Band 4: 



8 bits 



25 
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A series of band sizes were found experimentally of mask dimensions by powers of 2 
increments. The results are in the table of Figure 17. These values are used to 
partition the spectral content for quantization. This information together with the 
number of bits used in each band and the filter size measurement determines the 
5 common mask overhead of an image. 

The mask overhead based for this quantization scheme is tabulated in the table in 
Figure 18 for the corresponding image sizes. There is an additional overhead header 
included in this result for mask reconstruction on the decoder side. Note that the 
mask overhead is the same size for both grayscale and color images since it is 
10 generated from the intensity information. However, the percentage overhead in the 
color case is reduced by a factor of 3 since it can be distributed over 3 channels. 

After applying the inverse DCT on the quantized coefficients a low pass 
approximation of the original common mask is obtained. It is this result that is used 
to guide the region channels in classifying and processing the wavelet coefficients. 
15 The result indicated in Figure 19 is obtained by applying this technique to the first 
wavelet transformed level for a 256 by 256 Lena image. From Table 17, the overhead 
for this mask is 1021 bytes. 

Notice the high energy changes have been eliminated by low pass filtering (especially 
the changes that occur over small portions of the image). However, the mask does 
20 succeed in capturing most of the essential parts of the original image. This mask is 
common to all orientation at this level and used to partition the coefficients into 
region categories. 

A Second Quantization Approach 

There are some apparent problems with the underlying procedures used for 
25 quantization in the first approach. The first of these is that the coefficient band 
classification is fairly rigid and not that flexible. The filtering and quantization stages 
are tied together and fully independent. Affecting one parameter has consequences to 
the other. It is true that a logarithmic relationship exists in the spectrum that can be 
exploited to determine a reasonable filter size, but the quantization step must be 
30 tailored fi-om the filtering stage. 

26 
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DAC is currently implementing another quantization scheme for mask filtering that is 
much more dynamic and flexible. Quantization is no longer tied to the filtering stage 
and the coefficient banding technique of the previous approach. The current 
technique follows the actual content of the spectrum in forming quantization 
5 categories. In this way, coefficients are grouped into categories according to the 
energy gradient of the spectral distribution. The concept is illustrated in Figure 20. 

The rate of spectral decent together with position information within the spectrum can 
be used to quantize the number of bits used to categorize the data. This process can 
be controlled precisely. The procedure can be calibrated as required by studying the 
10 effects of changing image dimension in terms of an optimal quantization change. 

The previous quantization technique can be used to estimate the mask overhead for 
the new scheme proposed here. Consider the same base conditions presented for the 
15 banding approach applied to the 256 by 256 Lena image. The spectral filter size in 
this case is 44 diagonal rows. If 4 quantization intervals are used at 16, 12, 8 and 4 
bits each the mask overhead can be estimated. Let Ni be the number of coefficients in 
each interval, bi the number of bits used to quantize this interval, and Om be the total 
mask overhead. A rough estimate of the new packed size is calculated in Eq.III.9. 

■3 

25 (Eq.III.9) 
Om = 6x 16 + 39x 12+ 186 x 8 + 804x4 
Om = 5268 bits (-659 bytes) 

30 

Even though the actual overhead has not yet been determined, fi-om the insight gained 
in using the first quantization approach a saving of 30% is not unreasonable. This 
promises to be a substantial saving in the mask overhead. There are ways to utilize 
35 saved overhead. The first is to reduce the amount of mask information in the code 
stream. This approach will retain the same mask quality leaving extra space for 
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refinement information in the code stream. The second way to utilize the saved 
overhead is to retain a better quality mask in keeping the overhead size constant. 

The Common Mask and Multi-resolution Classification 

A translation technique is required to apply the region mask at other resolution levels. 
The simplest approach is to look at the significance of the mask coefficients in a small 
block to determine an appropriate value for the next resolution level. This process is 
repeated until a mask for each resolution is obtained. The most logical approach is to 
take the largest coefficient through to the next level. 

Mv(L*ifj = Max(Mv(L)2,^ , Mv(Lpj+i^ , Mvoft^i , Mv^js+i^i) (Eq.in. 1 0) 

Some heuristic approaches have been investigated for resolution translation. One 
approach is to make decisions based on the sum of the individual mask coefFicients. 
The upper level masks are affected by either enlarging or shrinking the encompassed 
coverage regions. However, there is currently no conclusive evidence confirming the 
validity of this approach and more experimentation is required. 

It is worth mentioning at this point, that DAC has begun investigating the wavelet 
kernel ROI mask formation technique of verification model (VM) 3.0 as a method of 
resolution translation for other transformation levels. Preliminary investigations 
indicate this technique introduces a gradient bias for regions of higher classification 
level in the translation. A combined hybrid approach may have some advantages. 

The rest of the common mask hierarchy obtained by using the simple 4 to 1 resolution 
translation techniques and observing the most important region as key is given in 
Figure 21. The complete set of region masks is used to classify regions of importance 
in the wavelet transform domain according to coefficient magnitude. 

Misclassification and the Common Region Mask 

It should be noted that while it is true that misclassification occurs in this technique, 
there is a tradeoff between code stream overhead and region channel accuracy. The 
original raw common mask size is 32 kBits for a 4 region mask. The basic uniform 
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quantization scheme compresses this amount to 8 icBits. The modified scheme 
suggests that mask compression of 5 or 6 to 1 can be achieved with a tradeoff between 
mask accuracy and code stream overhead. 

The low pass filtering and spectral quantization steps can be tailored to deHver similar 
5 quality masks for most images in a dynamic implementation. In addition the partition 
can be studied in regards to how the original grouping is affected by the fihering step. 
For the Lena image example the original partition is 10/15/25/50%. After filtering the 
coverage ratios used for bit budget calculation and the code stream distribution (see 
multiplexer in Ch.V.) changed to approximately 9/31/33/27%. There may be some 
10 advantages in adjusting the original partition such that the data coverage ratios are 
controllable in the final stages. 

The question of miscoverage must be addressed in all lossy region formation 
techniques. It is the subsequent modules that are affected by the misclassified 
information. Thus the interaction between the DCT region formation module and the 
15 subsequent sorting and packing modules must be understood in any attempt to obtain 
the best region channel results for this classification scheme. 

Handling Misclassification in a Region Formation Scheme 

Subsequent stages in the compression algorithm are already tailored to process 
20 information in a descending level of importance. Most wavelet implicit quantization 
procedures process the coefficients in bit plane order of importance. What this means 
in terms of any region formation scheme is that the information encompassed in each 
region is further partitioned by bit level processing. In the common mask approach, 
fijrther processing causes the most important misclassified coefficients to filter to the 
top at a faster rate than misclassified coefficients of lower importance. Effectively bit 
level processing can be considered as a second partitioning stage that deals with 
misclassification. 

There is a tradeoff between the accuracy of the original region partition and the 
additional bit level processing overhead placed on subsequent sorting stages. It is the 
final processing stages that implicitly accomplish the final coefficient partitioning. 
DAC is currently using two bit plane partitioning core technologies for both grayscale 
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and color images, and binary partitioning for bi-level images is now under 
implementation. 

Bit Plane Sorting for Common Mask Miscoverage 

The concept of bit plane sorting is a well know technique used to organize priority in 
a distributed list of data. ID sorting a list of information generates a map specific to 
that list in an optimal fashion. The largest values are mapped first as the list 
transposition routine progresses. The order of decent is the standard order used to 
classify information if considered fi^om the normal processing standpoint where no 
regions are involved at all. As a result, most of the misclassified information is 
considered for further classification and are thus biased to be included first in the final 
code stream order. The only difference is that now there is a region description 
overhead as well as the final mapping overhead of sorting the region of interest 
partitions in a hierarchical fashion. 

Currently DAC has not documented the exact overhead introduced by region 
processing in comparison to normal processing modes that do not use regions. But 
the final results obtained by using the common mask approach are encouraging to say 
the least. At this level the underlying routines do not care where the information to be 
sorted originated. The ordering technique will run to completion whether region 
channels are involved or not. The original misclassified coefficients generated in the 
first cut vnW filter through to the final ordering set by the sorting stage. There is an 
important tradeoff here between the sorting cost and the region overhead that must be 
considered in the final design implementation. 

Bit Plane Ordering for Common Mask Miscoverage 

DAC has developed a quad tree ordering technique (see EQW Ch.4) to track 
coefficient importance in both normal and region based processing modes of 
operation. The normal quad tree approach is to track the leading one position through 
each block under consideration in a recursive fashion. A map is generated in the 
process for each entry based on importance such that decoding follows the same 
order. At this basic level the quad tree order is an efficient method of processing the 
multi-resolution decomposition image information in an effective manner. 
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The tracking or positioning overhead introduced by mapping the data is generally 
smaller than in the ID sorting case. Retaining vertical and horizontal positioning 
information is highly advantageous (in most image processing applications) since the 
underlying technology is generally consistent in this medium. The ID sorting 
5 approach sacrifices the natural vertical correlation that may exit (although it can be 
partially recovered by bit level entropy coding in the sorting routine.) However each 
of the two sorting techniques can be used from an operational standpoint to produce a 
similar net effect. One may produce slightly better results than the next, but the 
control architecture is consistent. Each final ordering routine can be tuned to meet 
3 specific processing needs. The ID sorting approach may be a usefiil candidate for bi- 
level image processing, and DAC is currently investigating this possibility. 
The quad tree approach has currently been modified to operate in arbitrary region 
processing modes. The ID sorting routine tracks the importance in a list of 
information specific to the block under consideration. That is the region partition is 
distributed in ID hierarchical lists of information. The original preordering of the 
coefficients is captured in the common mask procedure. Thus the decoder map 
already exists for the inverse ordering. 

Initially DAC's region processing quad tree begins by building an information map 
for each region channel to guide the classification. The core engine has been 
modified to drops blocks not contributing to the final ordering in each region 
category. Currently the overhead comparisons for each sorting case have not been 
documented. The results will be available fi-om DAC once the region based quad tree 
design implementation is finalized (currently under implementation at DAC.) 
In each case, the misclassified coefficients obtained from the original preorder mask 
partition are considered again for ordering in the sorting stages "that follow. 
Effectively, the most important misclassified coefficients still have the opportunity to 
make important contributions to the final code stream as well as improving the 
reconstrtjction quality measurements in the decoded image. The difference between 
this approach and most efficient classification approaches is that the original 
processing order used of the efficient scheme is somewhat distributed through each 
region channels. However the misclassification in each region channel is forgiven to 
a certain extend the subsequent bit level sorting module. One of the prices that must 
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be paid for any automatic region classification scheme is that the optimal processing 
orders that exist in standard approaches must be considered in a different context. 
The original ordering concepts used for standard procedures must be extended to 
include other ordering techniques specifically tailored for operation in arbitrary region 
processing. 

General Region Processing Concerns 

There is one important hurdle that must be overcome in the design of any arbitrary 
region processing technology. There are times when rigid tile or block control is the 
optimal way to process image information. However, it is difficult to achieve 
arbitrary accuracy in a strict processing order. DAC's ID and 2D sorting techniques 
address this issue in design implementation. In our opinion arbitrary region 
processing channels can be defined in a still image environment. Careful study and 
analysis suggests that this is a viable approach not only for the DCT common mask 
approach but also for other arbitrary region formations techniques that have not yet 
been developed. Furthermore region classification may have interesting properties 
that can transferred to video processing environment where the region of interest 
concept and the object motion vectors can be considered together in a unified 
approach to region of interest processing. 

DAC's region technology in its original form was designed to include region 
processing capabilities from the start. All internal modules can be operated in both 
normal and region based modes of operation. The MUX architecture Cintroduced in 
Chapter V) is designed to meet operational requirements in both normal and region 
processing modes. It is possible to extend the notion of block based image processing 
to a level where a block of information is interpreted in an arbitrary way that depends 
only on the underlying core processing technology. In this way, blocks can be 
considered as individual processing units each fully capable of making a contribution 
to the final code stream. Furthermore control architecture can be developed to fully 
exploit the region processing capability. All region formation schemes, whether 
arbitrarily defined user regions of interest, automatic region formation schemes or any 
other specially designed region classification technique, can be processed in a 
common framework. The same overall processing architecture can be tailored to 
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operate in an optimal fashion for both normal and region processing mod^ of 
operation. 

Additional Processing Concerns 

One of the important concerns in the JPEG-2000 community regarding the 
5 development of a new compression standard is that it must be able to support ROI 
processing: There are a number of techniques that have been tested with some 
success. Some of these techniques are listed below. 

• Sequence based mode. When this mode of operation is used, each set of data 
encapsulated by a specific region is treated as a separate sequence. The sequences 

0 are coded independently. 

• Scaling based mode. In this mode of operation, the magnitude of each ROI 
coefficient is increased using a predetennined shift factor. If multiple ROIs are 
used with increasing levels of importance, multiple shift factors distinguish each 
ROI. The effect of this scaling technique is to force the ROI coefficients to be 

• encoded first. 

• ROI from the middle. In this mode of operation, a concern of coding an ROI in 
the middle of the bit stream is addressed. This is usefiil for progressive 
transmission of images where the user is given the opportunity to focus on a 
specific ROI before the complete set of image data arrives at the decoder. This is 
important for medical image processing. 

• ROI without sending mask information. This mode of operation is a special case 
for scaling based mode. The coefficients for each region are scaled such that the 
magnitude of each ROI exceeds that of the ROI of less importance. The decoder 
knows which ROI the coefficients belong to based on the magnitude. 

The ROI coding technique that DAC has developed generally falls into the sequence 
based category. But internal modules have been tailored to operate in both block 
based and arbitrary block sized processing units. DAC is currently studying the 
different mask and region formation techniques to determine the effect each ordering 
and partitioning approach has on the overall organizational stmcture required for 
region channel processing. The current implementation includes a scalable ROI mode 
of operation implemented without actually shifting the ROI coefficients encompassed 
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by the classification. Data scaling is handled in the MUX control architecture. 
Internally, the scaled version is a special case that can be implemented in MUX 
control. 

Preliminary study has begun for the last of the four categories cited above. ROI from 
the middle implementations can be addressed by using resolution a progressive 
processing order. This type of ordering is supported in the MUX architecture. DAC 
is currently considering a design implementation for client side region of interest 
requests in a resolution progressive architecture for client-server region interaction. 
ROI without mask information makes sense in theory. However it places a huge tax 
on the subsequent sorting and classification stages. Generally speaking, the overhead 
cost introduced by additional sorting is comparable to the size of a bitmap mask that 
could have been used to classify the data initially. The shift introduced to each ROI 
must be processed in the sorting stage. Eventually sorting must traverse from top to 
bottom. Coefficient scaling under tight control has some definite benefits in many 
cases. However it would be highly inefficient in a general region processing 
architecture because of the high sorting costs as a result of traversing all region 
channels. In a region formation scheme of 4 or more regions, the total bit level 
difference between all ROI is too large to process entirely. On the other hand, if the 
resolution of the overall partition is reduced (less bit levels), less room is available to 
form an accurate region classification. There is some common ground that must be 
explored in order to determine the optimal processing mix. 

Intra-Region Coding 

The segregation of image primitives into regions provides a basis for accessing the 
image contents. Compact representation of image primitives is achieved via intra- 
region coding. 

Independent coding units (ICUs) are the building blocks for intra-region coding. In 
RICS, intra-region coding with ICUs is designed with a set of objertives. 
(1) Full data independence. From the notion of ICU, the encoding and decoding of a 
region is done without referring to the data of any other regions. 
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(2) Allowance for multiple coding schemes. Each ICU may be coded using different 
coding schemes. New coding schemes can be added to the system (modular 
openness.) 

(3) High coding efficiency. Depending on the characteristics of the primitives in a 
5 given ICU, one or more efficient coding schemes can be selected to produce a 

code block with high compactness. 

(4) JPEG and JBIG transcodability. It is preferred to have a single JPEG 2000 
coding platform that also accommodates JPEG and JBIG modes of coding. 

(5) Scalability of code stream. 

(6) Error resilience. The ICUs are natural blocks for bit error localization. 

(7) Parallelism. There is no data dependence between ICUs and so the encoding and 
decoding of all ICUs can be performed in parallel. 

Intra-region coding in RICS involves the following steps. 

(1 ) Choosing a coding scheme. 

(2) Determining ICU structure. 

(3) Coding. 

(4) Pre-packing. 

In the rest of this section the first three steps are described. 
Choosing a Coding Scheme 

For coding an ICU, the current RICS design employs six categories of coding 
schemes. 

• Zigzag DCT Coefficients Packing (JPEG baseline scheme) 

• Embedded Quad tree (or similar schemes) 

• Embedded Zero tree (or similar schemes) 

• ID Progressive Sorting Schemes 

• Block Based Quantization 

• JBIG Coding Algorithms 
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There is no general restriction on what schemes have to be used for coding a given 
ICU. However, certain preferred embodiments may apply. For example, if DCT is 
used at the transform step, then the zigzag DCT coefficient packing might be 
preferred (but not restricted to). 

Choosing ICU Structures 

While being a logic concept, an ICU still has a geometric structure. The ICU 
structure is directly relevant to the coding scheme chosen. Currently, the RICS 
recognizes three types of ICU structures (Figure 22). 

Type 1 ICU Structure: 8x8 Blocks 

This type of ICU is designed exclusively for use in JPEG mode, although it is not 
necessary that DCT transform coefficients be coded this v^^ay. 

If there is no region defined in the image and the JPEG mode is chosen, the entire 
image is paved with type-1 ICUs. 

If a region R is to be coded using JPEG mode, the set of ICUs that covers R form a 
unique, minimal pavement for R, in the following procedure. 

Step 1. Let ;;mm be the first row from the top that has at least one pixel inR, the 
last row, Xmin the first column from the left, and Xmw the last column. 
Step 2. Let (xxaji,ymm) and (xm«,>'n«x) be respectively the top-left and the 

bottom-right comers of the bounding rectangle of /?. 

Step 3. Starting from the top left comer, pave the bounding box completely with 
type-1 ICUs. 

Step 4. Remove those ICUs that cover no pixels in R. 

After this procedure, the remaining ICUs form the minimal pavement of R (Figure 
22). 

Type 2 ICU Structure: Rectangles 

Type-2 ICUs are rectangles. Except for the JPEG and embedded zero tree schemes, 
any other intra-region coding scheme uses type-2 ICUs to pave and code a region. 
There is no general restriction on the dimensions of a type-2 ICU. It is usually 
36 



Copied from 105<lJemu?DiEi9»4t!£9 jRfiaitCiS) 



wo 00/49571 



PCT/CAOO/00134 



defined by the selected coding scheme. For example, both embedded quad tree and 
explicit block based quantization can use arbitrary sized rectangles, with some 
preferred embodiments (such as the 64x64 blocks used in EBCOT of VM 3.0 (A)). It 
should be noted that once a region is paved by type-2 ICUs, different coding schemes 
can be used for each ICU. 

For subband decompositions, no type-2 ICU is allowed to cross subband borders. 
Essentially, this means that type-2 ICUs support only various types of intra-band 
coding methods. 

Tvpe-3 ICU Structure: Pyramids 

The third type of ICU is the pyramid strurture, which is used by various inter-band 
coding methods, such as the embedded zero tree. In the current RICS design, there is 
one limitation on the region definition procedure for type-3 ICUs: a region must be 
specified in the top-down manner, fi-om lower resolution to higher resolution in the 
decomposition pyramid. That is, for a given region definition R in LL, every element 
in Jt defines a set of type-3 ICUs (i.e. in three spatial orientations, respertively). 

Let C^ijJ) be a wavelet transform coefficient at spatial location in the LL 
subband, C," {m,n) be a wavelet coefficient at spatial location {m,ri) in the subband 
at level / (1=1,2,...) in orientation r (r='J,2,3). Assuming that a 5-level wavelet - 
decomposition is used, the three type-3 ICUs defined by C^(i,j) can be represented 
in the following fashion. 
In orientation 1: 

[C\{iJ)} U {Ci(2/,2;), C:(2/ + l,'2y), Cj(2/,2y + l), C:(2; + l,27 + l)} U ... 
In orientation 2: 

{CliiJ)] U {C^(2/-,27).C,^(2/ + 1,27). Ci(2/.2y + l). C:(2/ + 1.27 + l)) U ... 
In orientation 3; 

{Cli},])) U {Cl{7j,lj),Cl{2i + \aj\Cl{2i,2j + \),Cl{2i + \,2j-^\)) U ... 

Notice that the primitive C^{i,j) is not included in any the three ICUs. Instead, the 
LL subband is coded differently, using the type-2 ICUs. 
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Coding Tvpe-l ICUs 

Type-l ICUs are defined exclusively for coding 8x8 DCT coefficient blocks. Coding 
of type-l ICUs follows the JPEG baseline algorithms. 

Coding TvDe-2 ICUs 

5 Once a region is decomposed into a set of non-overlapping type-2 ICUs, each ICU is 
coded independently, and the collection of codes for all ICUs is the coded 
representation of that region. In this subsection we describe three techniques that may 
be used for coding type-2 ICUs. 

Using Embedded Quad tree 

10 Embedded Quad tree Wavelet (EQW), Figure 24 is an efficient and fast method for 
type-2 ICU coding. This technique implements an embedded progressive sorting 
scheme in a quad tree-like structure. In contrast to inter-band coding methods, the 
EQW explores the intra-band self-similarity of the wavelet decomposition 
coefficients. The EQW-produced code-stream realizes scalability by pixel-precision 

15 (the scalability by spatial resolution is realized by the multiplexer.) 

Depth-First and Breadth-First Quad trees 

The EQW method can be used for both lossless and lossy coding. In lossless coding, 
the primitives to be encoded in the ICUs are coefficients of a reversible wavelet 
transform. In lossy coding, after the wavelet transform, each transform coefficient is 
20 represented in a fixed-point binary format - typically with less than 16 bits - and is 
treated as an integer. 

In each ICU, the maximum coefficient magnitude M is determined. Then a value N is 
found which satisfies the condition 2" <M < 2"*^ . The EQW works in a bit-plane 
manner: the initial bit-plane is set to 2'', followed by 2''"',2*'"^.. and so on. A 
25 binary significance map is produced for every bit-plane by comparing coefficients in 
power of 2 increments. 

Figure IV.3 illustrates the EQW encoding process on a single bit-plane. Two working 
lists are used. One is called the list of significant primitives (LSP). Each entry in 
LSP corresponds to an individual primitive in the ICU. The LSP is initialized as an 
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empty list. Another list is called the list of insignificant blocks (LIB). Each entry in 
the LIB is registered with the coordinates of the top-left comer of the block, together 
with the width and height information. Each LIB entry represents an individual 
primitive of width and height equal to 1. At the highest bit-plane 2'*', the ICU to be 
encoded is put into the LIB. A temporary list of insignificant blocks (TLIB) is used 
for refreshing the LIB after each bit-plane coding pass is completed. 
There are two methods used to add the four sub-blocks into the LIB that occur as a 
result of the quad tree decomposition. The first method, known as depth-first quad 
tree coding, is to insert the four sub-blocks into LIB at the position of their parent 
block. In this case, the four child blocks are evaluated immediately after the parent 
block. This rule is applied recursively until no more subdivision is possible. This 
means that control returns to the next entry in the LIB only after the present entry is 
completely encoded up to the highest resolution level. The second method, known as 
breadth-first quad tree coding, is to add the four sub-blocks under consideration to the 
end of LIB. In using the breadth-first process, all parent blocks at the same level will 
be processed first, followed by their respective children blocks. 
After all entries in the LIB have been processed, the entries in TLIB can be reordered 
according to certain pre-defined set of rules (this is an optional step.) Experimental 
evidence suggests that a higher PSNR can be achieved by using an effective 
reordering scheme. Then the LIB is replaced with the TLIB and the coding resumes 
at the next bit-plane. 

The next step is the refinement pass. The N*** bit is output for entries in the LSP. 
Then, the process resumes using the new LIB and a new bit-plane set to Error! 
Objects cannot be created from editing field codes.. 

EOW with VLC 

In order to use variable length coding (VLC), the standard EQW algorithm of the 
previous section is modified as follows. 

Step 1 . Initialization: output n satisfying the inequality 2" < max{|Cj, |} < 2"*' , set the 
list of significant primitives (LSP) as an empty list. Put the ICU to be coded into the 
LB. 
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Step 2. Bit-Plane Sorting: for each entry of the LIB, perform quad tree coding. 

If all primitives within the block are insignificant, output a 0 bit and add this block to 

the temporary LIB (TLIB). 

Otherwise, look at the significance of the four sub-blocks. Then according to the 
5 VLC table, output the corresponding VLC sequence (see the table in Figure 25). If 
significant, sub-blocks that are not a single primitive are inserted into the LIB at their 
parent position. If insignificant, they are added to the TLIB. If the sub-blocks are 
single primitives, they are added to the LSP. Sign bits are output only if units are 
significant. If insignificant, the units are added to the TLIB. 
10 Step 3. Reordering (optional). 

Step 4. Refinement: for each entry in the LSP, except those included in the last 
sorting pass, (i.e. with same ri), output the mh most significant bit of |Cj,.| . 

Step 5. Quantization Update: decrement n by 1 and go to step 2. 
Some other researchers have also proposed the embedded coding methods where quad 
15 tree is used to explore the intra-band similarity ""^l However, some major 

differences exist in how the quad tree operates. The work of [2] has suggested that 
empirically the quad tree decomposition should stop at size 1 6x16. However (as a 
result of our experience and experimentation) this size may not always be the optimal 
choice. Even without VLC, the efficiency of DAC's EQW is comparable with EZW. 
20 In fact, for a 2x2 block, if it is insignificant, one zero can represent four zeros. 
However, when a 2x2 block is significant, quad tree coding is inefficient. 

Using Block-Based Quantization 

The embedded quad tree (EQW, EQPC, etc.) is a special case of the more general 
block-based quantization techniques. Other block-based quantization algorithms, 
25 such as EBCOT, can also be very efficient for type-2 ICU coding. 

Using JBIG Alporithm.s 



' F. Kossentini et al. "Embedded Quadtree Predictive Codec", ISO/IEC/JTC1/SC29AVG1/N667. 
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It is noted that in all implicit quantization schemes using progressive bit-plane coding 
methods, including the embedded zero tree and the embedded quad tree approaches, 
at each bit-plane, the ICU coder is actually dealing with a bi-level image (the 
significance map.) Therefore, by adopting those efficient JBIG algorithms as ICU 
5 coders in a JPEG 2000 system, one can realize two modes of interplay with JBIG. In 
the following EQW is used as an example to illustrate this idea. 

In the first mode, the transcode mode, the JBIG ICU coder and EQW ICU coder 
provide a two-way transcode between JPEG 2000 and JBIG (Figure 26). In this 
10 transcode, the EQW procedure is called for bi-level type-2 ICUs, with the total 
number of bit-planes being set to one. Note that the sign bit coding in the EQW 
algorithm should be skipped. 

In the second mode, the embedding mode (Figure 27), the proper JBIG routines are 
called as the bit-plane coder in the EQW routine for certain bit-planes and for certain 
15 decomposition sizes (the quad tree decomposition stops at this particular size and the 
coding is handed over to the JBIG routine). Since the algorithms for bi-level data 
coding are also evolving with the compression technology, having an embedding 
JBIG mode reserved in the JPEG 2000 system can make the Standard evolve with its 
sister technologies. 

20 Codino Tvpe-3 ICUs 

Type-3 ICUs are defined to support various inter-band coding methods for wavelet 
transform coefficients. In particular, the methods of EZW, SPIHT, etc., are known to 
be efficient approaches for coding the type-3 ICUs. In testing these approaches, it is 
noted that the standard versions must be modified slightly to make them fit into the 

25 RICS architecture. Normally, these algorithms implement the ICU coding and 
multiplexing steps into an integrated procedure. A simple modification is required to 
bridge the existing architecture to the RICS architecture. The two steps must be 
separated, which is straightforward. Since the multiplexing stage in RICS can 
effectively simulate the performance of zerotree-based schemes (refer to Ch.V), this 

30 particular split in functionality has virtually no impact on the attainable coding 

Taubman, D. et al. "EBCOT: Embedded Block coding with optimized truncation", ISO/EC/ 
JTC1/SC29/ WGl/ N1020R. 
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efTiciency for standard schemes. Instead it adds more flexibility and openness to the 
existing architecture. 

Intra-Region Coding with Different Typ es ofinu 

Generally there is a need for using different types of ICUs for a intra-region coding. 
5 For example, in coding the wavelet transform coefficients, the LL subband is usually 
encoded differently than the other high-pass subbands. Because different sets of 
primitives may possess different statistical characteristics, providing several coding 
scheme choices in the intra-region coding module may offer a higher coding 
efficiency. In practice, the following combinations are useful. 
10 . In a wavelet decomposition, coding the LL subband with type-2 ICUs whereas the 
other subbands with either type-3 or type-2 ICUs. 

• Tiling an image into several regions, with some of the regions paved with type-1 
ICUs and others with type-2 ICUs. 

• Using different coding schemes in different type-2 ICUs. 

15 Intra-ReeioD Coding without Any TCTk 

Intra-region coding can be performed without specifying, any ICUs in a particular 
region. In this case, the entire region is considered an ICU. If the natural geometric 
shape of the region fits into any of the three ICU categories, the appropriate ICU 
coder can be used. If the region has a very irregular shape, then the ID coding 
20 method can be an efficient approach. Once the data for a particular region is scanned 
in a certain pre-defined order to form a ID stream, the following ID progressive 
sorting algorithms can be used to produce an embedded code-stream. 

ID Progressive Sorting Algorithms 

The EQW algorithm can be readily extended to ID cases where a binary partirion is 
25 used instead of a quad tree decomposition: 

Let L= {c.} be the ID list to be encoded, LSP the list of significant primitives, LIS 
the list of insignificant subsets, TLIS the temporary list of insignificant subsets. 
Step 1. Initialization: Output n satisfying the inequality 2" < max{jc, |} < 2"*' , set the 
LSP as an empty list. Put the set L into LIS. 
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Step 2. Bit-Plane Sorting: For each entry in the LIS, perform the binary partition. 
If all primitives within the subset are insignificant, output a 0 bit and add this subset 
to the TLIS. 

Otherwise, output the two bits that reflect the significance map. For those subsets, if 
5 they are not single primitives, insert them into the LIS at their parent position if they 
are significant, or add them to the TLIS if they are not. If the subsets are single 
primitives, add them to the LSP and output the sign bit if they are significant, or add 
them to the TLIS if they are insignificant. 
Step 3. Reordering (optional). 
10 Step 4. Rennement: for each entry in the LSP, except those included in the last 
sorting pass, (i.e. with same /i), output the nth most significant bit of |c,| ; 

Step 5. Quantization Update: decrement n by 1 and go to step 2. 

Processing with a Multiplexer Architecture 

15 The JPEG-2000 committee has been investigating many emerging compression 
technologies in order to define a new still image compression standard. The emerging 
standard will address both present and near future image compression needs. The task 
of selecting what technologies to include in the new standard is not easy given the rate 
at which technological advances occur in this area. The primary focus of DAC has 

20 been to develop a still image compression engine that addresses the issues set forth by 
the standards committe^. 

The concept of using a multiplexer (MUX) has been adopted in some standardized 
processes of information coding such as MPEG. In contrast, incorporating a MUX as 
an integral part of a still image compression engine is rather unique, and there has 
25 been very little research conducted in this area. However it will be shovim that the 
MUX concept is extremely useful in the development of a general purpose still image 
compression engine. 
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Background 

One of the concerns in developing a new Standard is the layout design of the encoded 
bit stream. A protocol must be developed to guide the new standardized compression 
procedures. Some of the important considerations in this area are listed below. 
5 • A well defined bit stream. All syntactic and semantic aspects of the encoded 
image format must be defined. The JPEG-2000 standard will encompass a wide 
variety of compression needs. One to one relationships will exist between fields 
that exist in the bit stream and functionality that exists in the core compression / 
decompression engine. The number of fields that exist in the bit stream is directly 
) related to the overall functionality encompassed by the new standard. 

• A highly degree of data accessibility. The bit stream must be organized in such a 
manner that it need not be decoded in its entirety to interpret the physical meaning 
of the compressed data. Highly accessible data requires that the encoded data is 
divided into logical units each corresponding to a specific part of the original 
uncompressed information. The way this was accomplished in the existing JPEG 
standard was to partition the original image into 8x8 blocks. Each block was 
compressed using the DCT. Quantization tables were developed for coding the 
coefficients such that the reconstructed image quality was degraded in a standard 
fashion. However, strict block based data accessibility may not be a suitable 
medium for many of the current image processing needs. 

• Dynamic re-orderability. Dynamic re-orderability is a data access issue relating 
to the degree that encoded information can be reorganized to suite user specific 
needs for an image under consideration. This type of data access is an important 
concern for many applications, especially in the medical field. 

• Built in error resilient. In any transmission medium there is always the 
possibility of incurring bit errors. The encoded bit stream must have a certain 
degree of error resilience to address this concern. 

Using a MUX to organize the encoded bit stream can address these issues. The fu-st 
step is to divide the encoded information into logical groups. Organizing the data in 
ICUs leads naturally into a high degree of accessibility. Each ICU is an independent 
unit. The size and position information is self-contained within each ICU. One of the 
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effects of using a muliplexed bit stream design is that it is inherently error resilient. If 
any logical unit incurs an error, only that unit needs to be recovered. If the encoded 
data is divided into an array of logical units in preparation for MUX encoding, 
dynamic re-ordering is much easier to achieve. 

MUX Design Overview 

There are five important considerations to address in the design of the MUX for still 
image compression systems. 

• Defining a working model. 

• Defining an optimal data processing order. 

• Capturing the data to be compressed in a suitable data structure. 

• Defining a data priority to be used for compression. 

• Packing the compressed data based on data priority and processing order. 

A Working Model for the MUX Discussion 

In preparation for the MUX discussion, consider that a multi-resolution 
decomposition hierarchy of data has been obtained for an image to be compressed. 
Also consider that there are three channels of data. Assume that the original image 
has been transformed from RGB to YUV for example. In addition to this, assume that 
the U and V channels have been down sampled by a factor of 4:1. This type of 
approach is common for lossy image compression. The data sets appear in Figure 18. 
The diagram illustrates a 4 level multi-resolution decomposition for the Y-channel 
and 3-levels for each of the UA'-channels. The wavelet transform MALLAT 
decomposition is used here for convenience. 

This is by no means the only model that exists for the MUX design. DAC makes the 
distinction between region and non-region processing modes of operation in their 
design implementation. A hierarchy of list structures is used to control optimal bit 
budget distribution and flow of both non-region and region bit levels through MUX 
channels. The discussion presented in. the next Section begins by introducing the 
concepts in normal processing modes. The normal MUX modes of operation are 
generalized later on to include specialized region and mixed processing modes. 
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Non-Region Data Processing Orders for the MUX 

Generally speaking, given a wavelet multi-resolution decomposition of data, 
coefficients of a fixed magnitude contained in a lower resolution level are more 
important to the image reconstruction procedure than coefficients of a similar 
5 magnitude at a higher resolution level. This leads to a natural ordering of the 
coefficients beginning at the lowest resolution level. The natural ordering concept is 
an important consideration in the design of the MUX. 

General Color Processing Order (Lossy Caset 

Given the multi-resolution data sets of Figure V.l, a natural processing order exists 
10 for each decomposition channel. This natural order is illustrated in Figure 29. This 
particular order reflects a level priority inherent to each data set. 

Another two orders are obtained by simply interchanging orientation data sets in a 
given level. The difference lies only in the convention used for ordering the detail 
information (i.e. LH, HL, HH). The 4-1-1 do-wn sampling relationship that exists in 
15 the color transform domain also exists in the wavelet decomposition data since there 
is one less transform level for the UA^ channels. This is one of the keen design 
considerations that can be exploited in many stages of the MUX implementation and 
encoding / decoding procedures. The data is organized into list structures to cover the 
natural processing orders that exist in the data to be compressed. 

20 Level Priority Processing Orders (Lossy Case) 

The individual orders of Figure V.2 can be combined into a single natural processing 
order for the example under consideration. The new order is obtained by interleaving 
the expressions of Figure 29 using the inherent decomposition level priority of the 
data. The result is illustrated in Figure 30. 

25 Notice that in this particular order, the 4-1-1 color down sampling relationship is 
maintained. In terms of the reconstruction process, the complete level 4 Y-channel 
information (i.e. LL4-HL4-LIL,-HHL,) is required together with the U and V channel 
low pass information (i.e. LL3 in each case). This is the information the decoder wall 
need first to reconstruct the low pass inverse color transform image at the lowest 

30 inverse wavelet resolution level. A "level split" on the U / V data channels at the 
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lowest inverse wavelet resolution level is used to exploit the natural relationships that 
exist for both inverse transform spaces. A similar order exists for the lossless case 
where flill data sets are used. 

Color Interleave Processing O rder (Lossy ) 
5 Another processing order that exists for color images in the wavelet transform domain 
is obtained by interleaving the color channel detail information. The order is 
illustrated in Figure 31. This type of ordering may be suitable for certain applications 
such as those that require the data to be encoded for a progressive download. As in 
the previous case, the down sampling relationships that exist are maintained in both 
10 color / wavelet transform spaces. 

Lossless Color Processing Orders 

As- in the lossy case, similar processing orders can be designed for lossless color 
image compression where full data sets exist for all color channels. In this case, 
corresponding orientation information from each wavelet channel is related in the 

15 original inverse color space. These data should be processed and grouped together in 
the wavelet transform domain and eventually the final bit stream. 
The natural processing order for lossless color image processing is illustrated in 
Figure 32. The color interleave order for lossless image processing is illustrated in 
Figure 33. Note that in each case, the full decomposition data set for each channel is 

20 maintained. The inherent relationships that exist in the color transform space are 
maintained in the wavelet transform space for ordering and reconstrurtion of the 
original image. 

Capturing the Data in a MUX List Structure 

There are a number of candidate algorithms under consideration by the IPEG-2000 
25 committee for optimal quantization of the wavelet coefficients. A common approach 
used to quantize the wavelet coefficients in a progressive manner is to use an optimal 
bit plane ordering technique. DAG has testing a number of bit plane ordering 
techniques in both one and two dimensional schemes. Bit plane ordering techniques 
are generally classified as using implicit quantization for the wavelet data sets. 
30 Coefficients are packed into the final bit stream based on bit level priority. The bit 
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plane processing technology will be used to introduce the data structure used inShe 
MUX design. However, the MUX processing discussion that follows is not restricted 
to bit level processing algorithms. 

The general ordering and processing relationships are useful for introducing the MUX 
5 control architecture at a basic level. However, the key to the operation and many 
benefits of the MUX technique is in the data structure design used to encompass the 
natural ordering concepts. The first step is to define a list structure for each level, 
channel, and orientation of data that exists in the wavelet transform space. A typical 
example of this type of structure is given in Figure 34. 

A general list structure for each level data set can be utilized in many ways. A 
lossless "prepack" of multi-resolution hierarchy information is useful for exploiting 
the many relationships that exist in the data taken in MUX list processing context. 
The data and information contained in each list is used to control how the final bit 
stream is organized. All data packed into the final bit stream must conform to this 
structure. New fields may be added, but the basic operation of the structure remains 
the same. Most implicit quantization techniques are implicitly decodable. In other 
words, minimal header information is required for each list. The MUX architecture 
organizes the lists or units of information into a bit stream that is scaleable in terms of 
bit precision and resolution and is controllable by the many MUX modes of operation. 
As an example, suppose that EQW (DAC's current two dimensional bit level 
processing design) is used to capture the leading ones and refinement bits at each bit 
level in each data set. Following the order given in Figure V.3, a multi-dimensional 
hierarchy of list structures is defined for use in the ensuing data processing stages. 
Within each list (as EQW progresses), the fields in the corresponding list structure are 
updated. The bit stream at each bit level is appended to the MUX "prepack" buffer as 
well as its corresponding size being added to the total in the bit packing information 
field. The bit level position where processing begins is put in the high bit information 
field. The total packed list size is placed in the total packed field. Finally if more 
than one internal processing scheme is used (e.g. NULL transform mode for bi-level 
information or other internal modes of operation), the scheme used for the list is 
placed in the scheme field. 



Copied from lC^I6W&fiW7(W-RR9£ 



wo 00/49571 



PCT/CAOO/00134 



Each MUX list is fully independent, implicitly contains a full set of statistical 
information available for determining the amount of data to pack for each list and is 
optimal for organizing the final compressed bit stream. 

5 Data Priority and MUX Control 

The information contained in each list structure is used to organize the final bit stream 
based on end user compression requirements. DAC has implemented numerous MUX 
modes o operation. Two common modes of operation are outlined in this Section. 

■<;i >nai to Noise Ratio (Sm.) Progressive Mode 

10 This mode of MUX places priority on the data such that the final ordering optimizes 
the SNR given a user specified image compression size. For the discussion that 
follows, assume that the user has specified the desired compression ratio in some 
manner. Assume that an optimal bit budget distribution scheme is in place. Thus the 
number of bytes that has been allocated to each color channel is known (more will be 

15 said on how to determine the optimal bit budget distribution for each channel using 
the MUX later). The information fields in the MUX list structures are used to 
determine the amount of data to pack for each list such that the SNR is optimized for 
the specified bit budget. A psuedocode explanation of this technique appears in 
Figure 35. 

20 Fields that refer to the MUX list structure appear in italics while other local variables 
appear in normal typeface. Basically the idea is to loop through each list checking to 
see whether the current processing bit level is equal to the bit level of the list under 
consideration. If the two bit levels are equal, the channel bit budget is decreased by 
the amount of data available in this particular list for the bit level in question (from 

25 MUX field pliPackinglnfofcCurBitLevelJ. The MUX field liCurBytesCount is 
incremented by the number of bytes available. Additional processing takes care of the 
remaining bits. The remaining bits will be considered in the channel bit budget on the 
next visit to this particular list. The MUX field cCurBitLevel is decreased by one 
such that it will again be enabled when the bit level in the main loop is decreased by 

30 one. 

49 



Copied from 103^P^1WW^W&^6) 



wo 00/4957] 



PCT/CA00/0ni34 



This procedure optimizes the SNR since it processes the largest bit levels in each list 
first. Thus information for the largest coefficients throughout individual 
decomposition data sets is sent first according the natural processing order outlined 
earlier. Once the channel loop exits, each MUX list will contain a field indicating the 
5 amount of data to pack into the final bit stream in each case. 

Resolution Progressive Mode 

The resolution progressive mode of bit budget distribution is a simple extension of the 
psuedocode implementation of the SNR mode of the previous Section. In order to 
obtain a resolution progressive mode, the 'Waveket Transform level' for loop is taken 
10 outside of the main BitBudget / BitPlane while loop such that resolution level is given 
priority. In this manner, a lower resolution image can be reconstructed in the inverse 
wavelet transform stage. The resolution of the ensuing image depends upon the 
number of levels required in the end user requirements. 

Data Packing Using the MUX Lists 

15 Composing the final bit stream is a simple matter once the bit budget distribution 
scheme has run to completion. The packing technique can follow one of the normal 
processing orders outlined eariier (e.g. color level priority or level priority color 
interleave) or another access technique as required by the end user. The amount of 
data to pack for each list has been determined in the bit budget distribution stage. If 
20 the data to be packed for a given list is greater than zero, then the total data size, 
scheme and high processing bit level are packed as header information into the final 
bit stream. If the data to be packed for a given list is zero (i.e. not required), then a 
smaller header is packed to indicate a zero length list. Currently a register type 
approach is being developed to handle zero length lists. All lists that have data 
25 available will set a bit position in the list register so that the decoder can determine 
which lists have made a contribution. Empty list will be flagged with a 0 bit in the list 
register. The list register is packed in the final code-stream. This technique will save 
the size of the list header overhead of empty lists. 
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Overhead Consideration of the MUX Architecture 

TTiere is a small overhead placed the compressed data by the MUX implementation. 
There are three header fields to pack for each MUX list. 

• The total bytes packed. 

5 • The scheme used to process the list. 

• The highest bit plane for data in list. 

Note that if there is only one internal scheme used to process the entire image, then 
the scheme field does not need to be packed for each list. 

In order to calculate the cost of using a MUX in terms of header sizes for each list, 
10 consider that a 24 bit color image of size 2048 by 2048 needs to be processed in a 
lossy fashion (i.e. the working MUX model of Figure 28). Assume that the size of the 
smallest wavelet transform level is 8 by 8. Thus there are 8 decomposition levels to 
consider in forming the MUX list hierarchy. Assume that the wavelet decomposition 
data has been converted to a 16 bit integer representation. Also assume that there is- 
15 only one internal processing scheme for the lists. Given these initial conditions, a 
worst case analysis is conducted to determine the header size for each list. The results 
are tabulated in Table V.l for the Y-Channel wavelet decomposition. 

From table in Figure 36, the total Y-channel header size for this particular image is 
443 bits. A similar calculation can be conducted on the UA^-Channels (assuming one 

20 less decomposition level for each) to yield 368 bits each. The original raw image size 
is 2048 X 2048 x 3 bytes. Thus the MUX packing overhead in this particular example 
is about 1 header bit for every 854 bits of compressed data. In the lossless case where 
each channels has the same number of decomposition levels (or correspondingly, in 
the 8 bit grayscale case), the overhead is still very small at 1 header bit for every 757 

25 bits of compressed data packed. 

The table in Figure 37 outlines the overhead size for square images with dimensions 
that are a power of 2 beginning at 16 by 16 and ending at 64k by 64k. Note that the 
overhead size if 44 bits in both the lossy and gray scale cases. The reason for this is 
that the wavelet transform is not used on the UA^ channels when the original image 
30 dimensions are 16 by 16. The Y-channels is decomposed once. A plot of percentage 
overhead versus image dimension is given in Figure 38. Both lossy and lossless cases 
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appear on the same plot. Note that, the image size can be quite small while still 
maintaining a relatively low overhead in terms of the total MUX list header sizes. 
Bit Budget control from the MUX Architecture 

One of the most difficult tasks that must be addressed in any compression scheme is 
scalability. The internal procedures must be developed such that the image under 
consideration can be compressed to a user specified size in an optimal manner. The 
task is compounded for color images since there are three channels to consider. The 
major problem is how to distribute the user specified bit budget between the three 
channels. 

Standard color transforms such as YUV or YIQ are normally used to redistribute the 
RGB color information prior to the multi-resolution decomposition stage of the 
compression procedure. These transformations concentrate most of the energy into 
one channel making them better suited for compression. Another color transform that 
is gaining popularity is the KL transform. This particular transform technique is 
based upon a principle component analysis (PCA) implemented on the original raw 
color information in order to determine the optimal color redistribution for the three 
RGB channels. Generally speaking, after completing a color transform stage, two of 
the three channels are down sampled by a factor of 4 to 1 for lossy compression. In 
doing this, the amount of data that must be compressed is reduced by a factor of 2 
with minimal loss in visual quality for the reconstructed image. 

The easiest bit budget distribution scheme to implement is one that follows the color 
transform dovm sampling ratios. If the transform under consideration is yUV-411, 
then 4 parts are allocated to the Y-channel for every 1 part allocated to the U. / V- 
channels respectively. However, the energy distribution of the wavelet transform data 
generally does not follow this strict down sampling guideline. Instead it varies from 
one image to the next. 

DAC has developed a dynamic bit budget control architecture that is based on the 
implicit information contained in each MUX list. In this Section, a bit budget 
allocation procedure will be introduced that can be used for optimal scalability in 
normal processing modes of operation. This technique is based implicitly on the 
amount of energy contained in each color channel of the wavelet decomposition data 
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sets. A data ratio concept is used together with the prepack information contained in 
the MUX list structures to determine the budget for each channel. 

Compression procedures that process wavelet coefficients in a bit plane order are 
implicitly tracking the energy contained in the decomposition. The most significant 
bits (MSBs) of the largest coefficients are packed first followed by refinement bits 
and the MSBs that are significant for coefficients at the next bit level. This process is 
repeated in a recursive fashion until the desired compression size is obtained. 

Let yij , Ujj and vy , be the total data available in a particular orientation level for the 
three color channels. TTien the total amount of data available in the Y-channel 
decomposition Yt is obtained by summing the individual totals. A similar procedure 
can be followed to obtain Ui and Vt. 



yN4 + S 2, yu 

N 3 

UN4 + X 2 "u (Eq.V.l) 



^ VN4 + 2 S Vii 



The next step is to calculate the pack ratios to be used for each channel of the wavelet 
decomposition hierarchy. The pack ratios are determined by taking the ratio of the 
two largest amounts of data to the smallest amount of data. Let the three pack ratios 
be denoted Ry , Ru and Rv. The smallest data size will be in either the U or V 
30 channels because of the down sampling step. Note that one of the pack ratios will be 
unity. 

Y, 

35 . R = 

" MIN(U,,VO 

. V, 
~ MIN(U,.VO 
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The pack ratios given in Eq.V.2 are used to determine the optimal amount of data to 
allocate for each color channel based on the user specified compressed file size. 
However, any additional overhead introduced into the final bit stream by the MUX 
must be taken into consideration. Let yhij, uhij and vhjj be the individual header sizes 
in each wavelet channel for a particular orientation level. Then the total header sizes 
for each wavelet channel are obtained by summing the individual totals. 

N 3 

1=1 ri 

N 3 

Uk = uhK4 + g 2, "hii (Eq.V.3) 

N 3 

v. = vhw4 + X 2 vhii 
1=1 j=l 



If doing ROI processing, the packing overhead introduced by the mask must be taken 
into consideration in determining the bit budget for each channel. Let the total pack 
size for the ROI mask be denoted as M,. The logical approach to take is to make an 
adjustment in the bit budget for each channel based upon the pack ratios of Eq.V.2. 
In this manner, the mask overhead is distributed proportionally to each wavelet 
channel. The unit adjustment factor U. is determined from the total mask overhead 
and the total pack ratio. 

Mt 

U. = (Eq.V.4) 

Ry + Ru + Rv ^ ' 



The adjustment size for each channel is determined by using the result of Eq.V.4 and 
the pack ratios of Eq.V.2. 

= Ry • u. 
U> = Ru-U. (Eq.V.5) 
V' = RvU. 54 
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A useful set of calculations that comes out of Eq.V.l to Eq.V.5 is the minimum 
compression bits per pixel (BppMin) and the maximum compression bits per pixel 
(BppMax). These values are calculated below (assume that the original raw image file 
size is Fi bits). 

BppM» = — • ( Y. + Ut + Vi + Yb + Ub + Vb + Y. + U. + V.) 

8 s (Eq.V.6) 

Bppwi. =— • (Yb + Ub+ Vb + Y. + U. + V.) VM / 



These expression represent bounding compression values available for the image 
20 under consideration. If all data is packed for each of the wavelet data sets, the result 
is BppMax. If just the header information is packed, then the result is BppMin. 
Suppose the user specifies a compressed file size of The unit bit budget Uhb is 
determined from the pack ratios and the total header sizes. 

25 

Fi-(Yh + Uh + Vh) 

Ubb = * — ^ (Eq.V.7) 

Ry + Ru+Rv 

30 

Now the optimal bit budget for each channel Ybb, Ubb and Vbb can be determined 
using Ubb, the individual pack ratios and the overhead parameters. 

Ybb = Ubb . Ry + Yb- Y. 

Ubb = Ubb. Ru +Uh-U. (Eq.V.8) 

Vbb = Vbb. Rv + Vh- V. 

40 

The expressions given in Eq.V.8 will always yield a near optimal bit budget for each 
color decomposition channel. The technique is based implicitly on the energy 
distribution map recorded in MUX lists and the total overhead associated with the 
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underlying process. Note that if processing with no ROI, then that variable falls out 
of the calculation. 

This concept can be extended to form pack ratios for each subband if more accuracy 
is required. A similar data ratio technique can be implemented to control the bit 
5 budget in that case. The necessary information is contained in the MUX list 
strurtures. However for most practical application, the cited distribution technique is 
sufficiently accurate. The user specified file size can be obtained within several bytes 
in implementing this technique. In addition, the maximum and minimum attainable 
file sizes are known prior to the final packing stage. 

10 A simple routine has been developed for distributing the final bit budget in the level 
where it is determined it will expire. The idea is quite simple. Given that the bit 
budget will expire on a particular bit plane in a certain transform level, then the bit 
budget is redistributed such that each orientation gets a proportional amount based on 
the amount of data that each orientation can take. The amount of data each 

15 orientation gets is determined by using ratios between orientations on the bit level 
under consideration. This is an important consideration since the net effect is 
approximately a 1 dB improvement in the PSNR. 

Using the MUX for Region Processing 

So far the focus has been on the design of a MUX control architerture for normal 
20 image processing modes of operation. The main highlights are the data processing 
orders, the MUX list structure, the bit budget control technique, and organizing the 
final bit stream. DAC has extended the MUX concept to include a variety region 
processing modes of operation. 

Region Processing Orders 

25 The processing orders introduced in the previous Section are specific to the normal 
processing modes where one or more ROI are not specified. DAC has implemented a 
region processing design based on the natural extension of the MUX concepts 
introduced in the previous Sections. Both lossy and lossless processing orders are 
again considered. 
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General Color Region Level Processing Order 

The region level processing order places priority on each ROI in a descending fashion 
followed by the normal level priority scheme of the original MUX design. Inherent in 
the MUX scheme is the general assumption that data of similar magnitude in each 
5 succeeding region is less important to the overall image reconstruction. A secondary 
priority key is placed on the wavelet transform level. Data of similar magnitude 
contained in a lower resolution level is more important to the reconstruction 
procedure than data at a higher resolution level. The general processing order for 
each color channel (assuming a 4 region priority scheme) is given in Figure 39 where 
10 Ri is the region level for i = 1,2,3,4. 

Region Level Color Processing Order (Lossy') 

The individual orders of Figure V.IO can be combined into a single processing order. 
The new order is obtained by interleaving the expressions making use of the inverse 
color and wavelet transformation relationship that is presented for the normal MUX 
15 modes of operation above. The result is illustrated in Figure 40. 

The intrinsic 4-1-1 color down sampling relationship is maintained for region 
processing MUX modes. In terms of the reconstruaion process, the complete level 4 
Y-channel information (i.e. LU-HI^-LHi-HHU) is required together with the U and V 
channel low pass information (i.e. LL3 in each case) for each region under 
20 consideration. 

Color Interieave Region Level Processing Order (Lossy') 

The color interleave processing order extends naturally to the MUX region processing 
modes. The order is illustrated in Figure 41. This type of ordering may be more 
suitable for certain applications that require the data to be encoded for a progressive 
25 download based on regions of importance. As in the previous case, the doAvn 
sampling relationships that exist for the inverse color and the wavelet space are 
maintained. 

Region Color Processing Orders HLossless) 

As in the lossy color region processing MUX modes, similar processing orders can be 

30 outlined for the lossless case where full data sets exist for all color channels in each 
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region. In this case, corresponding decomposition level orientation information in 
each channel is related in the original color transform space for each region under 
consideration. These data are processed and grouped together in the wavelet 
transform space and are ultimately put in the final bit stream by the MUX control 
5 architecture. 

The processing order for lossless color region processing is illustrated in Figure 42. 
Note that in each case, the full decomposition data set for each channel is maintained. 
The inherent relationships that exist in the inverse color and wavelet space are 
maintained for each region channel. 

10 DAC has developed both SNR progressive and resolution progressive MUX modes of 
operation for region processing. In addition, several specialized MUX modes were 
developed that may be useful in certain applications. The regions can be defined 
automatically by the technique described in Chapter in, or user defined ROI can be 
used to group the data. User defined ROI can have arbitrary shape or can be defined 
. 15 in terms of simple geometric elliptical or rectangular region primitives. 

One of the specialized MUX region processing modes can be introduced at this point 
since it fits into the context of the ordering discussion. 

Transparent Region Color Processing Orders fLcs-sy) 

This particular region color processing order is useful for applications that require 
20 transparent region channels that has little or no influence on the ordering of the data. 
Either of the normal SNR and resolution progressive modes are overlaid with a 
complete region channel description. The technique is implemented simply by taking 
the region index to the inner processing loop in the pseudo code implementation 
example of Figure V.8. This ordering technique may be usefiil in video processing 
25 applications where a complete region mask description is required together with a 
compressed fi-ame of information. The region description may be used to process 
subsequent fi-ames in the sequence. One of the processing orders that falls into this 
MUX mode of operation is illustrated in Figure 43. In this case, the high bit plane of 
each region list is processed first causing the region classification of all coefTicients to 
30 be transparent to the distribution and packing routines. 
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Transparent Region Color Processing Orders (Lossless) 

Corresponding transparent region processing modes exist for the lossless case. A 
typical example is given in Figure 44. As in the lossy case, the final bit stream is 
organized independent of the constraints of the region channels which may be useful 
5 in certain cases. 

MT7X List Structure for Region Processing 

The MUX organizational list structure concepts introduced for normal processing 
modes have been extended to include many region processing modes of operation. In 
this case, the total number of lists is a function of the number of ROI. If there are 4 
10 ROI, then there are 4 times as many MUX lists. However, the basic operation of the 
MUX control architecture is similar in each case: 

There are three high level modes used to categorize the MUX architecture. 

• Normal processing mode (ROI disabled). 

• Region processing mode (ROI only). 

15 • Mixed processing mode (normal and ROI enabled). 

The next Section outlines the operation of the bit budget and MUX controls for many 
of the region processing modes developed at DAC. The mixed processing modes are 
briefly discussed later in the Chapter. 

Bit Budget Control for Region Processing MUX Modes 

20 Transparent Region Level Color Mode 

Transparent region processing is mentioned above. In that particular ordering 
example, the region index is placed in the inner most loop in the bit budget 
distribution function. The processing index placement is basically a method of 
incorporating a transparent region layer of processing into the normal MUX modes 

25 outlined earlier in the Chapter. Processing and packing in this mode gives ROI a low 
priority. 

Distributing in this context ensures that the SNR progressive and resolution 
progressive MUX modes operate as they did before. The exception in this case is that 
there is a variable length region mask overhead that must be taken into consideration. 
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The bit budget distribution functions are calculated as they were for the normal hfex 
processing modes. Based on the file size (or resolution level) requirements, each 
color channel receives a proportional amount of the bit budget based on the ratio 
distribution technique used in the MUX architecture. 
5 The DCT auto regions and arbitrary or primitives user defined region types can be 
operated in this mode. In DAC's internal processing architecture the user must select 
the region coverage technique used to categorize the wavelet coefficients. DAC has 
implemented the wavelet mask down sampling technique of the JPEG-2000 VM, and 
it can be used in any of the region processing MUX modes. In this manner a common 
10 mask can be used in each orientation level or the VM mask down sampling technique 
can be used for individual orientation level mask coverage, for both automatic or user 
defined ROI in lossy and lossless MUX modes of operation. The number of region 
categories can be selected as 2, 3, or 4 channels. 

In order to test the transparent mode of operation the YUV-411 color transform is 
15 applied to a 24 bit 256 x 256 Glacier park mountain scenery image. The 9-7 kernel 
implemented in a lifting scheme is used as the wavelet transform. The DCT region 
formation technique of Chapter in is employed to generate the common masks for 
each wavelet transform level. The mean square error (MSE) is measured for a 
number of compression ratios. The result is illustrated in Figure 45 in a plot of MSE 
20 versus Bpp. The plot shows the break points in MSE for each automatic ROI used in 
the example. The rate of image degradation (i.e. MSE and PSNR) must be controlled 
precisely for each ROI. 

The MUX overhead is calculated based on the mask type and process selections. For 
25 arbitrary user constructed masks, the mask file can be loaded to guide the MUX. In 
the 4 ROI common mask case, the overhead is 0.5 Bpp in packing the entire mask. In 
the VM mask case, the overhead is 2.0 Bpp. If simple primitives are used to form the 
user mask instead of an arbitrary shape, the overhead is greatly reduced. The arbitrary 
mask overhead can be reduced with the addition of an entropy coding stage. The 
30 common DCT mask overhead is cited in Chapter m as approximately 1Kb (0.2 Bpp) 
for an 8 bit 256 x 256 image size. These are rough estimates that do not take the 
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header sizes into account. However, they do serve as a comparative guideline for the 
current discussion. 

Region Priority Level Color Mode 

In this mode of MUX operation, the pack ratios for each wavelet channel are 
5 determined as in the normal processing mode. However, there is one important 
difference. In this particular mode of operation pack ratios are determined for each 
region channel. This implies that there are 12 pack ratios for a 4 ROI process. The 
bit budget distribution is based on the overall wavelet channel pack ratios and the ROI 
data totals in each region channel. The distribution function also takes the mask and 

10 list header size information into account in the overall calculation. 

This mode of MUX operation is designed to distribute the MSE and the PSNR image 
reconstruction measurements in an approximately uniform manner for all region 
channels. A proportional amount of overall bit budget is allocated to each region 
channel based on the region and color channel pack ratios. Figure 46 illustrates the 

15 result of using this mode of operation for the Glacier park image. Notice how the 
quality of each region channel degrades in comparison to the others. Auto detected 
DCT common masks are used to generate the result. 

Absolute Region Priority Level Color Mode 

20 In this mode of operation, absolute priority is given to the region channels. The 
distribution function is the same in this case in that the bit budget is divided into 12 
according to total data ratio technique mentioned earlier. However the bit budget is 
distributed in a biased fashion giving the highest priority to the most important ROI. 
Some regions may not receive a budget based on the compression requirements. The 

25 bit budget is distributed region by region beginning at the most important ROI. Thus 
complete sections of the original image can be eliminated altogether if the bit budget 
is small. If only the important regions of an image need to be saved, this technique 
can be employed to partition the image. 

A plot of MSE versus Bpp for the absolute region priority color MUX mode is given 
30 in Figure 47 for the Glacier park image. Notice how the regions fade out very fast 
and sharply. This occurs when the region no longer has any budget allocated to it. At 
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that point, the region is basically invalidated in terms of any contribution to the 
reconstructed image. That portion of the image basically fades out. Note that the 
amount of quality of the fade out depends on the down sampling technique to translate 
the masks. The VM mask formation technique causes a very gradual fade out to 
occur. In the common mask approach, the effect is much more abrupt. 

Scaled Region Priority Level Color Mode 

In this mode of MUX operation degradation rate of each ROI can be controlled. Each 
ROT contributes to the final bit stream. However, instead of having each ROI degrade 
at a uniform rate (as Section V.5.3.2.), the quality measurements of each succeeding 
can be controlled. The initial bit budget for each region channel is calculated as 
before in that the specified compressed file size is split into 12 bit budgets to be 
spread between each of the 4 region channels. However after the initial split, the 
allocation amounts are changed heuristically based on a priority factor that is set for 
each ROI. For example, suppose it is decided that ROI 1 is 50% more important than 
ROI 2, ROI 2 is 30% more important than ROI 3 and ROI 3 is 20% more important 
than ROI 4 (the background). Using this assumption as a starting point, the amount of 
data allocated to each region channel is changed slightly. The net effect is to cause 
the quality measurements in each region channel to degrade at slightly different rates. 
A plot illustrating this effect is given in Figure 48. Note how the MSE break points 
for regions 2 and 3 have shifted slightly towards the left. 

The result of using this particular MUX mode illustrates an important property that 
should be available for any ROI processing technique. The problem with many ROI 
processing techniques is that it is difTicult to control how much each region 
contributes to the final bit stream. The technique outlined here can be used not only 
to implement this eflFect, but also to control the degradation for each region channel 
based solely on an importance factor associated with each ROI. 

Region Percentage Prioritv Level Color Mode 

In this mode of operation the user can specify a data percentage for each region based 
on the total amount available in each region channel. The distribution fijnction is 
slightly different in this case. The data totals are determined for each region channel 
along with the wavelet color channel totals. There are still 3 bit budget for each 
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region channel. In this case however, the amount of data to include for each region 
channel can be set by the user as a percentage of the total in each case. 
Currently when this mode, of operation is invoked in DAC's compression engine, the 
user can cycle through the operation until the desired size for each region is obtained. 
The data totals are displayed after each nin. Note also that in using this mode of 
operation, one or more ROI can be eliminated or faded heavily as in the absolute 
region priority MUX mode of Section V.5.3.3. 

A plot illustrating the result of using this mode of MUX operation is given in Figure 
49. Initially all available data is included for each region followed by equal 
decrements for each region channel. 

Using the MUX for Mixed Processin|g 

So far both normal (non-region) and ROI processing modes have been discussed. In 
addition to these modes of MUX operation, DAC has developed mixed mode 
processing capabilities. The modes of operation discussed for both normal and ROI 
processing are extended such that they can be run simultaneously for an image under 
consideration. A wavelet transform level partition is conducted based on the desired 
number of region and the number of non-region levels. 

Currently non-region levels can be defined for. processing lower resolution levels and 
region levels can be defined for process higher resolution levels. However, the 
implementation is not restricted by this distinction. It can be changed to regions over 
non-regions or to an interlaced combination of the two. The parameter that controls 
this distinction is termed the region start level. It can be set to any valid v/avelet 
transform level (it is set internally to -1 to disable region processing.) Thus the MUX 
technology can be used in exclusive non-region mode, exclusive region mode or a 
combination mixed mode of operation. 

The ordering techniques outlined for normal and region processing MUX modes will 
not be duplicated in this Section. The same concepts apply for mixed mode 
processing. The bit budget distribution functions work as they did before, but in this 
case they exist simultaneously. In some modes of operation there may be as many as 
1 5 bit budget definitions used to organize the final bit stream. The method used to 
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determine them has not changes. The same ratio technique that has its basis in the 
MUX list structure is used to determine the appropriate allocation in each case. 
There is one additional feature that can be exploited for mixed mode processing. An 
importance factor can be attached to the non-region levels. The net effect of this 
5 parameter is to taper the amount of data included for the non-region levels. In the 
current implementation, the non-region levels are processed first. Thus they 
considered first by the distribution function. The importance factor allows the user to 
decrease the amount of data included for the non-region levels by a certain 
percentage. The delta amount is considered in the bit budget distribution for the 
10 region levels. 

The region processing overhead is slightly smaller in this case. Depending on the 
region start level parameter, there will be less region header information required in 
the bit stream since there are fewer region levels. However, the header overhead for 
the lower resolution levels is quite small anyway. 

15 One example is given here to illustrate the operation of mixed mode processing (see 
Figure 50). In this particular case SNR progressive mode is used for the upper 3 
levels and scaled region priority on the bottom 3 levels. Thus the bottom level 
degradation rates have been adjusted to favor the most important region and 
attenuating in some fashion between the other 3 regions. Notice how the MSE is 

20 lower in this case with the same region overhead as the previous case of Figure 48. 

DCT Region Fo rmation as a Classification Scheme 

In observing this result, it is apparent that for the same region overhead, a better result 
is obtained in mixed mode. There are a number of reasons for this. The first is that 
the region channel coverage is not an exact overlay. All masks used to group or 

25 categorize data must deal with the down sampling issue at different resolution levels, 
and the tight overhead restraints of the compression channel. As it appears in the 
results presented here the DCT region detection mode can be applied to threshold 
wavelet coefficients to form region channels. The amount of data packed for each 
channel can be controlled by the MUX. In addition the initial data split used to create 

30 the partition can be controlled before the DCT mask procedure begins. 
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There is some region migration or intermingling of the coefficients that occurs at the 
boundaries where sharp changes occur in the original image. Some work was 
conducted in developing heuristic techniques to decrease the miscoverage in hot 
areas. And there is a benefit there. Furthermore the DCT low pass filtering stage 

5 affects the original priority bias used to partition the data. There will be a benefit in 
determining the optimal data split. By adjusting the original region partition (e.g. 
equally spaced regions), the plot of MSE versus Bpp can be set to partition the data in 
other ways. Inter-resolution coverage is another problem to consider. The mask 
generation technique of the JPEG-2000 VM addresses the miscoverage problem. 

10 DAC has completed some initial investigations of the VM down sampling technique 
and further testing is required. 

There is a trade off between accuracy of each region channel and mask overhead 
required for region channel operation and generation of the accuracy in the first place. 
The DCT approach shows much promise especially on the highest resolution level 
15 where the masks are the most accurate. One of the benefits of using this the DCT 
approach is that it can be translated to other transform levels in the frequency domain. 
Or alternatively, the new masks could be generated at a different wavelet level. 

Other Processing Modes 

There are many other modes of region processing operation that have not been 
20 presented in this document. The whole user defined region generation and MUX 
modes were not included. The processing modes introduced in this Chapter can be 
used in the same way to control each region channel. Simple primitives are cheap jn 
terms of overhead. One mode of operation supports full mask sets so that arbitrary 
mask can be loaded and sent with the data to form the region channels. The mask 
25 coverage technique can be selected as common masks (a smaller raw size that is up 
sampled at the decoder side) or VM masks, which retain the original image 
dimension. More time is required to study the effects of both techniques as well as 
other region formation schemes. 

Experimental results for the resolution progressive modes of operation are not 
30 presented in this document. These modes of operation are not currently available. 
However their implementation is not difficult. These MUX modes are currently 
under development. 
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The experimental results presented in this Chapter were obtained using DACs ID bit 
level sorting implementation. DAC is currently incorporating the 2D version based 
on EQW into the region processing channel. There may be benefits in retaining both 
techniques. For example, in bi-level image processing. Currently the EQW sorting is 
5 implemented for normal (non-region) processing MUX modes. Both sorting modes 
are available for lossless compression. There are numerous w/avelet kernels and color 
transforms for both lossy and lossless compression. In addition, the number of region 
channels can be set with a current maximum of 4. Primitives can be used as desired. 
The number of primitives is not restricted to 4 with region overlaps given to the most 
10 important region. 

Bit Stream Svntax 

DACs current bit stream structure is rather dynamic given the different types of 
organizational strategies that exist in the underlying core technology. Normal, region, 
15 and mixed processing modes in addition to arbitrary wavelet, color, entropy coding 
various sorting stages in both lossy and losses cases. 
Currently the core architecture can be divided into 2 categories. 

• Lossless compression 

Full data sets lossless color selection, lossless integer lifting scheme wavelet 
20 types, lossless sorting and packing stages, with additional entropy coding can be 
selected for completely lossless and slightly lossy (some of the color transforms 
are slightly lossy) image compression. 

• Lossy compression 

In the lossy case, color selections, more wavelet kernels / lifting schemes, lossy 
25 sorting and packing stages, with variable length coding can be selected 

In both lossy and lossless cases region and mixed modes can be used. Currently 
normal operational modes can be realized in an almost transparent fashion for many 
current compression schemes as well as our own internal EQW / ID sorting. The 
general form for the lossy and lossless header structures is given in the tables in 
30 Figures 52 and 53. 
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The current architecture is very flexible for introducing new technologies. Many of 
the modules are completely interchangeable. The exact bit stream syntax depends 
largely on the mode of operation selected by the user (i.e. regions, no regions, sorting, 
lossy/lossless etc.) Generally normal processing modes can be selected for all 
5 transform levels, or for any number of lower resolution transform levels. Region 
levels can be defined for all transform levels, or the higher resolution levels used in 
combination with the lower resolution normal levels. At some future date, the region 
and the normal levels can be mixed. 

Each unit list organized by the MUX has a header tag. This header tag carries the list 
10 size and the high bit plane processing level for the list. The current implementation 
uses 5 bytes per list. However bit packing is currently under implementation for tag 
headers. This will reduce this by a significant amount. Only 2 packing schemes are 
used at the global level so that no additional tag header information is currently 
required. A diagram illustrating the structure of the tag header is given in Figure 53. 
15 As other core technologies are added to the core engine or the core technology 
advances, other fields may be required in the tag headers. 

The basic structure for normal modes of operation is given in Figure 54. The diagram 
illustrates a lossy pack arrangement for a color image (YUV down sampled color 
transform assumed) according to the packing/processing order given in Figure 30 
20 above. In this case the code stream consists of the file header, followed by the header 
tags/data. 

The basic structure for region processing modes of operation is given in Figure 55. 
The diagram illustrates a lossy pack arrangement for a color image (YUV down 
25 sampled color transform assumed) according to the packing/processing order given in 
Figure 40 above. In this case the code stream consists of the file header, region 
header, regions description and finally the header tags/data. 

The basic structure for mixed processing modes of operation is given in Figure 56. In 
this case the first item in the code stream is the file header followed by the normal 
30 processing mode header tags/data. Next in line is the region header, region 
description followed by the header tags/data. Notice that there is a region start level 
parameter 'r' in this case. Since there is one less transform level for UA^ channels in 
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the lossy case, there is a processing shift in the UA^ channels. This makes reference 
to the "level split" notion discussed in Chapter V whereby the natural 4:1:1 
relationship that exists for inverse color/wavelet transforms can be maintained for 
internal processing and the final code stream. 

Transcodability 

In order to illustrate the generality and openness of RICS, two examples are shown to 
depict how the 'transcode* with bi-level images (JBIG) and JPEG can be handled in 
RJCS. 

JPEG transcode 

Figure 57 shows how a DCT-based code can be produced in RICS. With this 
transcodability, it will be straightforward to write a small conversion program to 
transcode an old JPEG file into a JPEG 2000 file. Although it is possible to transcode 
a JPEG 2000 file into a JPEG file (losing some scalability), it will probably be of very 
little use. 

VJBIG tran.scnde 

Figure 58 shows the execution path illustrating how a bi-level image such as a text 
document can be efficiently handled in the RJCS system. The NULL transform is 
applied (nothing has to be done). Generally binary text documents contain mostly 
high frequency energy. Multi-resolution decompositions will not necessarily be a 
suitable basis for efficient coding. In fact, current well-known techniques specialized 
for binary images are applied directly in spatial domain. While using a RICS system 
for processing a compound document with mixed grayscale/color/text information, a 
set of rectangle shapes suffice to enclose the text regions. This is roughly equivalent 
to run-length coding in existing binary image compression techniques to skip strings 
of zeros. The pixels within each rectangular region are then coded into a one- 
dimensional stream via 1-D algorithms or JBIG routines, as discussed above. 

Post Processin g 

For reconstructed images with wavelet based coding methods at low bit rate, de- 
ringing is usually a useful post processing procedure to improve (mainly) the visual 
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quality. Nonadaptive procedures, which apply the de-ringing fihering to all pixels 
without discrimination, have the following problems. 

• While they can successfully remove ringing artifacts, they may also wash out 
details of image. 

5 • They usually rely on too many parameters that are to be determined by human 

users. 

• The process is time consuming. 

The RICS system employs an adaptive de-ringing algorithm for post processing. 
Since the ringing artifacts usually appear around edges where sharp changes occur, 

10 this adaptive filtering process is applied only to edge areas. As the result, the post 
processing removes artifacts around edges and prevents fine details from being 
smoothed out by the filter. Compared with non-adaptive methods, the processing time 
is remarkably reduced since the fihering is applied selectively to pixels. The diagram 
of Figure 59 shows the algorithm of the adaptive post processing filter. First of all, the 

15 pixels of reconstructed image are classified into edge pixels and non-edge pixels. 
Edge pixels are enlarged into edge regions since the ringing artifacts do not occur 
right at the edges but around the edges. Then the artifact removal filter is applied 
only to the edge areas. 

We did some tests on this adaptive de-ringing procedure. For the filtering stage we 
20 used the filter described in the JPEG 2000 VM 3.0 (B). Figure 60 shows the result of 
the modified post processing. Artifacts are clearly visible in the first picture especially 
behind the cameraman's back and around the tripods. As we can see from the second 
picture, artifacts are removed, however, the details in the grass field and on the man's 
pants are also removed. In the last picture, the modified filter successfully removes 
25 artifacts and at the same time retains the details. Figure 61 shows the edge area .(in 
white color) that was used in producing the third picture. 

As the result of using this adaptive filtering, processing time is greatly reduced since 
less pixel are processed. The table in Figure 62 shows the performance comparison of 
the VM 3.0 (B) post filtering and the RICS adaptive filtering 

30 Edge Detection 
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To determine if a pixel is an edge pixel, the average power difference between the 
pixel and its neighboring pixels are measured against a threshold (bottom threshold). 
If the average power difference of a pixel is above this threshold, the pixel is marked 
as an edge pixel. The threshold acts as a high pass filter that filters out the pixels with 
5 low power levels (Band pass filters can also be used for edge deteaion by introducing 
a proper ceiling threshold). 

Edges are thickened to form edge areas after the edge pixels are deteaed. With the 
right choice of the threshold, this modified filtering process can be applied to 
reconstructed images with compression ratios as low as 8:1. Since for most wavelet 
JO coded images, no artifacts can be visually detected for compression ratio lower than 
8:1, our results suggest that this filter can be applied to most reconstructed pictures 
regardless of its compression ratio. 

Parameters Optimization 

One of the problems with Shen's filter is that there are many parameters associated 
15 with the filter, which makes the filter difficult to use. Reducing the number of 
parameters will improve the usability of the post processing filter. 
There are three different potential ftinctions implemented in VM 3.0 (B) for 
controlling the post processing: Quadratic Truncation, Huber, and Lorenzian. Our test 
results shown that the potential functions Huber and Lorenzian do not outperform the 
20 Quadratic Truncation in any of the 62 cases for PSNR measurements. Furthermore, 
visually Lorenzian creates many noticeable ghost images and shadows, and Huber 
does remove artifacts successfully however it blurs the image more than Quadratic 
Truncation. 

Computation Speed 

For evaluating the three potentiar functions, the following computation needs to be 
calculate/7x(n-l)/2 times for every pixel (where « = filter size x 2- \). 

Lorenzian: 1 float division, 2 float multiplication, 1 float addition and 1 

float 

log operation 
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Huber: 



1 comparison, 1 float multiplication or 2 float multiplication 



with 1 



addition 



Quadratic Truncation: 1 comparison and 1 float multiplication 

5 In our tests, Quadratic Truncation is the fastest and it performs the best for artifact 
removal with the least degradation to image quality. Therefore, the other two 
potential functions will not be used. 

Threshold (y parameter in the potential function) 

The y parameter is examined using the quadratic truncation potential function (default 
10 value is 16), given by: 



We can see from the equations that only pixels with similar neighboring pixels are 
affected when varying y. When y is decreased, there will be more y^ as the result of p 
15 and the intensity of similar pixels output image will be more uniformed resulting a 
more blocky looking type of image (less gradual) for regions, with similar intensities. 
Experiments were performed with 2 images (one grayscale and one color image). The 
table in Figure 63 shows the result for the color image. 

We can see from the results that the higher the y, the worst the image quality. 

20 However, the difference in image quality is not very significant. The pictures look 
very similar and there are no significant changes that can be detected. However as 
predicted, the blocky looking clouds (see Figure 60) are seen when y = 8 or below. 
Therefore by setting y lower we can increase the image quality by some small 
amounts without changing the visual quality of image. Therefore y = 12 was chosen to 

25 be the default value for y parameter in the potential function. 



e = arg min ]^ p{X, - X j) 



where, p = { 
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Filter Lenp th 

Filter length (F) determines of the size of samples collected for pixel estimations. It is 
obvious that the larger the filter length the longer the processing time. The default 
length given by VM3.0 (B) is 9 pixels. Experiments are performed with varying filter 
lengths to examine how the length affects the image quality (MSE and PSNR) and the 
artifact removal ability. Two pictures are used in this experiment (one color and one 
grayscale). The results are shown in the following tables in Figure 64 and 65. 

As we can see from the experimental results, the small the length, the better the image 
quality and the shorter the processing time. However, artifacts are not completely 
removed for length = 5 (or smaller) at high compression ratios. The default value for 
the filter length is chose to be 7 pixels to increase image quality and to decrease 
processing time without degrading the ability for artifact removal. 

Constraint 

The value of constraint (ThI) affects the fiker as described in the equations below. 
y = x+c{d ,ThI) 
d= x-x 

c{d,Thl) = sign(aO * Max(0, ahs{d) - Max(0, 2*(ahs{d) - ThJ))) 
From the equations above, we can see that as Thl increases, more pixels will have 'd' 
as the output for c{d.ThJ). Therefore, as Thl increases, more estimate. value of x will 
be used as the final estimate and the image quality will decrease while the smoothness 
of an image will increase. In order to keep the image quality as high as possible, Thl 
should be kept as low as possible. However, artifacts might not be properly removed 
if Thl is set too low. Three images were tested with different level of constraints in 
the experiment. The test results are given in the tables in Figures 66 to 68. 

These results show that image quality decreases as constraint increases. However, 
some artifacts (very small, can only be seen when the picture is enlarged) are not 
completely removed when C is equal to 4 or less. In order to achieve the best image 
quality while successfully remove artifacts, the value 8 is suggested to be the default 
value. 
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Iteration 

It is obvious that processing time is directly proportional to the number of iterations. 
However, artifacts might not be removed successfully if the image is not processed 
enough times. Figure 69 clearly shows that image quality decreases as the number of 
5 iteration increases. However, when image is enlarged, small artifacts can be still 
visible for Rl. In order to achieve the best image quality while successfully remove 
artifacts, 2 iterations is suggested to be the default setting. 

Mask Shape 

Different shapes of masks are discussed in Shen's paper, however, only the one with 
10 the shape of a + sign was implemented in VM3.0 (B). Different masks were 
experimented and the shape of mask is best left unchanged since simplifying the mask 
degrades the performance of the filter and complicating the mask increases processing 
time greatly. 

Using the Modified Post Processing Filter 

15 The adaptive de-ringing is implemented based on the post processing filter in VM3.0 
(B). Few parameters have been eliminated and default values have been established to 
increase the usability of the post processing filter. 
To use the modified post processing filter: 

Usage: post2 -s width height bpp -i infile -o outfile [-1 mask lower threshold] 
20 [-W mask width] |-t thresh] [-f fjength] (-c constraint] [-r iterations] 

Default values will be used for the parameters if the parameters are not specified in 
the command line. The default parameters are: 

mask lower threshold: lower threshold for edge detection [30] 
mask width: width of mask [12] 
25 thresh: y parameter in the potential function (for estimation) 1 1 2) 

fjength: filter length J7) 

constraints: constraints using in the clipping fiinction [8] 
iterations: number of iterations (2] 

The modified post processing filter performs well with the above default parameters, 

30 however, these parameters can be changed at the command line if the user washes. 

The post processing filter can" now be applied to pictures with 8:1 compression ratio 
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with very small increase in MSE. The modified post processing filter seems to 
improve image quality for compression ratio beyond 11:1. 

In practice, there are two ways for the user to play with the settings: on the encoder 
side and on the decoder side. 

5 Decoder end 

The decoder will have the fiill control of the post processing and the of the associate 
parameters. The encoder will have no control on how the images will look at the 
decoder end. Accordingly, there will be no addition to file header since no post 
processing parameters will be included. 

10 Encoder end 

When necessary, the encoder can also predetermine the post processing filter 
parameters. Post processing filter parameters will be stored in the image header, and 
the image will be restored according to these parameters. In this setting, the user at 
the encoding end knows exactly how the image will look at the decoder end. 
15 However, adding these parameters in the header will increase the size of the 
compressed image. 

Total of 6 parameters will need to be packed in the header: mask lower threshold, 
mask width, estimation threshold, filter length, constraints and number of iterations. 
The table in Figure 70suggests range limits on the parameters to reduce the header 
20 size. 

In total, three bytes will be needed to include these parameters in the header. 
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WE CLAIM 

1 . A region-based method for encoding and decoding digital still images to 
produce a scalable, content accessible compressed bit stream comprising the 
steps: 

decomposing and ordering the raw image data into a hierarchy of multi- 
resolution sub-images; 

determining regions of interest; 

defining a region mask to identify regions of interest; 

encoding region masks for regions of interest 

determining region masks for subsequent levels of resolution; and ' 

scanning and progressively sorting the region data on the basis of the 
magnitude of the multi-resolution coefficients. 

2. An apparatus for the region-based encoding and decoding of digital still 
images that produces a scalable, content accessible compressed bit stream 
comprising: 

a means of decomposing and ordering the raw image data into a hierarchy of 
muhi-resolution sub-images; 

means of determining regions of interest; 

means of defining a region mask to identify regions of interest; 

means of encoding region masks for regions of interest; 

means of determining region masks for subsequent levels of resolution; and 

a means for scanning and progressively sorting the region data on the basis of 
the magnitude of the multi-resolution coefficients. 

3. A region-based system for encoding and decoding digital still images that 
produces a scalable, content accessible compressed bit stream and comprises 
the steps: 

decomposing and ordering the raw image data into a hierarchy of multi- 
resolution sub-images; 

determining regions of interest; 

defining a region mask to identify regions of interest; 
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encoding region masks for regions of interest 

determining region masks for subsequent levels of resolution; and 

scanning and progressively soning the region data on the basis of the 
magnitude of the multi-resolution coefficients. 

A method for encoding and decoding digital still images to produce a scalable 
content accessible compressed bit stream comprising the steps: 

decomposing and ordering the raw image data into a hierarchy of multi- 
resolution sub-images; 

setting an initial threshold of significance and creating a significance index; 
determining an initial list of insignificant blocks; 

forming the list of significant coefficients by encoding a significant map usine 
a quadtree representation; 

recursively reducing the threshold values and repeating the encoding process 
for each threshold value; and 

transmitting refinement bits of significant coefficients. 

An apparatus for encoding and decoding of digital still images that produces a 
scalable, content accessible compressed bit stream comprising: 

a means of decomposing and ordering the raw image data into a hierarchy of 
multi-resolution sub-images; 

means for setting an initial threshold of significance and creating a 
significance index; 

means for determining an initial list of insignificant blocks; 

means of forming the list of significant coefficients by encoding a significant 
map using a quadtree representation; 

a means of recursively reducing the threshold values and repeating the 
encoding process; and 

transmitting refinement bits of significant coefficients. 

A method of decoding digital still images to produce a scalable, content 
accessible compressed bit stream comprising the steps: 

decoding the bitstream header; 
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determining the initial threshold values and the array of initial significant 
pixels, insignificant bits and wavelet coefficients; 

decoding the significance maps; 

modifying the significance lists and decoding the refinement bits for each 
threshold level; 

reconstruct the wavelet coefificient array; 

perform the inverse wavelet transform; and 

reconstruct the image. 

7. A method of transmission of digital signals that creates a scalable, content 
accessible bitstream comprising the steps; 

pack the most significant bits of the largest coefficients first followed by 
refinement bits and the most significant bits that are significant for coefficients 
at the next bit level; 

repeat this process in a recursive fashion until the desired compression size is 
obtained; 

calculate the pack ratios to be used for each channel of the wavelet 
decomposition hierarchy by taking the ratio of the two largest amounts of data 
to the smallest amount of data; 

determine the optimal amount of data to allocate for each color channel based 
on the user specified compressed file size; and 

if performing region of interest processing, consider the packing overhead 
introduced by the mask when determining the bit budget for each channel. 
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Magnitude vs. Diagonal Position 
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Y-chaiuid: LL,-HL,-LI1,-HH,-HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH, 
U-chaiuiel: LLj- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LHi-HH, 
V-chaiuiel: LL3- HL3-LH3-HH3-HL2-LH2-HH2-HL,-LH,-HH, 
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Yll4-Ull4-Vu^-Yhl<-Ylh4-Yhh<-Uhl<-Ulh4-Uhh4-Vhi^-Vlh4-Vhh<- 

YhL3-YlH3-YhH4-UhL4-UlH4-UhH4-VhL3-VlH3-VhH3- 
YhL2-YlH2-YhH2-UhL2-Ulh2-UhH2-VhU-VlH2-VhH2- 
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YLW-ULL4-VLU-Ym.,- UHL--Vm^-YLH-.-ULH--VLH-.-YHH4-UHH4-V„H4 

Yhl3-Uhu-Vhu-Ylh3-Ulh3-Vlh3-Yhh3-Uhh3-Yhh3 

YhL2-UhL3-VhL2-YlH2-UlH2-VlH2-YhH2-UhH2-VhH2 

Yhli-Uhli -Vhli -Ylhi-Ulhi-Vlh) -Yhhi -Uhhi-Vhhi 



FIGURE 33 



Data Type: MUXLIST 



Parameters: 



liTotBytesPacked - long integer total bytes packed into the data buffer for this list 
cScheme - character processing scheme used for data contained in this list 

cHighBit - character highest bit-level where data processing begins for this list 

•pucMuxBuff - pointer to unsigned character buffer where data for each 

bit-level is packed for this list 

Fields for MUX: information for packing after list processing is complete. 

pliBitPacklnfol 1 6] - pointer to long integer number of bits packed into the data buffer 

at each bit-level for this list 
liCurBytesCount - long integer current byte count used for bit budget distribution 
when packing this list 

for packing this list 



cCurBitLevel 
cRemainingBits 



- character remaining bits to be packed at a given bit-level when 
data to be packed is not evenly divisible by 8 for this list 
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CALCULATE Channel BitBudget // determine optimal bit-budget for each color channel. 
INITIALIZE Channel CuirentBitPlane // highest bit plane that exists in each color channel. 
INITIALIZE liCurBytesCount and cRemainingBits FOR each MUX list 

FOR Each Color Channel // process each channel separately. 

WHILE Channel BitBudget > 0 AND Channel CurrentBitPlane >=0 

FOR Each Wavelet Transform level // beginning at lowest resolution level. 

FOR Each Orientation Set of Data // according to lossy case natural processing order. 
IF cCurBitLevel NOT EQUAL to Channel CurrentBitPlane 

CONTINUE 
ELSE 

SET BitLevelBytes to pliBitPacklnfolCtaama CurrentBitPlane] » 3 

SET RemBits to pliBitPadcInfo\Claimcl CurrentBitPlane] & 7 

IF S\im{cRemainingBits, RemBits) >= 8 
INCREMEhJT BitLevelBytes by 1 
DECREMENT cRemainingBits by 8 - RemBits 

ELSE 

INCREMENT cRemainingBits by RemBits 
ENDIF 

IF Channel BitBudget >= BitLevelBytes 

INCREMENT liCurBytesCount by BitLevelBytes 
DECREMENT Channel BitBudget by BitLevelBytes 
DECREMENT cCurBitLevel by 1 

ELSE 

INCREMENT liCurBytesCount by Channel BitBudget 
SET Chaiuiel BitBudget to 0 
ENDIF 
ENDIF 

IF Channel BitBudget EQUALS 0 

BREAK 
ENDIF 
END FOR 
END FOR 

DECREMENT Channel CurBitPlane by 1 
END WHILE 
END FOR 
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Y-channel: R,(LL4-HL,-LH4-HR,-HLj-LH3-HH3-HLj-LH2-HH2- HL,-LHrHH,), 
Rj(LL,-HL.-LH4-HH4-HL3-LH3-HH3-HL2-LH2-HHrHL,-LH,-HHi). 
R3(LL,-HL4-LH,-HH4-HL3-LH3-HH3-HLj-LH2-HHj- HL,-LH,-HH,), 
R,(LL4-HU-LH4-HH.-HL3-LH3-HH3-HL2-LH2-HH2-HL,-LH,-HH,). 

U-channel: R,(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,). 

R2(LL3- HLa-LHa-HHj-HLj-LHj-HHj- HL,-LH,.HH,). 
RjOIo- HL3-LH,-HH3-HL2-LH2-HH2- HL, -LH, -HH,). 
R4(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LHrHH,). 

V-channel: RiO-Lj- HL3-LH3-HH3-HL2-LH2-HH2- HL,-LH,-HH,), 
R2(LL3- HL3-LH3-HH3-HL2-LH2-HH2- HL,.LH,-HH,), 
RjOXj- HL3-LH3-HH3-HL2-LHJ-HH2- HL,-LH,.HH,). 
RiCLLj- HL3-LH3-HH3-HL2-LHJ-HH2- HLi-LHi-HHi). 
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Rl (VLL<-YHL<-yLH4-YHH4-Uu3-Vuj- 

Yhu-YlH3-YhH3-UhL3-UlH3-UhH3-Vhl3-Vlh3-VhH3- 

Yhu-Yuo-Yhh2-Uhl2-Uue-Uhh2-Vhl2-Vlh2-Vhh2- 
Ym.i-YLH.-YHHrUHi.,-ULHi-U™,-V„Li-VLHrVm,). 
RjCYlw-Yhu-Yum-Yhh^-Ulu-Vuj- 

Yhu-Yuo-Yhh3-Uhl3-Ulh3-Uhh3-Vhl3-Vlh3-Vhh3- 

Ynu*YLH2"YHH2"UHL2"UtJi2-UnH2-VHL2-VLH2-VHH2* 

YjjLi -Ylhi -Yhhi -Uhli -Ulhi -Uhhi - Vhli • Vlhi - Vhhi ), 
R3(Yll<-Yhl4-Ylh4-Yhh4-Ulu-Vlu- 

YhL3-Yud-YhH3-UhIJ-UlH3-UhH3-VhL3-Vuo-VhH3- 

YHu-Yun-YHHj-UHLz-Uuc-UHHj-VHu-Vuo-Vinn- 
Yhli -Ylhi -Yhhi -Uhli -Ulhi -Uhhi -Vhli -Vlhi-Vhhi ). 
R<(Yll4-Yhl4-Ylh4-Yhh4-Ulu-Vll3- 

Yhl3-Ylh3-YhH3-Uhl3-Ulid-UhH3-Vhl3-VlH3-VhH3- 
YhL2-YlH2-YhH2-UhL2-UlH2-UhH2-Vhu-VlH2-VhH2- 

Yhli -Ylhi -Yhhi -Uhli -Ulhi -Uhhi -Vhli -Vlhi-Vhhi)- 
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Ri (Ylu-Yhli-YlwYhh^-Uuj- Vuj- 

Yhu-Uhij-Vhu-Yu13-Uuu-VlH3-YhH3-UhH3-Vhh3- 
Yhu-UhL2-Vhl5-YlH2-UlH2-Vuc-YhH3-UhH2-VhH2- 

YHLi-U„LrV„LrYLH.-ULH.-Vu„-Yra,-Um,-Vm,), 

R2(Yll4-Yhu-Ylh4-YhH4-Uuj-VlL3- 

Ymj-UhL3-Vhu-Ylh3-Uud-Vuo-YhH3-UhH3-Vhh3- 
YMJ-UHL3-VKj-YLm-ULH2-VLH2-YHH2-UHH2-VHH2- 

Yhli -Uhli - Vm., -Ylh, -Ulhi - Vm, -Yh„, -Uhhi - Vhhi ). 
R3(Yli^-Yhl4-Ylh4-Yhh<-Uuj-Vuj- 

Yhu-Uhl3-Vhu-Yud-Ulh3-Vuo-Yhh3-Uhh3-Vhio- 

YhU-UhL2-VhL2-YlH2-UlH2-Vijc-YhH2-UhH2-VhH2- 

Yhli -Uhli - Vhli -Ylhi -Uuh -Vu,, -Yhhi -Uhhi -Vhhi ), 
R«(Yll4-Yhm-Ylh4-Yhh<-Ulu-Vuj- 

YhL3-Uhu-Vhu-YlH3-Uud-Vud-YhH3-UhH3-VhH3- 
YHL2-UHL2-V„U-Yuo-ULH2-Vuo-YHm-UHH2-VHH2- 

Yhli-Uhli-Vhli-Ylhi-Ulhi-Vlhi-Yhhi-Uhhi-Vhhi). 
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R 1 (Y lw-Ull4- Vu.,-Yhl-.-Ylh4-Yhh-,-Uhl4-Ulh4-Uhh4-Vhl<-Vlh4- Vhh4- 

YhL3-Ylh3-YhH4-Uhl4-UlH4-UhH4-Vhl3-Vuc-VhH3- 
YhL2-YlH2-Yhh2-Uhu-UlH2-UhH2-VhL2-Vuc-VhH2- 

Yhli -Ylhi -Yhhi -Uhli -Uuji -Uhhi -Vhli -Vlhi -Vhhi ). 

R2(YLL4-ULU-VLL4-YM^-YLH4-YHH4-UHL4-ULH4-UHH4-VHL4-VLH4-Vmi4- 
YhL3-Ylh3-YhH4-UhL4-UlH4-UhH4-VhL3-VlH3-VhH3- 
YhlS-YlH2-YhH2-UhL2-UlH2-Uhh2-Vhu-Vuc-VhH2- 

Yhli -Ylhi - Yhhi -Uhli -Ulhi -Uhhi -Vhli -Vlhi - Vhhi ). 
R3{Yll4-Ull4-Vll4-Yhl4-Yl„4-Yhh4-Uhl<-Ulh4-Uhh4'Vhl4-Vlh4-Vhh4- 

YhU-YlH3-Yhh4-UhL4-UlH4-UhH4-Vhl3-VlH3-VhH3- 
YhL2-YlH2-Yhh2-UhL2-Ulh2-UhH2-VhL2-VlH2-VhH2- 

Yhli t Ylhi -Yhhi -Uhli -Ulhi -Uhhi -Vhli -Vlhi -Vhhi ), 
R4(Yll4.Ull4-Vll4-Yhl,-Ylh4-Yhh4-Uhl4-Ulh4-Uhh4-Vhl4-Vuj4-Vhh4- 

YhL3-Ylh3-Yhh4-UhL4-UlH4-UhH4-VhL3-Vuo-Vhh3- 

Yhl2-Ylh2-Yhh2-Uhu-Ulh2-Uhh2-Vhl2-Vlh2-Vhh2- 
Yhli -Ylhi -Yhhi -Uhli -Ulhi -U„„, - Vhl, -Vl„, - Vh„, ). 
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Yu.4-Yhl4-Ylh4-Yhh4-Ull3-Vijj- 

YhU-YlH3-YhH3-Uhl3-UlH3-UhH3-Vhu-VlH3-Vhh3- 
YhLJ-YlH2-Yhh2-UhL2-UlH2-Uhh2-Vhl2-VlH2-VhH2- 

Yhli-Ylhi-Yhhi-Uhli-Ulh,-Umi,-Vhli-Vuii-Vhhi (R1, R2, R3 , R4). 
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YLW-ULU4-VLU-Y„,^-YuH4-Ym4-UHL^-ULH4-UHH.-V„L4-VL„.-V„„4- 

YhL3-YlH3-YhH4-Uh1/4-UlH4-UhH4-VhL3-VlH3-VhH3- 

Yhl2-YlH2-Yhh2-Uhu-UuH2-UhH2-VhL2-VlH2-VhH2- 

YHLrYLHrYHH,-U„LrULHrUHHi-VHLi-VLHrVHHi (Rl. R2, R3. R4) 
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MSE vs. Bpp for Glacier Park Image 
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MSE vs. Bpp for Glacier Park Image 
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MSE vs. Bpp for Glacier Park Image 
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MSE vs. Bpp for Glacier Park Image 
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Ima^e Width 


10 


(64-1024 columns) 
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Image Height 


10 
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Wavelet Transfomi Levels 
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Level >= 1 


Region Channels 


2 


1.2.3,4 


Mask Type 


2 


None, User, Auto (DCT) 


Pack Raw Mask Flag 


1 


Read in file for mask 


Mask Procedure 


3 


DCT Common Mask, Raw Common Mask with VM Translation, 


Down Sampling Type 


2 


Heuristic DCT Down Sampling 


Region Start Level 


4 


No Regions to H wavelet levels 


Lossy Flag 


1 


Lossless 


Color Transform Type 


2 


2 Internal, YIQ (full data sets) 


Sort Type 


1 


1D,E0W 


Wavelet Kernel 


2 


Lifting Scheme 
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Compressioo 


VM 3.0 (B) 
time (sec) 


RICS 
Adaptive 
time (secj) 


Compression 
ratio 


VM 3.0 (B) 
time (sec) 


RICS 
Adaptive 
time (sec) 


3.2170 


7.310 


1.743 


3.709 


22.392 




9.0140 


7.260 


1.702 


8.889 


23.023 


8.151 


40.210 


7.325 


1.482 


10.000 


22.442 


5.288 


48.763 


7.170 


1.462 


11.428 


22.232 


5.809 


56.492 


8.262 


1.402 


13.333 


22.462 


5.878 


66.550 


8.703 


1.382 


16.000 


22.482 


5.629 


81.143 


6.769 


1.362 


20.001 


22.522 


5.489 


103.716 


7.171 


1.175 


26.666 


22.382 


5.568 


130.168 


7.170 


1.222 


32.000 


22.422 


5.338 


169.019 


8.310 


1.132 


40.002 


22.352 


3.896 


253.337 


7.130 


1.052 


79.987 


22.282 


3.465 
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Compression 
Ratio 


PSNR (Y=8) 


PSNR(r=12) 


PSNR (7=16) 


3.709 


31.329611 


31.177579 


30.511348 


4.000 


31.300452 


31.154471 


30.490295 


5.333 


31.157031 


30.994670 


30.357978 


8.000 


30.630383 


30.474574 


29.862022 


8.889 


30.401649 


30.216825 


29.646734 


10.000 


30.067778 


29.909074 


29.375402 


11.428 


29.691347 


29.546710 


29.056083 


13.333 


29.228454 


29.103298 


28.645196 


16.000 


28.407724 


28.276664 


27.908937 


20.001 


27.508670 


27.413509 


27.134810 


26.666 


26.385419 


26.307257 


26.084268 


32.000 


25.333808 


25.280108 


25.111808 


40.002 


24.696655 


24.660851 


24.509536 


79.987 


22.179274 


22.152787 


22.093103 


100.004 


21.592595 


21.568550 


21.525140 



FIGURE 63 



31/35 



Copied from 1034©siigrpTfiiS(Bi4?t(F^£eQS) 



wo 00/49571 



PCT/CAOO/00134 



Compression 
Ratio 


PSNR (F=3) 


PSNR (F=5) 


PSNR (F=7) 


PSNR (F=9) 


PSNR(F=11) 


1.583 


35.514592 


34.206867 


33.687448 


33.409849 


33.223806 


].600 


35.516204 


34.201801 


33.674565 


33.408178 


33.220934 


2.000 


35.314595 


34.080237 


33.538998 


33.301232 


33.090358 


2.667 


34.957612 


33.765719 


33.291270 


33.011713 


32.832969 


4.000 


34.002797 


32.997735 


32.564248 


32.338257 


32.200522 


8.000 


31.089735 


30.554027 


30.324381 


30.196463 


30.106473 


8.889 


30.572145 


30.142348 


29.935209 


29.796335 


29.725453 


9.999 


30.131952 


29.795631 


29.604549 


29.472401 


29.394763 


11.429 


29.629346 


29.390456 


29.261471 


29.13131 


29.064764 


13.334 


28.854873 


28.701271 


28.602145 


28.527964 


28.470308 


16.000 


27.870101 


27.760478 


27.671142 


27.632609 


27.57%55 


19.999 


27.031806 


26.910914 


26.857123 


26.835492 


26.794267 


26.662 


25.962692 


25.994100 


25.999759 


26.019787 


25.991909 


40.010 


24.159134 


24.178761 


24.173011 


24.191618 


24.182049 


80.020 


21.654340 


21.683920 


21.731559 


21.761676 


21.774576 
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PSNR (F=5) 


PSNR (F=7) 


PSNR (F=9) 


PSNR (F=ll) 


3.709 


31.656960 


30.740420 


30.806267 


30.511348 


30.298514 


4.000 


31.634763 


30.735842 


30.791791 


30.490295 


30.278695 


5.333 


31.436411 


30.558898 


30.635121 


30.357978 


30.141580 


8.000 


30.819587 


30.065789 


30.126713 


29.862022 


29.680912 


8.889 


30.537579 


29.839074 


29.879854 


29.646734 


29.474362 
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30.170406 


29.541099 


29.603527 


29.375402 


29.205815 


11.428 


29.740858 


29.186726 


29.255123 


29.056083 


28.880456 


13.333 


29.255225 


28.752598 


28.839294 


28.645196 


28.494299 


16.000 


28.391715 


27.982943 


28.075866 


27.908937 


27.788081 


20.001 


27.479117 


27.162369 


27.270921 


27.134810 


27.012976 


26.666 


26.329007 


26.079132 


26.175590 


26.084268 


25.999266 


32.000 


25.279960 


25.088255 


25.181387 


25.111808 


25.044080 


40.002 


24.649334 


24.486896 


24.568579 


24.509536 


24.446939 


79.987 


22.195351 


22.099853 


22.127757 


22.093103 


22.066329 


100.004 


21.598939 


21.520845 


21.561808 


21.525140 


21.501686 
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1.583 


29,654831. 


33.409849 


35.271851- 
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